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1. Introduction. The theory of the motion of a projectile is usually based on the 
assumption that the projectile possesses complete rotational symmetry, in the sense 
that rotation through any angle about the axis of symmetry leaves the projectile as a 
whole unchanged in outward form and in internal mass distribution. This assumption 
of symmetry is suitable in the case of an ordinary shell, but is not suitable in the case 
of a finned projectile which possesses symmetry of a more restricted type. The purpose 
of the present paper is to investigate the effect of various types of rotational and re- 
flectional symmetry on the aerodynamic force system acting on a projectile. According 
to the basic aerodynamic hypothesis’ the aerodynamic force system depends on the 
instantaneous motion of the projectile. The form of this functional dependence is re- 
stricted by the symmetry of the projectile. Our purpose is to investigate this restriction 
mathematically. 

The instantaneous motion of a projectile is described by choosing some base-point O, 
fixed in the projectile, and specifying the velocity (u) of O and the angular velocity 
(w) of the projectile. The aerodynamic force system, exerted by the air on the projectile, 
is described by the equipollent system consisting of a force (F) acting at O, together 
with a couple (G). Then the vectors F and G are functions of the vectors u and w. They 
are also functions of the air density (p) and the local sound velocity (c), but these func- 
tional dependences will be suppressed as a matter of notational convenience. Choosing 
any rectangular coordinate axes Ox,x.x; , we may express the functional dependence of 
force system on motion by the six equations 

F,, = fm(Ur » Uz» Us » Wi y We , Ws), 


(1.1) 


Gn _ Jm( Uy » U2 ’ Uz , @ » We; ws), (m - 1, 2, 3), 


Um, Om, Fm, Gm being components along Ox,, . 

2. Types of symmetry. Let us consider projectiles possessing one or both of the 
following types of symmetry (Fig. 1): 

(¢) n-gonal rotational symmetry about an axis A; 

(72) reflectional symmetry in a plane P. 
Symmetries may be described in terms of “covering operations’ which leave the pro- 
jectile as a whole unchanged. We shall be concerned only with the aerodynamic force 


1Received March 15, 1948. The theory contained in this paper was developed under contract with 
Aberdeen Proving Ground, and is now published with permission. 
2K. L. Nielsen and J. L. Synge, Q. Appl. Math. 4, 201-226 (1946) 
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system and not with the dynamics of the projectile. This means that we are interested 
in its external form, but not in its distribution of mass. Consequently all remarks about 
symmetry are to be interpreted relative to the external form only. The covering opera- 
tion for n-gonal rotational symmetry about an axis A is a rotation about A through 
an angle a where a = 27/n, n being an integer greater than unity. (An ordinary bomb 
with four fins possesses tetragonal symmetry, n = 4.) For certain formal purposes, it 
is convenient to include the case n = 1, although this actually implies no symmetry. 
The covering operation for reflectional symmetry in a plane P is a reflection in P in 
the ordinary geometrical sense, i.e. each particle of the projectile is replaced by another 
particle at a position such that the line joining the two positions is bisected orthogonally 
by P. The covering operation for n-gonal symmetry may be carried out physically; 
that for reflectional symmetry cannot. 

However, we may view the question of symmetry in a different, but equivalent, 
way. We may apply the covering operations, not to the projectile, but to a set of axes. 
Consider first the case of n-gonal symmetry about an axis A. Choose a set of rectangular 
axes Ox,x.2; with origin on A. Apply the covering operation to these axes, leaving the 
body unmoved. Let the new axes be Oz{xi13 . The two sets of axes are equivalent in the 
following sense: Any general statement about the connection between aerodynamic force 
system and motion must have the same form no matter which axes are used. Otherwise 
put, it is impossible for an experimenter to determine, by aerodynamic measurements, 
which of the two equivalent sets of axes he is using. This means that the six aerodynamic 
functions in (1.1) are the same for two sets of equivalent axes. If primes indicate com- 
ponents along Ox{xr3x3 , we must have 

Fre = fm(ui , U2, Us , WI , w2 , Ww), 
(2.1) 
Gre = Gm(Ui , US , UZ , Wi , W2 , 5), (m = 1, 2, 3), 


J, and g,, being the same functions in (1.1) and (2.1). 

The same remarks apply in the case of reflectional symmetry in a plane P. In this 
case, we choose the origin O on P, and reflect Ox,x.x; in P to obtain Oxjxjx3 . However 
some care must be exercised in this case because the orientation of the axes is changed. 
This introduces no difficulty in the case of the polar vectors u and F, but must be taken 
into consideration in the case of the axial vectors w and G. This point will be discussed 
below in the appropriate place. 

On account of the ease with which the notation of complex variable lends itself to 
the treatment of rotations in a plane, we find it convenient to use the following notation: 


U= Uy + ue, ® = w, + lw, F=F,+7F,, G =G, + 74., (2.2) 


with similar definitions for the primed quantities. Let us rewrite (1.1) and (2.1) in the 


form 


F _ flu, u, Us » W; @, Ws), F; ica fa(u, U, Us » @, @, Ws), 
(2.3) 
G = glu, U, U3 7 W, @, Ws), G; = gs(u, U, Uz 7 W, @, @3), 
F’ a fw’, u’, uz ’ w’, a’, w3), F; — fs(u; u’, U3 ’ w’, a’, w3) 
(2.4) 


_ = a? — —p of ‘= 
G’ = glu »U, U3, w’, a’, w5), G3 gs(u’, U,U3,0, a’, 5), 
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where the bars denote complex conjugates. Here the functions f and g are complex 
functions of complex arguments, and f; and g; are real functions of complex arguments. 
Since we have changed the arguments in passing from (1.1) and (2.1) to (2.3) and (2.4), 
the functional forms f;, and g; are different in the two sets of equations. No confusion 
will arise from this, because we shall not refer to (1.1) or (2.1) again. 

3. Consequences of n-gonal rotational symmetry. We suppose that the projectile 
has n-gonal rotational symmetry about an axis A. We choose the origin O on A, the 
axis Ox, along A, and the other two axes in the plane perpendicular to A, but otherwise 
arbitrary except for the condition that they must of course be at right angles. We consider 
some state of motion and the corresponding aerodynamic force system. These are 
described respectively by the complex and real components (u, us , w, w3) and (F, F; , 
G, Gs) on the axes Ox,2.x2, . Now apply to the axes the covering operation, i.e. a rotation 
about Ox, through an angle a = 2x/n. The same motion and the corresponding force 
system may also be described by the complex and real components (u’, us , w’, 3), 

F’, FS, G’, GD on the new axes Ox{xixj . The transformation from one set of com- 
ponents to the other is easily seen to be as follows: 
u’ = we”. Us = Us, wo’ = we, @, = Os, 

(3.1) 
F’ = Fe"*, F; = F;, 7’ = Ge"*, G3; =G;. 


3 


If we substitute these in (2.4) and then compare with (2.3), we get the following equa- 
tions, which must be satisfied for arbitrary values of (u, us , w, ws): 


e'*f(ue**, We'*, us , we, we *, ws) = fu, U, Us , w, @, ws), 
ia ( -ia — ia ta — sa = / — = 
e** g(ue , Ue, Us , We , we", ws) = glU, U, Uz , W, W, Ws), 
(3.2) 
fs(ue**, Ve'*, us , we **, we'*, ws) = fa(u, U, Us , w, @, ws), 
g,(ue~, We'~, Us , we", WE~, ws) = gs(U, U, Us , wW, @, ws). 


We now make an important assumption, namely, that the functions f, fs , g, gj may 
be expanded in power series in the transverse components U, , Us , @ , We , OF, equiva- 
lently, in power series in the complex quantities u, U, w, . We shall write these power 


series in the form 


fF = f(u, U, Us , w, @, ws) = > Pigrs(Ugz , 3)WUw'@’, 
pare 
G = gu, U, Us , w, @, w3) = + Grors(Us » W3)WUw'@', 
pare 
(3.3) 
al — _— (3 — — 
F, = fa(u, U, Us , w, @, ws) = D> FS? (us , ww’, 
pars 
-y — — (3 —@ 1 
G; = g3(U, U, Us , W, @, ws) = > Grors(Us , w3)UU' wo’, 
pers 
the summations covering the ranges p, g, r, s = 0, 1, 2, --- . Here the functions 


Fyers(Us ] Ws), Gyors(Us , 3), FS?) (us ’ Ws), GS) (Us ’ ws) 


are the same functions for equivalent axes. 
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We note that the coefficients F;;,, and G,,,, are in general complex, and certain 
relations must be satisfied by them in order that f,; and g; may be real. We have f; = f; , 
and so 

3) —q T—s8 (3) _ —r is 
Dd F2. (us , www" = DY FE (us , ws)uo'w 
sare ere 
(3.4) 


bs Fe. (Us ’ w3)UW'U'w'@’. 

pars 

This holds for arbitrary values of u,; , Uz , w; , w. , and may be treated formally as an 
identity in u, U, w, w. A similar identity arises from g; , and so we have 

r1(3) wH(3) ¥(3) 7V(3) 

Frare(Us , Ws) = Fi52-(us , ws), Goare(Us » Ws) = Giyar(Us , Ws), 


(3.5) 

(p, q,7,8 = 0, 1, 2, ++), 
We see that F;,;,(us , w3) and G,);,(us , w3) are real; this result is important later in the 
paper. (The summation convention for repeated suffixes does not operate, here or else- 


where. ) 
We now substitute the series (3.3) in (3.2), and obtain the following result from the 


first of (3.2): 
. x Prars(Us , @3)WU'w' oD = Zz. exp [ia(l —p+q—rt+8)|Foo.(Us , w3)UU'w'’. (3.6) 
pars pars 


This may be treated formally as an identity in u, U, w, a, and we obtain 


er coe 7 
F jors(Ug , @3) SIN gal —p+¢e—-7r +4) = @, (p, g, 7, 8 = 0,1, 2, ---). (8.7) 


Thus all coefficients F,,,,,(u; , w3) must vanish except those for which p, q, r, s are such 
that 


ce 
sin5a(l1-—p+q—r+s)=0. (3.8) 


Since a = 27/n, we see that the only surviving coefficients F,4,,(U3 , 3) tn the case of n-gonal 
rotational symmetry are those for which p, q, 7, 8 satisfy one of the equations 


p-qtr-s—-l=mn (m = 0, +1, +2, -:-) (3.9) 


(We recall that p, g, r, s are of course non-negative integers.) Similarly, the only surviving 
coefficients Gyor,(Us , ws) in the case of n-gonal rotational symmetry are those for which p, 
q, T, 8 satisfy one of the equations (3.9). 

A similar argument applied to the third and fourth equations of (3.2) shows that 
the only surviving coefficients F'53),(us , ws) and G,2;,(Us , 3) in the case of n-gonal rotational 
symmetry are those for which p, q, r, 8 satisfy one of the equations 


Pp—-qtr—s=mn (m = 0, +1, +2, ---). (3.10 
The results are most easily interpreted if we introduce P and Q defined by 


P=prtetr, Q=aqts. (3.11) 
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We note that as p, q, r, s cover the range 0, 1, 2, --- , P and Q cover the same range. 
For a given value of P, there will occur in the power series (3.3) the combinations 
> > > 2 Deas P 
ur ” ee ue, tee ia wo, (3.12) 


and for a given value of Q the combinations 


—, -_ ao fas = —2—2 — _ —_ 
uw. u°~"S, u?-“e’, tee ww", a”. (3.13) 


The result regarding surviving terms may be displayed as follows: 


7 (3) 
Fyare(Us ’ ws) I pere(tl3 ] ws) 
ry(3) 
Goers(Us ’ ws) Tpqre(tls ? ws) 
survive only if survive only if 
(3.14) 
P—-Q=mnr+1 P—Q= mn 
m= 0, +1, +2,::: 
(P=ptr, Q=q+:) 
Since, for small cross velocity and cross spin, the significance of the terms in the 


series (3.3) decreases with increasing degree, it appears advisable to arrange the terms 
in order of increasing degree. The degree of a term with coefficient having subscripts 


P, q; r, SIS 
D=pt+qtr+s=P+Q. (3.15) 


To get a surviving coefficient F,,,, or G,,,, for a term of degree D, we must satisfy 
the first of the following equations and one of the second: 


P+Q=D, P—Q=mn+1 (m = 0, +1, +2, ---). (3.16) 


Equivalently, 
P= 5 (D + 1+ mn), O= 4 (D — 1 = mn). (3.17) 


rt =. ~ * (3 ¥(3) . . ‘2 
lo get a surviving coefficient F,,,, or G,,,. for a term of degree D, we must satisfy 
the first of the following equations and one of the second: 


P+Q=D, P—Q=mn (m = 0, +1, +2, ---). (3.18) 


Equivalently, 


P= 5 (D + mn), Q = 5 (D — mn). (3.19) 


The equations (3.16) and (3.18) form a sieve by means of which we reject certain 
terms from the series (3.3), the existence of the rejected terms being in fact inconsistent 
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with the assumed symmetry. It is difficult to describe in words the detailed consequences 
of the operation of this sieve. Some statements will be made below, but the most rapid 
way to see the results is to consult the diagrams (Figs. 2a-n). These show by heavy 
dots the surviving terms for symmetries ranging from digonal to octagonal. 

To interpret these diagrams, let us take Fig. 2e as an example. This refers to the 
cross force F and cross torque G in the case of tetragonal symmetry. P and Q are plotted 
as rectangular Cartesian coordinates; since neither P nor Q can be negative, only the 
first quadrant is shown. The two families of lines are those given by (3.16) for D = 
0, 1, 2, --- and m = 0, +1, +2, --- . The dots are those intersections of these lines 
which occur at integer or zero values of P and Q. Lines, such as P + Q = 2, which give 
no such intersections, are not shown. 

We may analyze the diagram by proceeding through increasing values of the degree 
D. There is no term of zero degree, nor indeed of any even degree; this is consistent 
with the general Theorem I, given later. 

For terms of the first degree, we have D = 1, and an intersection occurs at P = | 


and Q = 0. Thus, by (3.11), we have surviving terms with 


p=l, 7] 0, y = @, g = @, 


and 


Consequently there are linear terms in the expansions of F and G of the forms 


F coo ts » Ge) + Pool tls , Gs)e; 
(3.20) 
G jo00(Usg , @3)U + Gooro(Us , &3)H, 
respectively, and no other linear terms. 
For terms of the third degree, we have D = 3, and intersections occur at 
(P = 2,Q = 1) and (P = 0, Q = 3). Thus we have surviving terms for the following 


values of p, 9, 7, 8: 
2 | 0 0 
P = 2, QO= ] 
2 0 0 if 2 


0 Oo 2 1 J 


0 38 0 o|] 
0 2 O hOB 

-P=0 Q=3 
»o 10 2f{? , d 


0 0 0 3 


Thus the following terms of the third degree occur in the expansion (3.3) of F, with 


similar terms in the expansion of G: 
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F, (Us ; ws; uu + Fyi10(Us , Ws )Ullw + Foi 20(Us ’ ws ier 


+ F500: (Ug ’ ws)U@ + Pyoii (Ug » Ws Ua + Poo (us , @3)0 @ 
(3.21) 


2 is F os00(Us ’ w;)u + Foo (Us ’ w;)u @ + Foio2(us ’ w3)itar 


, 3 
+ Foo03(Ug » W3)@ - 


There are no other surviving terms of the third degree in these expansions. 

In the case of the diagrams for the axial components F; , G; , the lines shown have 
the equations (3.18). The reader should have no difficulty in interpreting the other 
details of the diagrams. Let us now make some general statements with regard to the 
consequences of n-gonal rotational symmetry. 

First, let us consider consequences of evenness or oddness in the order n of the 
n-gonal rotational symmetry. The following results follow easily from (3.17) and (3.19): 

[. In the case of n-gonal rotational symmetry of EVEN order (digonal, tetragonal, 
hexagonal, etc.), the series (3.3) for F and G contain only terms of ODD degree, and the 
series (3.3) for F, and G, contain only terms of EVEN degree. 

Il. In the case of n-gonal rotational symmetry of ODD degree (trigonal, pentagonal, 
etc.), the lowest EVEN degree occurring in the series (3.3) for F and G is (n — 1), and 
the lowest ODD degree occurring in the series (3.3) for F; and G, ts n. 

Secondly, let us consider the survival of terms of low degrees (D = 0, 1, 2). 

If D 0, then P = Q = O. In the case of F and G, we have to satisfy the second 


of (3.16) in the form 0 = mn + 1. If n # 1, this cannot be satisfied by any m in the 
range 0, +1, +2,--- .Ifn = 1, it is satisfied by m = —1. Thus we have the result: 

III. The absolute terms in the series (3.3) for F and G survive only in the degenerate 
case n |, which corresponds to no real symmetry at all. 


In the case of F; and G, we have to satisfy the second of (3.18) in the form 0 = mn. 
This can always be done by taking m = 0. Hence we have the result: 
[V. The series (3.3) for F; and G; always have absolute terms. 


Now let us consider terms of the first degree, D = 1. Then either 
P = 1, Q = 0, (3.22) 
or 
P = 0, Q=1. (3.23) 


In the case of F and G, we have to satisfy the second of (3.16). Corresponding to (3.22) 
and (3.23) respectively, this gives 
mn = 0, (3.24) 
and 
mn = —2. (3.25) 
Equation (3.24) is satisfied for any n by taking m = 0. Equation (3.25) can be satisfied 
only if n = 1 or n = 2, in which cases it has solutions m = —2,m = —1, respectively. 
Noting that (3.22) gives 
p=1, q=0 r = 0, s = 0, 
p = 0, q = 0, r=1, s = 0, 
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and (3.23) gives 
p= Q, g= l, rp = @, 3 = @, 


p= 0, q= 0, r = QO, s =], 


we may state the following result: 
V. For any n, the series (3.3) for F and G contain linear terms of the form 
F’ ,o00(tes ; W3)U + Poa Uz , 3) W, 
(3.26) 
G 000! U3 ] W3)U = Goo10( U3 ? W3)W. 


If n > 2, these are the only terms of the first degree; if n = 1 orn = 2, there are additional 


line ar te rms of the form 


P nanttie , W3)U + Foo01 (Us ’ W3)@, 
(3.27) 


Go100(Us ; w3)u + Gooo1 (Us » G3 )@, 


respectively. 
To discuss terms of the first degree in the series for F; and G; , we go back to (3.22) 
and (3.23), and combine them with the second of (3.18). We get respectively, 


mn = l, (3.28) 
or 
mn = —1. (3.29) 
These can be satisfied if and only if nm = 1. Thus we have the result: 
VI. In the series (3.3) for F, and G, , linear terms occur only in the degenerate case 
ft = | 
Finally let us consider terms of the second degree, D = 2. Now we have three alter- 
natives 
P=2 Q=0, (3.30) 
P= 1 Q = (3.31 
on 
P=0, Q=2. (3.32) 


In the case of F and G we have correspondingly from the second of (3.16): 


mn = lI, (3.33) 
Ol 
mn —1], (3.34) 
or 
mn = —3. (3.35) 
The first two can be satisfied if and only if n = 1. Equation (3.35) can be satisfied if 
and only if n = 1 orn = 3. Thus we have the result: 
VII. The series (3.3) for F and G contain terms of the second degree if and only if 
n l orn = 3 (no symmetry or trigonal symmetry). If n ¥ 1, all terms of the second degree 


involve UW, Ue, @ only. 
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To discuss terms of the second degree in the series for F; and G; , we use (3.30), 
3.31), (3.32) with the second of (3.18). We get the following corresponding equations: 


mn = 2. (3.36) 
or 

mn = 0, (3.37) 
or e 

mn = —2. (3.38) 


Equation (3.37) can be satisfied for any n; the other two equations can be satisfied if 
and only if nm = 1 orn = 2. Hence we have this result: 
VIII. If n > 2, the series (3.3) for Fs and G, possess terms of the second degree of the 


following forms only: 


F root , @, UU oo Fionn (Us 9 w3)Ua@ + Foir0(Us ; w; uw + Foi (us ’ Ws ww, 
(3.39) 


Giioo(Us , ws)uu + G3}, (Us , ws)U@ + Gorio(us , ws) + Goorr(Us , ws)wo. 
Ifn = l orn = 2 (no symmetry or digonal), all terms of the second degree occur. 
The following result may be added as an immediate deduction from the second of 
3.16): 
[X. Except in the degenerate case (n = 1), there occur in the series (3.3) for F and G 
no terms making 


ptr=qtes. 


This means that, in the diagrams of Fig. 2 which refer to F and G, we find no dots on 
the main diagonal, P = Q. 

4. Consequences of reflectional symmetry. Let us now suppose that the projectile 
has reflectional symmetry in a plane P. If there is also n-gonal rotational symmetry 
about an axis A, we shall suppose that A is contained in P. 

Let us choose the origin O in P (and on A if there is n-gonal rotational symmetry). 
Let us choose the axis Ox; in P (and along A if there is n-gonal rotational symmetry). 
Finally, let us choose the axis Ox, in P, so that Ox, is perpendicular to P. 

[he covering operation is a reflection in P. Applied to the axes, this operation 
changes Ox,2.2, into Oxjxixri , where Ox! , Ox} coincide with Ox, , Ox; respectively, and 
Ox! is Ox. reversed. 

Consider any state of motion of the projectile. This may be described by the com- 
ponents (uw; , Us , Us , @; , 2 , #3) along Ox,2.x2, or by the components (uj , uz, Us , 
wi , wh , wf) along Oxixixi . It is understood that a component of angular velocity in 
a given direction is positive if it corresponds to a positive rotation about that direction. 
Although we may be accustomed to defining a positive rotation as a right-handed 
rotation, that definition must not be used here, because we are considering two sets 
of axes with different orientations. We accept the definition that a rotation is positive 
if it corresponds to a cyclical rotation of the axes. Thus ws is positive if it tends to turn 
Ox, into Ox, , and w$ is positive if it tends to turn Ox{ into Ox; . It is then easy to see 
that when we reflect the axes in accordance with the covering operation described 
above, and describe one single state of motion by components along the two sets of 
axes, the two sets of components are related by the following formulae of transformation: 
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“= Uw, i = —U, Uu;= Us; 
(4.1) 


.— Pe 
so “Ms 2 = We; Wz = —@3. 


To this state of motion there corresponds an aerodynamic force system, which may be 
described by the components (F, , F2 , F; , G; , Gz , Gs) along Ox,x.2, or by the com- 
ponents (F{ , Fi, F3 , Gi, G2, G3) along Ox{aiar4 . As in the case of angular velocities, a 
component of torque is positive if it tends to produce a cyclical rotation of the axes. 
Thus G; is positive if it tends to turn Oz, into Ox, , and G3 is positive if it tends to turn 
Ox; into Ox; . We have then the following formulae of transformation analogous to (4.1): 


R= Fh, Kh=-F, Fh= Pr, 
Gi = -G,, mZ2= G., Gs = —G;. 


The formulae (4.1) and (4.2) are most easily remembered by the following rule: 
When the axes are reflected, polar vectors (velocity and force) undergo ordinary re- 
flection, but axial vectors (angular velocity and torque) undergo ordinary reflection 
followed by a reversal of sense. 

In the complex notation of (2.2), the transformations (4.1) and (4.2) read 
“ue = 4, un = Us, , w’ = —o, 03; = —w;3, 

: (4.3) 
F’ = F, FP; = F;, G’ = —G, G; = —G;. 


The above formulae refer to the effects of reflection of axes; the idea that the pro- 
jectile is symmetric is not involved. This idea we now introduce, and substitute from 
(4.3) in (2.4). We see that reflectional symmetry in the plane Ox,x; implies the following 


identities: 


fu, u, us , —®, —w, —ws) = = f(u, U, Us , w, @, ws), 
gU, U, Us , —®, —w, —ws3) = —G(U, U, Us , w, @, ws), 

. <. ~ 7 i (4.4) 
f3(U, U, Us , —@, —w, —w3) = fa(u, U, Uz , w, @, ws), 
g3(U, u. Us 9 —O, —@W, —ws) = —93(u, u, U3 7 W, a, Ws). 


Here a bar over a functional symbol implies that the sign of 7 is changed throughout, 
i.e. both in the arguments and in any complex coefficients which occur in the function. 
Since f,; and g; are real, it is of course a matter of indifference whether we use bars for 
them or not. 

We now assume that the functions f, g, fs , g; admit series expansions as in (3.3). 
Substitution in (4.4) gives 


Zz; _- 1)" F yers(Us ’ — ws) Uw" = p F yqra(Us ’ w3)U Uo", 

pers pare 

} ® (- Ly: . "Goer (Us ; —w;)0u'o'w* =_— = > Gyore(Us , 3) U'o'w", 

pars pare (4.5) 
= (4.5 

DY (-1)*F2..us , —ws)Pua'w’ = >> FS (us , ws)u'a'w’, 

pare pers 

> (-1)'"'42,(us, , —ws)u'o’w’ = — D> G2. (us , ws)u'a'a’. 


pare pars 
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Since these may be treated as identities in u, U, us , w, @, w; , we deduce 


Feva(tts » —ts) = (—1)"**F (tts , Ws), 
Grars(Uzs » —W3) = —(—1)"*"Gyere(Us » a), 
FS (us, —ws) = (—1)"**F2),(us , ws), (4.6) 
7?) (us, —ee) = —(—1)'*G ats , a), 


(p, g, r, 8 = 0, 1, 2, ---). 


We note that the exponent (r + s) is the degree in cross spin of a term in the ex- 
pansions (3.3). 

To show the implications of (4.6), we introduce the symbols FE and O to represent 
real functions of us; and w; which are respectively even amd odd in w; . We may then 


show the results as follows: 


(r + s) even (r + s) odd 
| and . E+ i0 O+ iE (4.7) 


a and G, 


nell O + ik E + i0 

In making the above deductions from the assumption of reflectional symmetry, we 
have made no use of the existence or non-existence of n-gonal rotational symmetry. 
In fact, the consequences of the two types of symmetry are essentially different in 
nature. As we saw in Section 3, n-gonal rotational symmetry eliminates terms in the series 
(3.3). Reflectional symmetry does not eliminate any terms; it merely gives us informa- 
tion with regard to odd or even functional dependence on axial spin w; in the coefficients 
in (3.3). But in the special case discussed below in Sees. 6 and 7, reflectional symmetry 
does enable us to eliminate further terms from the expansion G; . 

5. Case of vanishing cross spin (w = 0). Let us now consider the case where the 
projectile has no cross spin. This may be either an instantaneous condition in a general 
motion, or it may be a permanent condition enforced by constraints on the projectile. 
Thus, if the projectile is mounted in a wind tunnel in such a way that its axis is fixed 
in direction, then the cross spin vanishes. The projectile may be completely fixed, or 
it may be free to rotate about its axis. We note that, for wind tunnel discussions, the 
free air stream is to be used as frame of reference. 

For vanishing cross spin, we have w, = w, = 0, and sow = w = 0. Thus all terms 
disappear from the series (3.3) except those for which r = s = 0. In fact, the series read 


flu, U, Us ’ 0, 0, Ws3) = } F ,q00(Us ’ ww iu’, 
pe 


glu, u, Us ; 0, 0, ws) = >» G,.00(Us ; w;,)wu', 
pa 


(5.1) 
f(u, U, U3 ; 0, 0, @s) = ) 2 Feeho(Us , w3 uu", 
Pa 


D> G2ro(us , ws)Uu. 


p@ 


g3(u, U, us , 0, O, ws) 
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We shall suppose that the projectile has reflectional symmetry in the plane Oz,7; 
ang also n-gonal rotational symmetry about the axis Oz; . 
Since r + s is even, it follows from (4.7) as a deduction from reflectional symmetry 


Fo =E+i0, F®,=E+ 0, 
(5.2) 


ome 3) . 
G200 = O + iE, F5q00 = O + iE. 


Turning to (3.11), we note that P = p and Q = q. Figure 2 shows the coefficients which 
survive in the case of n-gonal rotational symmetry, but now we may read p instead 
of P and q instead of Q on the axes of the diagrams. As already noted in Theorem IX, 
on the main diagonal (p = q) there are no surviving coefficients in the series for F and 
G (except in the degenerate case of no symmetry (n = 1), which is not shown in the 
diagrams). But in the series for Ff; and G, there are terms on the main diagonal in all 
cases. The coefficients are of the forms 


(3 (3) ie 9) 
Fyp00(Us ’ ws), Gs po0(Us » Ws). (5.3) 


From the reality of F; and G; , it follows (as seen after (3.5)) that the above coefficients 
are real. When this information is compared with (5.2), the truth of the following state- 
ment is obvious: 

X. For a projectile possessing reflectional symmetry and n-gonal rotational symmetry, 
the coefficients (5.3) are respectively real EVEN and real ODD functions of ws . 

6. Case of vanishing axial spin (w; = 0). We next consider the case where the pro- 
jectile has no axial spin (w; = 0), but may have cross spin. This may be an instantaneous 
state of motion for a fin-stabilized projectile in free flight, or a permanent state enforced 
by constraints to prevent axial spin for a projectile mounted in a wind tunnel. 

Since an odd function of w; must vanish with w, , it follows from (4.7) that, when 
w; = 0, the coefficients take the forms: 


(r + s) even (r + s) odd 


F and ) poll k(us) ik(us) (6.1) 


pars 


G and «el ik(us) k(us) 


pars 


where k is used to denote real functions of u; ; k does not, of course, represent a single 
function. These results are consequences of reflectional symmetry, which assigns odd 
or even dependence on w; , and of the condition of vanishing axial spin, which eliminates 
either the real or the imaginary part for any given coefficient. 

Since the coefficients in the expansions of F and G are initially to be regarded as 
unrestricted complex numbers, the terms which survive n-gonal rotational symmetry 
also survive after application of the above two assumptions. The character of the re- 
striction placed on the coefficients in the series for F and G due to vanishing axial spin 
may be summed up in the following statement: 

XI. For a projectile possessing reflectional symmetry and n-gonal rotational symmetry, 
and having no axial spin, the coefficients Pyar. , Gpare are respectively pure real and pure 
imaginary functions of us; for (r + s) an even integer. If (r + s) is odd, the situation is 


reversed. 











1949] AERODYNAMIC SYMMETRY OF PROJECTILES 357 


Further, it is easily seen that the assumption of vanishing axial spin does not eliminate 
any terms in the expansion of F, . From (6.1) it is seen that these coefficients F;3), are 
pure real or pure imaginary functions of u; depending on whether (r + s) is an even 
or odd integer. 

The fact that a projectile possesses reflectional symmetry and no axial spin, when 
combined with the result (3.5), does eliminate additional terms in the expansion of the 
axial torque G; for the projectile. By (3.5), G53?,(us , #3) are real, and so, in particular, 


G)p--(Us , 0) are real. But, by (6.1), GSS?,(ug , 0) are pure imaginary. Hence 
Joe r(Us , 0) = O. (6.2) 


Thus we have this result: 
XII. For a projectile possessing reflectional symmetry and n-gonal rotational symmetry, 
and having no axial spin, the series (3.3) for G, contains no terms for which p = q and 


Let us examine this result, proceeding by increasing degree D. 
If D 0, we have the coefficient Gooo0(usz , 0), and it vanishes by (6.2). 
If D 1, we have the coefficients 


Gyooo(us , 0), Goreo(ts , 0), Goor0(Us , 0), Gi . 


Obviously, (6.2) tells us nothing about these coefficients. However, as we saw in Theorem 
VI, these coefficients occur only in the degenerate case n = 1. 


(3) 


If D = 2, the subscripts of G7,.(u; , 0) may be written as follows: 


2000 
0020 
*1100 


1001 
(6.3) 
0110 


*O011 
0200 


0002. 


Of these, those marked with a star vanish by (6.2). We are left with 6 terms. Combining 
this with Theorem VIII and (6.1), we have the following result: 

XIII. For a projectile possessing reflectional symmetry and n-gonal rotational sym- 
metry, and having no axial spin, the terms in the series (3.3) for the axial torque G; may be 
described as follows: 

a. There is no absolute term. 

b. There are linear terms only if n = 1 (degenerate case). 

c. Ifn = 1 orn = 2, there are six terms of the second degree, with subscripts 
as shown by the unstarred entries in (6.3); if n > 2, there are only two terms 
of the second degree, which may be written 
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k(us)(ua + iw) (6.4) 


where k is a real function. 
The last deduction follows from (3.37), since, by (3.5), 


Gi3,(us , 0) = Gero(us , 0), (6.5) 


and, by (6.1), Gido(us , 0) is real. 
We note, that for n >, 2, the expansion of G; reads 
G; = k(us)(uw + Uw) + (terms of higher degree). (6.6) 
7. Case whére cross spin and axial spin both vanish (w = 0, w, = 0). We shall 
now bring together the situations discussed in Secs. 5 and 6 by supposing that the 
axial spin w, and the cross spin w vanish simultaneously. We discuss here only the axial 
torque G; . 
As a particular case of (6.2), we have 
Grrio(us , 0) = O (p = 0, 1, ---). (7.1) 
This means that when w = w; = 0, for a projectile possessing reflectional and n-gonal 
rotational symmetry, we are to delete the main diagonals in the charts of Fig. 2 for G; ; 
these coefficients no longer survive. Under these circumstances, let us seek the surviving 
coefficients of lowest degree in the series (3.3) for G; . It is then a question of finding 
the smallest non-negative integer D to satisfy (3.18) with the additional condition 
P # Q, since the main diagonal has been eliminated. 
The equations to be solved read 


P+Q = D, P—Q= mn (m = 0, +1, s+), 


P # Q. 
These imply 
D? = (P + Q)’ = (P — Q)? + 4PQ = m’n’ + 4PQ. (7.3) 


The possibility m = 0 is excluded, since by (7.2) it would give the forbidden relation 


P = Q. Hence the smallest D is given by taking m* = 1 and either P = 0 or Q = 0. 
This smallest value of D is n, and the solutions giving the smallest D are 


cee F Q = 0, m = 1; P = Q, Q =n, m= -l. (7.4) 


Let us summarize as follows: 

XIV. A projectile possesses reflectional symmetry and n-gonal rotational symmetry. 
The cross spin and the axial spin both vanish. Under these conditions the series (5.1) for 
the axial torque G; starts with terms of degree n, and the series is of the form 


Gz = —tk(us)(u” — u") + terms of degree n + 2. (7.5) 
To complete the proof of this statement, it follows from (7.4) that 
G; = Girolus , Ou" + Gorro(us , Ou” + terms of degree n + 2. (7.6) 


The coefficients are pure imaginaries by (4.7) and complex conjugates by (3.5), since 
G; is real. Thus, in terms of a real function k(u;), we may write (7.6) in the form (7.5). 
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The angle of yaw (6) is the angle between the velocity of the point O and the axis 

of the projectile. Thus 
tan 6 = |ul/us. (7.7) 
Let @ denote the angle between the cross velocity vector u and the axis Ox, , which it 
will be remembered lies in a plane of reflectional symmetry of the projectile. We shall 
count @ positive when the direction of u is generated from Ox, by a positive rotation 


(Fig. 3). Then 


u=|ule™. (7.8) 
Combining (7.7) and (7.8), we have 
u = Us; tan 6-e"”. (7.9) 
Substitution in (7.6) gives 
G; = 2uzk(us) tan" 6 sin n@ + terms of order tan”*? 6. (7.10) 


To summarize: 

XV. The formula (7.10) shows explicitly the dependence of the principal part of the 
axial torque G; on the yaw 6 and the angle 6 for a projectile possessing reflectional symmetry 
and n-gonal rotational symmetry, in the case where the projectile has no cross spin or axial 
spin. ‘ 

As was to be expected, G, is periodic in @ with period a = 22/n, the angle of the 
covering operation. 

The formula (7.10) holds even in the degenerate case n = 1, when we have no rota- 
tional symmetry, only reflectional symmetry. 

8. Stability. An interesting question concerning the generation of axial spin arises 
in respect of (7.10). Suppose that a projectile possessing reflectional symmetry and n- 
gonal rotational symmetry has no cross spin or axial spin (w = w; = 0). Then G; is of 
the form (7.10), from which we see that G; = 0 if either (a) the cross velocity lies in 
the plane of reflectional symmetry Oz,2,(@ = 0), or if (b) there in no yaw (6 = 0). 

Let us consider the case in which the projectile is yawing (6 # 0) so that the vanishing 
or non-vanishing of the axial torque depends only on @. We shall then have G; = 0 
for 6 = mx/n, where m = 0, +1, +2, --- . For these values of @, the cross velocity 
vector lies in a plane of symmetry. If the projectile has n-gonal rotational symmetry 
and one plane of reflectional symmetry (containing the axis of rotational symmetry), 
then it necessarily has n planes of reflectional symmetry. The angles between these 
planes are 7/n. 

We must observe that, for a given projectile, the axes Ox,x.x3; are not uniquely 
determined. Choosing Oz; along the axis of the projectile in the general sense of flight, 
and deciding to use right-handed axes, we have still a choice of 2n consecutive directions 
differing by an angle +/n for the Oz, axis. For a given motion, a change from one such 
choice of Ox, to the next alters 6 by /n, and changes the sign of sin né@ in (7.10). Since 
G; is not altered, such a change reverses the sign of the function k(us). We must there- 
fore be careful not to attach importance to the sign of k(u;) until we are sure what 
axes we are using. 

Suppose then that we have made some definite choice of Ox, , as in Fig. 3. Consider 
various directions of the cross velocity u. As we rotate u, we get successive positions of 
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vanishing axial torque, namely, whenever u crosses a plane of symmetry. These will be 
alternately stable and unstable, in the sense that if the vector u is displaced slightly 
from one of these positions, the axial torque consequent on this displacement will tend 
to decrease or to increase the angular displacement between u and the plane of sym- 
metry originally coincident with u. The criterion for stability is that the generated 
axial torque should rotate the projectile so as to bring the plane of symmetry from 
which uw was displaced back into coincidence with u. This requires a positive G, for a 
positive increment in 6; in other words the criterion for stability is 


dG,;/dé > 0. (8.1) 


Referring to (7.10), we see that we may make the following statement: 

XVI. If k(us) > 0, the projectile is stable with respect to axial spin when yawed in 
the planes of symmetry @ = 2mx/n (m = 0, 1, 2, ---), and unstable when yawed in the 
planes of symmetry 0 = (2m + 1)x/n (m = 0, 1, 2, ---). If k(us) < 0, the situation is 
reversed. 

We have no a priori knowledge about the sign of k(u;). Suppose the projectile is to 
consist of a body of revolution and fins as shown in Fig. 1, the broken lines being included 
so that the projectile has reflectional symmetry. Let us for the sake of definiteness take 
Oz, in a fin. 

In the degenerate case (n = 1), where there is no symmetry except reflectional 
symmetry (Fig. 4a), we expect stability when the fin is downstream. We note that u 
represents the cross velocity of the projectile relative to the air, so that —wu represents 
the cross wind relative to the projectile; we may expect that stability occurs when 
@ = z, as shown by the vector marked S in Fig. 4a. Thus, according to the above 
criterion for stability, in order that (8.1) holds for 6 = w we must have k(u;) < 0 when 
n = 1. Figure 4a also shows the relation between the angle @ and stability. The shaded 
portion of the plane is a region of stability in the sense that if the cross velocity vector 
u once enters the shaded half-plane it can never leave and tends to settle on the position 
S. If w is initially in the unstable region (right half-plane) the projectile will rotate so 
that u leaves the unstable region by the shorter route. These statements are not mathe- 
matical deductions; they represent what, from general experience, we would expect to 
happen in the case of a projectile with a single fin. 

In the case n = 2 (Fig. 4b), it would generally be agreed that the projectile would 
turn so as to set the face of the fins at right angles to the cross wind. That is, taking 
Ox, in one of the fins, the two positions for u given by @ = 2/2, 32/2 are positions of 
stability and the positions 6 = 0, z are positions of instability. In order that this situation 
may result from the above criterion of stability, it is necessary that k(u;) < 0 when 
n = 2. 

In both the above cases, it is interesting to note that the positions of stability are 
those for which the cross velocity u is coincident with the bisector of the angle between 
adjacent fins. The question then arises as to whether this holds true for projectiles 
having more than two fins. This of course could be answered if there were an a priori 
method for the determination of the algebraic sign of k(u;). In the absence of such a 
method, experiment alone can decide. To explain the point at issue, let us discuss the 
ease of a bomb with three fins (n = 3), as shown in Figs. 4c, d. 

Let us choose the axis Oz, in one of the fins. Two cases can arise, k(u;) > 0 and 
k(us) < 0, if we omit from consideration the particular case k(u;) = 0. 
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Let us consider first the case k(u,;) > 0. From Theorem XVI we infer that the pro- 
jectile is in a position of stable equilibrium when 6 = 0, 27/3, 41/3, i.e. when the cross 
velocity vector u lies in one of the fins (Fig. 4c). Also we see that there is unstable 
equilibrium when 6 = 2/3, x, 52/3, i.e. when the cross velocity vector u bisects the 
angle between two fins. The shaded regions in Fig. 4c indicate regions of stability, in 


ae a cae 
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Fig. 1 
Diagram showing cross-section with n-gonal rotational symmetry, with and without reflectional sym- 


metry. (If dotted portions are included, we have reflectional symmetry; if they are omitted, we have no 
reflectional symmetry). 
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the sense that, if wu lies in one of them, the projectile will rotate so as to bring u into 
coincidence with one of the fins, and there will be an oscillation of u (in the shaded 
region) about that fin. 

On the other hand, if k(us) < 0, the situation is reversed as shown in Fig. 4d. Now 
the stable positions of equilibrium are those for which u bisects an angle between fins. 

We have then the following interesting question, capable of answer by wind tunnel 
experiments: If a three finned projectile is mounted in a wind tunnel so that it is free 
to turn around its axis, are the stable positions of equilibrium those in which a fin 
lies in the plane of yaw (Fig. 4c), or are they those in which the plane of yaw bisects 
the angle between two fins (Fig. 4d)? The same question may of course be asked for the 
case of four fins, or indeed for the case of any number of fins. It would indeed be in- 
teresting to know whether the situation discussed for n = 1 and n = 2 continues for 
higher values of n. It may be that for higher values of n the screening effect of the body 
of the projectile becomes of significance, and no general statements can be made covering 
all projectiles which have the same value of n. 
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‘The heavy dots indicate surviving terms in the expansions (3.3) for various types of rotational symmetry. 
(F, G) are cross force and cross torque respectively. (Ff; , Gs) are axial force and axial torque respectively. 
For P and Q, see (3.11). D is the degree of the term. 
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Fig. 4c. n = 3, three fins, k(us) > 0. 








Fig. 4d. n = 3, three fins, k(us) < 0. 


Fia. 4 


Stability of a projectile with respect to axial spin. S indicates equilibrium directions of u (cross velocity 
of projectile) which probably correspond to stability. J indicates equilibrium directions of u which 
probably correspond to instability. Shaded regions are regions of stability. 
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ON THE STABILITY OF THE SUPERSONIC FLOWS PAST A WEDGE* 


BY 
G. F. CARRIER 


Brown University 


1. Introduction. When an ideal’ compressible gas approaches a symmetrically oriented 
infinite wedge with uniform supersonic speed, there exist for many wedge angles and 
upstream velocity magnitudes two alternative steady flow configurations. Associated 
with one of the flow fields is a plane shock which is attached to the wedge vertex and 
downstream of which is a uniform subsonic flow. The alternative solution includes a 
“‘weaker’’ attached shock which is more sharply inclined to the wedge surface. Down- 
stream of this shock there also exists a uniform flow field with a larger velocity than 
that of the strong shock solution. Usually, in this latter case, the downstream flow is 
supersonic but there is comparatively small range of “‘inlet’’ conditions for which it is 
subsonic. The weak shock is usually observed in experimental work although with 
sufficient care one can observe the strong shock configuration.” 

It seems of interest, then, to see what a perfect fluid theory will predict in regard 
to the stability of these two flows when they are subjected to small time-dependent 
disturbances. 

In the opening sections of this paper we shall be concerned with the interaction of 
small plane disturbances and plane shock waves. The results of both this plane wave 
investigation and an independent alternative procedure will then be used to show that 
both configurations are stable. 

2. The wave equations. The fundamental equations which define the motion of an 
ideal compressible fluid are those implying the conservation of momentum, mass, and 
energy. A convenient form in which they may be presented is the following [in Eq. (3) 
we have already used the fact that p = p.RT): 


(V - grad) V + 8V/dt + po‘ grad p = 0, (1) 
V - grad (In po) + A(In p.)/dt + div V = 0, (2) 
(V - grad + 0/dt)(p/po) = 0. (3) 


Here V, p, and pp are respectively, the velocity, pressure, and density, at any point in 
the field of flow,® and y is the ratio of the specific heats. We shall be interested in the 
superposition of a small disturbance on a known steady flow. It is therefore convenient 


to write 
V = U, + a,(grad ¢ + curl E), (4) 
p = pl + 7q), (5) 
Po = px(1 + p). (6) 


*Received Sept. 6, 1947; revised March 1, 1948. 
1We imply vanishing viscosity and heat conductivity, and the validity of the ideal gas law, p = pRT. 


2A. Kantrowitz has successfully produced the strong shock solutions. 
*The subscript zero merely distinguishes the density from the variable without subscript which will 


soon arise. 
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Here, the subscript 2 characterized the state parameters associated with the steady 
flow downstream of the shock. The disturbance is characterized by the functions, ¢, 
q, p, and E. No justification is needed for the foregoing representation of the disturbance 
velocity except to note that any vector can be expressed in this manner. In fact, we may 
(and shall) still require without loss of generality that div E should vanish identically. 
Using now a perturbation process to linearize Eqs. (1), (2), (3), with regard to the 
disturbance parameters,’ we obtain the following linear differential equation 


Ay — d’y/dt = 0, (7) 
dE/dt = 0, (8) 
dS/dt = 0. (9) 


It is implied in the manipulations which lead to these equations that 
q = —dg/dt. (10) 


In Eqs. (7) to (10) (and henceforth) ¢ represents the time multiplied by a, ; M is the 
Mach number | U, |/a, ; A is the Laplace operator; d/dt = 0/dt + M,-grad; and S = 
a= 6. 

It is of interest to interpret these equations from a physical point of view. Equation 
(7) implies that the function g represents an acoustic wave which propagates according 
to a conventional wave equation. The vector E characterizes that part of the velocity 
vector which implies a non-vanishing vorticity and, according to Eq. (8), drifts with 
the stream. Finally, S, the linearized entropy, also drifts with the stream. The fact 
that these latter quantities propagate in this manner is, of course, well known [4]. 

3. The boundary conditions at the shock. In this section we shall be concerned 
with the formulation of the conditions which must apply at the downstream side of 
the shock when the steady flow field is disturbed in the immediate neighborhood of 
the shock. We shall be interested only in those situations where the upstream field 
remains undisturbed (i.e. either the disturbance originates in the downstream field or 
sufficient time has elapsed since its formation so that the upstream field is again free 
from any perturbations). Since, for a stationary shock, the Rankine-Hugoniot condi- 
tions’ uniquely determine the state variables on the downstream side in terms of the 
upstream state variables, no time dependent changes in the field can occur unless we 
admit oscillations of the shock locus which have a magnitude consistent with that of 
the disturbance functions. At any time ¢t, then, the shock is displaced from its steady 
state position by an amount® y = y(y, z, t). In order to determine the downstream 
state variables we must apply the Rankine-Hugoniot conditions to a shock moving 
with a velocity y, and oriented with regard to its equilibrium configuration at angles 
Wy, ,¥.. Here, we have implied the omission of higher order terms’ in y, , ¥; , --- . This 
linearization is essential if the procedure is to be consistent with that of Sec. 2. 


‘Details of such a process may be found in [4]. The most naive approach, however, gives the same 
result as a rigorous analysis. 

5See [1],- [2], [3]. 

*Here, we choose rectangular coordinate axes as follows: the axis z is normal to the plane of the shock, 
the plane-z = 0 is normal to both U; and U;. 

7The angles, of course, are actually arctan yy, , °°: 
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Before writing down these jump conditions, the following notation must be intro- 


duced.” 

U, = iu, + jv; U, = iu, + ju, ; 

u*, v¥, w¥, are the upstream and downstream components of velocity in- 
’ , . 


* y* 
Wi, 


/ 
/ 
/ 
| 
arctan v, 
vn. . aimee State Shock Position 
uy 
u ™~ 
‘ . hi / Steady State Velocities 


7___—_ Disturbed State Velocity 


/ La Steady State Position of Shock, x=AAy 


/ 
/ 
! 

x 


ANAS ANA ANA VAY es 








Fia. 1 


The coordinate system depicted here is applicable only to the work of Sec. 5. The earlier work 
utilizes z = 0 as the shock locus. In this diagram, of course, the z components of velocity are 
not visible. 


*i, j, K will denote the conventional Cartesian unit vectors. 
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stantaneously normal to and tangent to (in planes z = 0 and y = 0, respectively) the 
displaced shock locus; 5u2 , v2, dW: , 5p2 , Sp2 , are the changes in U2, °** , Po, associated 
with the shock displacement. From Fig. 1, it is readily seen that (after linearization 


with regard to y, , ¥.) 


us = u, — 1, ; ve = 9, + my, , (11) 
uF = U2 — rep, + be, vy = v2 + wep, + bre, (12) 
w* = u,¥, , wk = uv, + bu.. (13) 
The Rankine-Hugoniot conditions may now be written in the forms 
pi(ur - y:) = (p2 + 5p2)(ut — v.), (14) 
ie pilus — vi) = Po + (p2 + 5p2)(ut = vi); (15) 
yR(T, — T2)/(y — 1) + (uk — ¥.)°/2 
(16) 
— = es S / * 2/6 
= 4 (dp2/p2 — 5p2/p2) + (uF — wv.) /2 
(y — lps 
yf = oF ., wr = we. (17) 


If we now substitute Eqs. (11) to (13) into Eqs. (14) to (17), linearize with regard 
to all disturbance parameters (i.e. du. , 6p. , --: , ¥,), and utilize the identities’ p,u,; = 


Pola, *** , we obtain 


bu atus — (vy — lui ( Vo ) 
= = ee (1 — pi/pro)\ Ve —# (18) 
a2 G2 — U2 pi/Pa)\ We + > Vy , 
iV. 
“ = (1 — pi/p2)(ti/a2)¥, (19) 
6 Jo 5 
= = (1 — p;/po2)(u;/d2)¥. (20) 
dp Ur 2a + (y — 1)(us — WU) ; ( Vo ) 
SP — 2, Se gp ' = (21 
Ds Ge ou ( p:/P2)\ We + a Vy ) 
op op F U; — Us Ve 
——= -—(y — 1) — (1 - | —= ) (22 
Aaa (y —_" ( Pr/ Po\ We FO vy 22) 


The five quantities on the left sides of these equations must be associated with 
the values of grad ¢ + curl E, — dg/dt, and S, on the plane x = 0. In fact, if we now 
take y = 1.40, and define m = uz/az , n = v2/a, , and use the fact that’ u,u. = a®” — 
(y — 1)03/(y + 1), we can write the boundary conditions in the form 


°These identities, of course, are the steady state Rankine-Hugoniot Conditions. 


Here a* is the critical velocity. 
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[grad » + curl E],-, = (1 — p:/po)[iai(v, + nv,) + jay, + kKasy.], (23) 


[—de (!) = (1 - pi/ po)as(y, + ny,), (24) 
[S].-0 = (1 — pi/pr)as(yr + ny), (25) 
where 
a, = (2 + 4m’)/3(1 — m’), (26) 
a. = (5 + m’)/6m, (27) 
a; = —m(5 + m’)/3(1,— m’), (28) 
a, = —(1 — m’)/3m. (29) 


4. The interaction of a plane acoustic wave and a shock. Let us now consider the 
hypothetical situation where the semi-infinite space + > 0 is unobstructed so that a 
uniform field exists downstream of the shock located at x = 0. Furthermore let us 
consider solutions of Eqs. (7), (8), (9), subject to the boundary conditions (23), (24), 
25), which are of the form™ 


g = Af(t — & — wy — 2) (30) 
E = (iB + jC + kD)f(t — ox — my — £2) (31) 
S = Ff(t — ox — ny — <2) (32) 
v = Gf(t — ny — £2). (33) 


This form implies that the disturbance consists of a plane acoustic wave with wave 
fronts tx + ny + ¢z = const. and of similar vorticity and entropy waves. In order 
that these waves obey Eqs. (7), (8), (9), one sees by a simple substitution that 

2+ 7? +? — (1 — mt — nn)’ = 0 (34) 
and 

1 — mo — ny = 0. (35) 

Here, we have used the notation of the foregoing section in writing d/dt = 0/dt + 
m 0/dx + nd/dy. The boundary conditions at x = 0 are given by Eqs. (23). (24), (25), 
and in this case (supplemented by the condition div E = 0). They take the form 











[ =€ 0 f -n» -a(1—n) | 4 | 
—7 —¢ 0 o Qs B 
-—¢ ” —o 0 Q> Cc 1|=0 (36) 
0 a n t 0 D 
(l—m-—n 0 0 0 a;(1 — nn) IL G _] 





Jt will turn out that these solutions are certain eigenfunctions. We shall interpret them later. 
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and 
F = a,(1 — n7y)G. (37) 
It is readily shown that this implies B = 0, and that Eq. (36) can be reduced, with 
the aid of Eq. (34) and (35), to the form 


| 0 we 4 + 8m’ 
ie ‘i 5+ m 


| 0 ; ' = 0, (38) 
| | 0 —2m | 
where » = [1 — (yn? + ¢°)(1 — m’)/(1 — ny)’]’”. (39) 
The solution of Eq. (38) is 
w= —m + (1 — m’)/(5 + m’)"™, (40) 
and implies that 
£ = (1 — ny)[+ (5 + m’)'? — 2m/(1 — m’))]. (41) 


Before attempting to interpret this class of solutions, we must establish some supple- 


mentary facts. Using the relations” 


Uy, = a®*? — (y — I)e2/(y + 1), (42) 
(a,/a2)” = 1 + (y — 1)(m* + n’)/2, (43) 

anc 
Qai/(y +1) < ui top = Ww +03 < 2ay/(y — 1). (44) 


(the first inequality is implied by the second law of thermodynamics) we may conclude 
that the range of physically possible flows is given by 
0 < n* <a, (y — 1)/2y < m’* < 1. (45) 


This in turn implies that the bracket in Eq. (41) is always negative for such flows. 
Finally Eq. (40) implies that 1 — u° never vanishes or becomes singular; Eq. (39) then 
implies that (1 — nn)* never vanishes or becomes singular, and hence (1 — m7) is always 
positive. These conclusions require that — be always negative [see Eq. (41)] and we 
have established that no waves of the form given in Eqs. (30) to (33) can be downstream 
moving waves. That is to say, a negative — implies a wave moving in the negative x 
direction. 

The waves implied by the foregoing analysis can now be interpreted in the following 
manner. A plane acoustic wave with wave fronts & + ny + fz = ¢t impinges on the 
shock at x = 0 and reflects not as an acoustic wave but as an “entropy wave” S and 
a “vorticity wave” E. In general, of course, an acoustic wave moving into the shock 
would reflect as a combination of entropy, vorticity, and acoustic waves, and hence 
the foregoing solutions are those eigen-solutions for which the amplitude of the acoustic 
reflection vanishes. 


Here, ao is the stagnation acoustic velocity of the steady downstream field. 
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If one wished to consider the general reflection problem, the following procedure 
would be indicated. Eq. (30) would be replaced by 


aed A, fit — &2 — q — $2) + A,f(t — &r — wy — $2), 


where for given » and ¢ the values of &, and & are the two roots of Eq. (34). In this prob- 
lem the incident wave is represented by A,f(t — --- — ¢z) and thus A, and the function 
f are known. The boundary value problem of determining the reflection is treated as 
before, but now the equation replacing Eq. (36) is non-homogeneous and thus there 
are not eigenvalues to be found but rather one dbtains the coefficients A,/A, , C/A, , 


5. The wedge flow stability problem. Let us consider the foregoing field of flow 
where now, however, we insert a solid boundary along the plane my — nx = 0 (ie. 
along the direction of the vector U,). We are interested in the existence within the 
region above this plane (i.e. the wedge surface) of oscillations which occur without 
the aid of external excitation. The boundary conditions which now apply are those 
already formulated for the boundary x = 0, the condition that the velocity component 
normal to the wedge surface vanish at the wedge, and the requirement that the wave 
be outgoing. 

Let us now consider that an acoustic disturbance which vanishes outside of a region 
R (yo < y < y;) exists in the field. This disturbance may be considered as being com- 
posed of plane waves.'* Some of these waves are moving away from the origin. Within 
a finite time these waves leave the region R and neither they nor their reflections from 
shock and wedge can return. Any waves moving toward the origin reflect appropriately 
from wedge and shock and, within a finite time, reach the origin. Here the nature of 
the combined reflection is not evident. In fact, if the region R is to be disturbed for 
more than a finite time, the disturbance introduced at the origin must be self-sustaining 


once it is initiated. It is evident, then, that we must look for an eigenfunction of our 
problem in a form that defines the behavior near the origin. In the following, we shall 
not worry about convergence in the large. The solutions considered will be of a type 


such that the leading terms of the expansions determine the eigenvalue. This will turn 
out to be such a value that whether the function is defined over the complete region 
or not, the solutions desired do not exist. We shall first consider those flows for which 
the down-stream field is supersonic. 

It is convenient at this stage to use a coordinate system such that the x-axis is rotated 
about the z-axis into the upper wedge surface. Using this rotated coordinate system 


see Fig. 1), we may write solutions of Eas. (7), (8), (9), in the form 


“ f(x ,y) exp fic[t — Mx/(M* — 1)]} cos az, (46) 


_— A a ee er * oe = - 
E {iF (y) + jG(y)} sin az + kH(y) cos az} exp {it Mu var — ;|} (47) 


\ 


t — x/M — d\y/M(M* — 1)]} cos az. (48) 





S S(y) exp {te 


Here f(x, y) must satisfy 


f.. —f,,/(M’* —-1)+ kf =90 (49) 


8Pither a discrete set or a continuum of such waves may be needed here. 
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where k* = c’/(M® — 1)? + a’/(M* — 1), and E and S satisfy Eqs. (8), (9) independent 
of the dependence on y of F, G, H, S(y). We may also write 

v = g(y) exp {ic[t — Mry/(M? — 1)]} cos az. (50) 
Finally, defining 6? = M’ — 1, r? = 2’ — 6’y’, 6 = arc tanh 6y/z, solutions of Eq. (49) 
may be written in the form, 

f,(z, y) = J,(kr) cosh v6. (51) 

These functions are such that f, vanishes at y = O and is non-singular at the origin if 
the admissible values of v are governed by Re(v) > 0. This latter condition will be 
critical in determining the stability of the flow. Using Eqs. (46) to (48) and (50), the 
boundary conditions Eqs. (18) to (22) and the condition div E = 0 may now be expressed 
in the following form. 


icM 


fy, ¥) — x52 — 7 SOu wy) + Ay) — aGy) = Arg) + Big’(y), (52) 
fQy, y) + aF(y) + i“ H(y) = Ag(y) + Bg’(y), (53) 
—af(rky, y) — “ Gy) — F’(y) = Asg(y) + Bsg’(y), (54) 
] fan 7 fQy, y) — Mf.Qy, y) = Ayg(y) + Big’(y), (55) 

Sy) = Asg(y) + Bsg’(y), (56) 
= Fy) + Gy) — oHlly) = 0. (57) 


The number A, , B; , given here can readily be deduced from Eqs. (18) to (22) and are 
functions of two parameters” say m and n. 
Let us now consider the function 


=) 


f= daJ..,(, y) (58) 


B=0 
as defined in Eq. (51). This choice for f and the associated choices of F(y), --- S(y) 
are to satisfy Eqs. (38) to (43). We can anticipate that such functions will be appropriate 
only for certain eigenvalues v. Let us assume, for the moment that v + 0. Then we 
must write 


2 


F(y) - x bud, eu(kéy), 
Gy) — te Cyd y+,(kéy), 


where & = (1 + \” — M’)™”. 


“Actually, in Eqs. (18) to (22), they are given as functions of u1, U2, V2, a2 °** , but each of these 
can be written as u; = u;(m,n),°**:. 
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If we now substitute these functions into Eqs. (52) to (57) and segregate the equalities 
in the coefficients of J,,,(k — y) for each yu, we find that the equality in the coefficients 
of J,-,(kéy) constitutes a homogeneous condition for the eigenvalue. The remaining 
equalities constitute non-homogeneous systems from which the a; , b; , --- , forj > 0 
may be evaluated. 

If one carries out the foregoing process, the characteristic equation for v can be 


written (after algebraic reduction) 


8 sinh (vy — 1) A, B, 
= 0. (59) 

M cosh ( — 1) A, B, 
where 6 are tanh (8/d). It is also implied that the lowest order non-vanishing co- 
efficient b, corresponds to n | whereas the lowest order coefficient in H is d, . Thus 


the boundary condition at y = 0 can be satisfied only if H(O) = 0, i.e. only if Re(v) > 0. 
It can readily be shown that B, < 0, By, < 0, B/X < MB./BB, < 2/8, for all wedge 


flows.’ This leads to the fact that Re(v) < 0 for all such flows. It is therefore evident 
that the type of disturbance defined here does not exist. 
Before these considerations can be interpreted as implying the non-existence of 


eigenfunctions of the type sought (i.e sustained oscillations), one more detail must be 
supplied. It is quite possible that any such eigenfunction might have a branch point at 
the origin which was not of the type admitted by the Bessel function expansion.” For 
example, one might expect that, near the origin, the wave function behaves like 
r’(In r)h(@). That this is not the case can readily be shown. Without algebraic detail, 
in fact, the procedure is the following. We consider a solution of the differential equation 
7) in the form of an expansion whose leading term is r’(In r)h(@), and define appropriate 
functions E and S consistent with this expansion. Then, as with the Bessel function 
expansion, the equalities among the coefficients of the leading terms which are implied 
by the boundary conditions lead to a transcendental equation for the eigenvalue v. Again, 
as one can readily verify, Re(v) is negative’’ for physically possible wedge flows and 
thus this type of eigenfunction is not admissible. As a matter of fact, it is not even 
certain that the expansion considered here leads to a function which is defined throughout 
the field. If not, this consideration is superfluous but in any event an oscillation of this 
type cannot exist. 

At this point we must extrapolate these results and assume that no other type of 
branch, point will lead to admissible solutions. One can readily obtain the general result 
that for a branch point of type r"g(r), where g(r) is singular at the origin, Re(n) will 
be negative. Therefore, using this extrapolation (or, if the reader prefers, this conjecture), 
we conclude that the flow configurations associated with a supersonic downstream field 
are stable. This conclusion is, of course, based on the fact that no self-excited oscillations 
‘an exist which vanish identically near the origin (the plane wave results imply this) 
ind that no non-singular solutions obey the boundary conditions in the neighborhood 
of the origin. This, of course, is in agreement with experimental evidence. 


'sSee appendix for details. 
'6This possibility was first suggested to the author by 8. 8. Shu. 
\7Tn fact, the same »; are associated with this solution as were found for the preceding expansion. 
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The subsonic case is more easily disposed of. The definition 8B? = 1 — M’, r? = 
x’ + By’, 6 = arc tan By/2z, will lead to a set of differential equations and their solutions 
which are analogous to those found by the corresponding procedure in the supersonic 
vase. Again, a set of Bessel functions multiplied this time by trigonometric functions in 
6 will form a product series solution. However, in order that the waves be outgoing in 
the subsonic field, we must now use the Hankel functions H}°’. These are singular 
for all v, of course, and again this type of series is not admissible. Therefore, with an extra- 
polation of the same type we used before, we must conclude that the subsonic field is also 
stable. 

In view of the somewhat unusual nature of our results, it seems worth while to 
summarize briefly without recourse to the physical picture. In the supersonic (subsonic) 
case, once we have demanded solutions of Eqs. (7) in the form f(x, y)e'“', we have a 
hyperbolic (elliptic) second order partial differential equation with constant coefficients. 


The homogeneous boundary conditions require g, = 0 on @ = 0, Lig, ¢,) = O on 
6 = @ , an outgoing wave, and the overriding condition that ¢ be non-singular at the 


origin. Here, L represents a linear combination of the argument functions and their 
derivatives with regard to r, ¢, is the normal derivative of ¢. Actually, the first three 
conditions given here are sufficient to find a unique set of eigenfunctions in either case. 
The satisfaction of the “overriding’’ condition (or the lack thereof) then determines 
the admissibility of these eigenfunctions and in our case is not satisfied. Thus the tenta- 
tively chosen functions are not suitable and hence no sustained oscillations exist. 

6. Some general comments. It seems of interest now to speculate on the causes of 
the fact that the weak shock is usually observed in experimental work. Each of the 
flow configurations will be altered (as indicated in the fore-going) by any externally 
generated disturbances. In the supersonic case, such disturbances can get into the field 
only through the boundary layer and even then must propagate downstream within 
their Mach cone. On the other hand, disturbances originating at the trailing end of 
the finite wedge may go upstream in a subsonic stream and the forced oscillations of 
the strong shock may be sufficiently great to “‘push it over’ into the other flow. Essen- 
tially the same point of view may be stated in the form: A disturbance in the weak shock 
field has a very small probability of getting near the origin where the overall behavior 
is determined. In the strong shock field, large disturbances are continually fed towards 
the origin, and unless such disturbances can be minimized sufficiently, the transition 
to the weak shock flow will occur. 

7. A remark on a previous investigation. In a recent paper [5], J. I. 
treated this same stability problem also using the small time-dependent-perturbation 
theory. He, however, considers a sustained excitation proceeding from the upstream field 
toward the downstream field. Although it is true that random upstream disturbances 
might initiate downstream perturbations, the proper investigation must not admit 
sustained excitation but must consider only the eigenfunction perturbations. We must 
therefore consider the physical argument as presented in that paper as incorrect for 
the problems which are usually of physical interest. Furthermore, in the boundary 
conditions at the shock, the author utilizes the fact that the shock is wavy but neglects 
the fact that it is moving. That is, he omits the terms corresponding to y, in our Eqs. 
(18) to (22). Hence, the actual calculations cannot even be used to find the response 


kaka} 


of the downstream field to an upstream excitation. His conclusions, which imply in- 


Levinson has 
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stability of all strong shocks and many weak shocks, are consistently contradicted by 
experimental evidence.* 
APPENDIX 
The inequalities which imply Re(v) < 0 may be developed as follows 


5+ m* — (2+ 4m*)n’ 


be = Gy 30 — mM’ 
- 5 + m* 
B, = —mn 3(1 — m)’ 
A= n/m 


8 =n’ +m’ — 1 = M —-1. 
We shall use the inequalities 

1 > m* > 1/7, n> >1— m’. 
We have, then, 





R _ at 8m )n> — (5 + m’)(1 — _m’) (7m* — 1)(1 — m’) 0 

6.M(1 — m’*) 6M(1 — m’) / 
VB.B,) —28 = 8. (4 + 8m*)n® — (5 + m*)(1 — m*) — 26°m*(5 + m’) 
rs h Xd 23°(5 + m®)m’? 


6 (7m* ~ Iti = m’) 
a 26°n?(5 + m’) *o 


Hence MB,/XB, > B/X. 
Furthermore 


MB. hr (4 + 8m’*)n® — (5 + m*)(1 — m’*) — 2n*(5 + m’) 


By, BB. 2(5 al m yn? 


h —(6n* + m + 5)(1 —m Rw 


8° 2n°(5 + m’) 


and the inequality is completed. 
We may now define (—/B,) — sinh a and (—8B,) = cosh a, so that Eq. (59) 
reduces to 
sinh [(v — 1)@) + a] = 
or 


v0. = 0 — a. 


*Note added in proof: Certain discussions have occurred while this paper was in proof which indicate 
that the analysis of the foregoing problem is not complete. Although no eigenfunctions exist which are 
exponential in time, we must still investigate how an arbitrary initial disturbance grows with time. This 
problem has been solved for certain special cases and a continued treatment of this stability question 


will appear as Part II of this paper. 
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However, we have shown that tanh a > tanh 6 , so for tanh a < 1, a > 6 and v is 
negative. For tanh a = Q > 1, we have a = 1/2 In (1 + Q)/(1 — Q) and Re(a) = 1/2 
In (1 + Q)/(Q — 1) = arc tanh 1/Q. 


Since we have shown Q™* > 8/) we again have Re(v) < 0. 
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HOMOGENEOUS HARMONIC FUNCTIONS* 
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1. Introduction. We consider, in the following, solutions of the Laplace equation 
G:z + Oy + Oe: = 9 (1.1) 
which are homogeneous of degree n in the coordinates, that is, such that 
g(Cz, Cy, Cz) = C"¢(a, y, 2) (1.2) 


for any constant C. The case n = 0 is of special interest and is considered at great 
length. As will be explained in the following paper, with proper change of variables, this 
case finds an important application in the supersonic conical flows. 

The relation (1.2) prescribes the variation of ¢g along lines through the origin; this 
relation may therefore be used to reduce the number of independent variables in (1.1), 
converting it to a differential equation in two variables. Transforming (1.1) to spherical 
coordinates, we shall use (1.2) first to eliminate 


R = (2? + y” + 2°). (1.3) 
thus converting (1.1) to a proper differential equation over the unit sphere R = 1; and 
then to eliminate z, converting (1.1) into a differential equation in the plane z = 1. 

For n 0 it is shown that by introducing isothermal coordinates along the unit 


sphere we may transform (1.1) into a two-dimensional Laplace equation. In this case 
the solutions of (1.1) can therefore be expressed in terms of analytic functions of a single 
complex variable. ; 
For positive or negative integer n the solutions of (1.1), (1.2) may be expressed 
similarly in terms of analytic functions of a complex variable by either one of two 
methods: the first consists in differentiating the solutions for n = 0 with respect to 
x, y, 2 and using inversion; the second, for positive integer n, utilizes Euler’s identity 


for homogeneous functions of degree n: 


oO. 0 4 2) 14) 

ra—+y-; z—le = ne. (1. 
ax | 4 ay dz/* ” 

The case n 1, also of interest from the point of view of conical flows, is discussed at 


length, and will be applied in the following paper. 

2. Harmonic functions of degree zero. In spherical coordinates the Laplace operator 
may be written in the form 
dy , 2 dy 


9 = Tm 3,(¢); (2. 
V'e = op? + Rar + B20): seit 
where B, is the second differential operator of Beltrami along ‘the sphere R = const. 


*Received April 26, 1948. 
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through the point in question. It will be recalled’ that for any surface S, the operator B, 
can be defined by means of 


— ; i dp / y. 9 
Bg) = lim | a ds/A (2.2) 


A-—0 


where, as shown on Fig. 1, the integration of the normal derivative is carried out over 





the boundary curve C of an area A on S, the normal to the curve being tangent to S, 
and the limit is taken as the area A shrinks to a point. The form (2.1) is obtained by 
applying Green’s theorem to a volume bounded by two spheres of radii R, R + AR 
and a conical surface through C, dividing by AR and allowing AR to approach zero. 

If y is homogeneous of degree n, then along t'x:e unit sphere R = 1, Eq. (2.1) reduces to 


B,(¢g) + n(n + le = 0. (2.3) 


In particular, for n = 0, we are led to 


over the unit sphere. 
° ’ ° ” . bd 
To solve (2.4) (along any surface S), introduce along S ‘‘isothermal’’ coordinates p, 
q, that is, coordinates such that 


ds* = X(p, q)[(dp)* + (dq)’]. (2.5) 
It is known’ that then the operator B, reduces to 
Bo(~) = (Ppp + a)/d.- (2.6) 


‘W. Blaschke, Vorlesungen ueber Differentialgeometrie, vol. 1, 2nd ed., J. Springer, Berlin, 1924, pp. 
114-116. 
*Jbid., pp. 121-122. 
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Thus Eq. (1.1) is reduced to 


Por + Po = 9, (2.7) 
whose general solution is 
¢ = f(p + iq) + g(p — 79) (2.8) 
or, if ¢ is real, 
¢ = Re [f(p + 19)] (2.9) 


where f, g are general analytic functions of their complex arguments. 
[t will further be recalled that if one system of isothermal coordinates p, q is available, 
the general system of such coordinates is given by 


p’ + iq’ = h(p + 7q) 08, 


where / is an arbitrary analytic function of the complex variable p + iq; this is in 
complete agreement with the fact that B, = 0 reduces to (2.7), and that under the 
change of variables (2.10), Eq. (2.7) is transformed into a similar Laplace equation in the 
coordinates p’, q’. 

Geometrically (2.5) can be interpreted to mean that the surface S is mapped con- 
formally on a plane in which p, q are rectangular Cartesian coordinates. The term 
“isothermal” is tied up with the fact that Eq. (2.4) is satisfied by the temperature in 
steady state conduction problems over a thin conducting shell of uniform thickness 
and thermal conductivity, provided no heat is lost over the faces. In particular, in each 
isothermal coordinate system, the curves p = const, g = const may be viewed as a family 
of isothermals and their normal trajectories. 

To render the solution obtained non-vacuous a particular set of isothermal coordinates 
is required. It will now be shown® that a set of isothermal coordinates on the sphere 


S:2°+y+2 =1 (2.11) 


is obtained by projecting a point P on S from the point Q: (0, 0, —1) of S on the negative 
z-axis onto P’ on the equatorial plane z = 0; the x, y coordinates of P’ in z = 0 con- 
stitute the parameters p, g of P (see Fig. 2). A projection of P on the tangent plane 
z = 1 at Q’: (0, 0, 1) is, of course, also possible, leading to values of p, g which are double 
those obtained by using the equatorial plane z = 0. 

Introducing in z = 0 the polar coordinates r = OP’ = tan y and w = tan ‘y/z, 
we have for the element of length de in that plane 


do® = dx’ + dy’ = dp’ + dq’ 


dr’ +r dw 


sec’ y dy’ + tan’ ¥ dw’. (2.12) 


On the sphere S, using the angle Q‘OP = 2y.= +6 and w as spherical coordinates, we 
obtain for the element of length 


*Jbid., p. 168. 
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P(x,y,2) 
2 
z 
Q(0,0-1) fe) Q! 
x 
Fia. 2. 
ds- = d(2y)’ a sin* 2y day 
= 4(dy° + sin” W cos’ W dw’). (2.13) 
Comparing (2.13) with (2.12) we obtain 
ds’ = 4 cos*y(dp* + dq’); (2.14) 
and since 
tany=r=p+q, (2.15) 


it will be concluded that. (2.14) is of the form (2.5), hence that p, g are isothermal co- 
ordinates on the unit sphere. 
The components of the vector QP on Fig. 2 are x, y, 1 + z where z, y, z are the co- 


ordinates of P; hence the components of the vector QP’ are 


‘ & Es 1 9 : 
i+s’ i+2’ ‘ (2.16) 


Thus the coordinates p, g of P’ are given by 
Y h 


while the complex argument of (2.9) becomes 


SS 





+- 
+ iq = ' — (2.18) 


x 
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From the above it now follows that the general (real) harmonic function which is 
homogeneous of degree 0 in 2, y, z takes on over the unit sphere the values which form 


the real (or imaginary) part of 


; ; xt wy 
{(p + 7q) = (z+) 2.1 
fet+iMm=AT5, (2.19) 
In order to obtain a representation for this function at arbitrary points of space, we 
replace x, y, z above by their ratios to R, thus obtaining 
p{ 2/ + ty 5) = (2+) - gz) (2.20) 
“\ 1+2/R "“\R+2z — F 
where 
‘ xt wy 
Z= — 2.21 
R+z ( ) 
It has thus been shown that any real harmonic function which is homogeneous of 


degree 0 is the real part of an analytic function of Z. 

As further examples of isothermal coordinates it is evident that an equally good set 
of such coordinates on the unit sphere can be obtained by projecting it from any one 
of its points onto the equatorial plane normal to the radius to that point. Thus if we 
project (2.11) from the point (0, —1, 0) onto the plane y = 0 we are led to the complex 
variable 

i «io (2.22) 
R+y 
in place of (2.17), as another possible variable in terms of which the harmonic functions 
in question may be expressed. The quantity Y must be expressible in terms of the 
variable Z in accordance with (2.10). Indeed, we may verify that 


. i—Z , ,~Y-1 
_ == —_—o—_-— ‘ 
) Tao? Z=(-i y (2.23) 
and this agrees with the fact that the only conformal mapping of a plane on itself (in- 


cluding the “point at infinity’ as a “point”’) is given by a linear fractional transforma- 
tion. Similarly we obtain from a projection from (—1, 0, 0) on x = 0 the variable 


_ 


eft StF he ; 
tt ewieke @ mies f (2.24) 
\s a further example consider (see Fig. 2): 
p’ + iq’ = log (p + ig) = logr + iw = log tan (6/2) + tw (2.25) 


where (6, v) are spherical coordinates on S. This isothermal coordinate system consists 
of small circles (circles of latitude) with Q as pole and meridians through Q. The function 


¢ = log tan 0/2 (2.26) 
corresponds to the temperature solution of a point source at Q’ and a point sink at Q. 


In space the harmonic function (2.26) corresponds tothe potential of a uniform positive 
charge distribution over the positive Z-axis, and a similar negative distribution over the 


negative z-axis. 
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The (p, qg)-net given by (2.17) consists of circles tangent to the directions of the x 
or y-axes at Q. The isothermal curves correspond to temperature flow due to a doublet 


at Q. 


3. Alternative derivation of the above results. Let » be a homogeneous harmonic 
function of degree zero. Then Euler’s identity (1.4) yields 


ro, + yg, + zg, = 0, (3.1) 
and when applied to ¢, , ¢, , g¢. which are of degree — 1, 
UPse TV YOu TF Bax = HPs» 
LOry + YO + 2s = — Gr; (3.2) 
LWOre TF YOus + 2: = ~“: - 


Solving the first two Eqs. (3.2) for ¢,, , g,, and substituting into the first en we 
can express ¢g,, in terms of x- and y- der rivatives of g as well as in terms of g, . Elimination 
of the latter by means of (3.1) and utilization of 


x0/dx + yd/dy = rd/dr (3.3) 
leads to 
2° Ges = Tre + 2re, ; (3.4) 


substitution into (1.1) and introduction of cylindrical coordinates r, w, then yields 
— wr r 2r l — 
Pez + Oy + (27g, + r¢,-)/e =\1+ a)¢,- + : r+ 3), + Pow = 0. (3.5) 


For z = 1 this reduces to 


(r* + Y)Orre + (2r° + T)o, + Pow = O (3.6) 


Ke (1 + 7°)” ai + par = 0. (3.7) 


By introducing the variable s 


or in a modified form, 





0= ape “ia ies 
we may reduce (3.7) to 
Pa t+ Poo = 9, (3.9) 
leading to the introduction of the analytic functions 
f(s + iw). (3.10) 
For z other than z = 1, Eq. (3.8) is replaced by 
$= log = + FT Ry = log | Z|, $9.11) 
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leading to 


s+ tw = log Z. (3.12) 


As a further alternative procedure‘ we may start with spherical coordinates R, 0, w® 


writing (2.1) in the form 


2 1 l $e) 1) s| 7 
Pre t pee + pe (a 630) * sin? 0 du | ~ (3.13) 


and for ¢ of degree zero, convert to (3.9) by introducing 








s = log tan 6/2; (3.14) 
its general solution is exprssible either in the form 
¢ = Relf(se'*)], (3.15) 


or in the form (3.10). 


4. Harmonic functions of integer degrees. The Euler identity (1.4), when applied to 
a homogeneous harmonic function ¢, of positive integral degree n, yields 


ae eer ee 

o> tor + Yay +2 az)” (4.1) 

thus expressing ¢ in terms of the harnonic functions ¢, , ¢, , ¢, of degree n — 1. By 
successive applications of (4.1) we are led to the identity: 
( 9 oa 9 y 

On = nl \" Or ty oy vs az/ °"? (4.2) 


where the subscript denotes the degree; this shows that ¢, is a polynomial in z, y, 2 
with coefficients which are harmonic functions of zero degree, so that each one is ex- 
pressible as the real part of an analytic function, say of Z, as in (2.20). These analytic 
functions, however, are not unrelated to each other. We shall illustrate this for n = 1. 

Suppose that ¢ is a real harmonic function of degree 1, then ¢, , g, , g, are harmonic 


of degree zero, and hence 





yg, = Re(U), yg, = Re(V), yg, = Re(W) (4.3) 
where U, V, W are analytic functions of Z. These functions must satisfy the three 
relations obtained by equating the cross derivatives: 

ae.) _ de) le) _ He.) lo) _ He.) “ 
oy ax’ dz dy ’ Ox dz ; 
Thus the last equation yields 
aw)\ _ (22) 
Re | = Re ae) (4.5) 


‘H. Bateman, Partial differential equations of mathematical physics, Cambridge University Press’ 
1932, pp. 356-357. 
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Considering U as a function of Z, we obtain 


ou) = (ae az), 
pe = RAO az)” 


now from Eq. (2.21) 


OZ r+ ty (3 ) Z 
——<— a- ~ Ps l — <—_ aa? 
dz (R+2)\R + R’ 
hence there follows 
au) = (mt z) 
Re 9 Re IZ R 
Similarly, considering W as a function of X and noting that 
OX =P x ay 7 r 
ox |: ia OY a R 
we obtain 
(aw) Ff ‘4 
Re =) = Re dX R/ 


Upon equating of (4.8), (4.10) and cancellation of R we obtain 
dW ‘) (wu /) 
v > = ft l 
RAS a ae p 
whence 
dV dU 


axe =a 2to 


where C is a real constant. Dropping the latter, there results 
TT = [ XdZ 


Z aX dW. 


Recalling (2.24) it will be noted that 
dX oe X) dZ dZ 2dZ 
— «== 6 of a j= = = : 7“ " r- = — o> 
X = 1-—-Z 1+Z (2 =— 1)’ 


hence (4.13) becomes 


Similarly (2.23) yields 
dY _ _2idZ 
y 2#+1’ 
and from the second equation (4.4) we now obtain 


pg A I cle I aces 
Vue Z dW. 
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(4.6) 


(4.7) 


(4.8) 


(4.9) 


(4.10) 


(4.11) 


(4-12) 


(4.14) 


(4.15) 


(4.16) 


(4.17) 
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The first relation (4.4) may now be verified from (4.14), (4.17). 
Application of (4.1) for n = 1 now leads to 


- Rel 2 [aw +2 [aw + aw | 








1 2. Z Zz 
(4.18) 
= s-— i ff ow, 28+ @ [68 . 
If the constant C in (4.12) and a similar constant of integration D occurring in the 


derivation of (4.17) had not been dropped, it would have been necessary to add terms 
—Cw, —Dw to u, v. However, the first relation (4.4) shows that C = D = 0. 

In a similar manner, by utilizing (4.1) for n = 2, expressing ¢,, , ¢:, , ¢:. in the form 
(4.18), then substituting in (4.4), it is possible to express ¢g, in terms of a single analytic 
function of Z and its proper double quadratures. 

Another method of handling homogeneous harmonic functions of integral degree n 
is by differentiation of functions of degree zero, and spherical inversion. We shall illus- 
trate this first for n = 1. 

Inversion of (2.20) leads to 





g-, = Rel[f(Z)/RI. (4.19) 
Upon differentiating (4.19) with respect to z, changing sign, and recalling (4.7), there 
results 
— Rel DZ 1D) | 
9-2 = Re| R + RR |: (4.20) 
Further inversion of (4.20) yields 
¢: = Re[RZf'(Z) + 2f(Z)] (4.21) 


for real harmonic functions of degree 1.* More generally, differentiating (4.19) n times 


iin «tt 2 : #2 |} (4.22) 
OZ 


for harmonics of negative integral degree, and by inversion of (4.22) we obtain a similar 


with respect to z yields 





expression for harmonics of positive integral degree. 
In particular, choosing 


f(Z) = 1 (4.23) 


the above procedure, in a familiar manner, leads to spherical harmonics with axial 
symmetry involving Legendre polynomials: 


3” 1 /yon+ 
P-(n+1) = an" (2) = const. P,(cos 6)/R"*’. (4.24) 


*To put (4.21) in a form resembling (4.18) replace R by z + (x — ty)Z. 
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Similarly, by choosing 
f(Z) = log Z, (4.25) 


we obtain axially symmetric spherical harmonics involving Legendre functions of the 


second kind: 
G-+1) = S E log +. | = const. Q,(cos 6)/R"*’. (4.26) 
Again, by letting | 
f(Z) = 2", (4.27) 
f(Z) = Z” log Z (4.28) 


we are led to harmonics involving the associated Legendre functions. 
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LINEARIZED COMPRESSIBLE FLOW* 


BY 
H. PORITSKY 
General Electric Company 


1. Introduction. The solutions of Laplace’s equation indicated in the preceding paper 
(referred to in the following as I) find ready application in the study of linearized 
compressible flow. This linearized approximate theory of high-velocity flow is based on 
the assumption that velocity components induced by an obstacle placed in a uniform 
stream of high velocity W, are small compared to W, . Under this assumption the 
velocity potential ¢ of the disturbed flow (that is, of the flow induced by the obstacle) 
can be shown to satisfy the differential equation: 


Prz + Puy + (1 _ M*)¢., - 0, (1.1) 


where subscripts denote derivatives, and M is the Mach number of the flow, that is, 
the velocity of the stream W, , assumed in the direction of the z-axis, divided by the 
velocity of sound at the pressure of the undisturbed flow. The same equation is also 
satisfied (under the assumption in question) by other thermodynamic gas quantities, 


for instance by the pressure p which is given by 
P — Do = —Wopow = —Wopoe.: ; (1.2) 


as well as by the velocity components u, v, w, given by 


u= ¢., v= g,, w= ¢.. (1.3) 
Equation (1.1) holds both for M@ < 1, that is when the flow is subsonic, as well as 
for .\4 > 1, when the flow is supersonic. In the former case, by introducing coordinates 
t,¥,2 = fe’, 8° = (1 — M”’) (1.4) 
one transforms the Eq. (1.1) into the Laplace equation 
bt ¢ = 0 1.5 
Cua TF Cur T Cee’ _— ° ( 5) 
On the other hand, to obtain such a transformation for the supersonic case, it is neces- 


sary to carry out the following change of variables: 

Z, yy, 2 = tz’, Bp? = M’ —-1 (1.6) 
in which real values of z correspond to pure imaginary values of the variable z’. The 
solutions of I thus have to be re-examined in view of this substitution. 

Of much recent interest are supersonic “‘conical’’ flows, that is, flows in which the 
velocity components u, v, w, the pressure p, etc. remain constant along straight lines 
through a point (the vertex or origin). Such flows are produced by conical obstacles, 
that is by solids whose boundary consists of straight lines through a vertex. To satisfy 
the assumptions of the linearized equation (1.1), the conical obstacle must lie close to 
the z-axis or it may be a thin fin containing the z-axis and inclined at a small angle of 
attack to the main flow. 


*Received April 26, 1948. 
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It will be recognized that for supersonic conical flows, u, v, w, p — Po are homo- 
geneous solutions of Eq. (1.5) of degree zero, while the velocity potential is a solution 
of Eq.. (1.5) of degree 1. 

The theory of conical flows is due essentially to A. Busemann,' where in the bibliog- 
raphy are given references to earlier papers. Among more recent studies of conical flows 
are these by H. J. Stewart,” R. M. Snow.* 

While most of the following is devoted to homogeneous solutions of Eq. (1.1) for 
M > 1 of degree zero, solutions of degree n are also discussed. Examples are given both 
of conical flow and conical flow modified by superposition of other solutions. 

The supersonic case M > 1 will be understood throughout the following. 


2. Conical supersonic flow. To adapt the solutions of I to the solutions of Eq. (1.1) 
for M > 1, in view of Eqs. (1.5), (1.6) the following replacements must be made: 
z—2z/B1, (2.1) 


R=(*?+y42)?-vW4+y -—2/8)? =([2 -—B(2e’+y)'7/m, (2.2) 


(x + wy) 


P zt wy , : 
Sah 78. gy —— FP Pee (2.3 
R+z 24+ —-Be ty)” 2.3) 
From Eq. I. (2.20) one obtains for the general solution of (1.5) the real part of 
it Bi(x + ty) | 
ZB) = | ‘er Kr a (2.4 
m ) d 2 + (2° - Br)” ) 
It will also be convenient to introduce the variable 
Z B(x + ty) : is in 
ee kt eS =t+ in = pe (2.5) 


. ae 2\7172 

a + [z (a +y)) 
and this can be used as the argument of f in place of Z; here &, » denote the real and 
imaginary parts of ¢; p its absolute value, w its argument. Thus, one is led to the solu- 


tions 


.. (2.6) 


Each velocity component u, v, w, for conical flows can be represented as the real part 
of an analytic function of «. 

Turning to the relations I. (4.4) between the velocity components u, v, w—these 
relations are equivalent to the vanishing of the curl of the velocity vector—we must 
keep in mind that while I. (4.15) yields a solution of the equation 

aU _ OW (2.7) 


Oz Ox’ 





1A, Busemann, [nfinitesimale kegelige Ueberschallstrémung, Schriften der deutschen Akademie der 
Luftfahrtforschung 7B, 105-122 (1943). 

2H. J. Stewart, The lift of a delta wing at supersonic speeds, Q. Appl. Math. 4, 246-254 (1946). 

3R. M. Snow, Aerodynamics of thin quadrilatera :wings at supersonic speeds, Q. Appl. Math. 5 


417-428 (1948). 
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for the complete functions U, W (that is, not only for their real parts) in the old co- 
ordinates, on account of (2.1), it must be replaced by 


pa Br Z = l oe 8 Pe +1 . ) 
u= 4% | 7 oT = 4 ain (2.8) 


new coordinates. Likewise I. (4.17) is replaced by 


Le) 
~] 
«+ 

— 


to furnish a solution of (5 


tA aie ok 2. ; 
v=5/ + an =F! —— ae (2.9) 


Taking real parts of (2.8), (2.9) to obtain u, v, one obtains 


ee 5 (| dW + | |, (2.10) 


€ 


where bars denote conjugates. 
For conical flow it is sufficient to determine the velocity components in a plane 


z = const., for instance in the plane z = 1. For z = 1, e€ reduces to 
B(x Bre*® Pee 
_ Be +o) = —s-5573 = pe”. (2.11) 
i+ {1 — B(x’ + y’)] 1+ (1 — 6’) 


The radical in Z or € vanishes along the Mach cone 


r -_ (x + y’) _ 2” 8”. (2.12) 
Its section by the plane z = 1 yields the “Mach circle” 
r= 1/8, z= 1 (2.13) 


Inside the Mach circle, the radicand in (2.11) is positive, and with the positive radical 
« may be considered to be a map of the Mach circle on the inside of the unit circle 


p | with radial direction angle w preserved, but with a radial distortion. 
A geometric representation of this transformation is shown in Fig. 1 where the 
point P in the plane z = 1 is transformed into the point P, in the «plane which has 


been placed on the z = 1 plane. This is done by drawing a sphere of radius 1/8 through 
the vertex or origin O, with center on the z-axis, projecting P by means of lines parallel 


to the z-axis on P,; and P{ on this sphere, then projecting these from O back on z = | 
to P, and P} respectively. Indeed, note from Fig. 1 that 


sin a = r/(1/8) = r8, 
(2.15) 
a _ P,0’ 


hoes = 0'0 = p. 
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(2.16) 
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one obtains from (2.11) and (2.15) 


rB 
— oe 2.1 
°-T4+ 0-378)” Oi) 


2.12), as well as the obvious fact that O, P, P{ , P, all lie in the same plane 


Kquation 4-14), 
through OO’, leads to (2.11). 

On the Mach circle, that is for r8 = 1, in the plane z = 1, the points P, , P{ coincide, 
and likewise P, , P3 coincide and lie on the unit circle p = 1. For P inside the unit circle, 
two poi P, , Pj are obtained in the eplane and these correspond to reciprocal values 
of p; the esult from (2.17) by using a positive and negative radical. 

If P lies outside the Mach circle, the construction of Fig. 1 fails; now the radicand 
in (2.11) is negative and one obtains 

bre 
é 7 8 mer; (2.18) 
1+ 267° — 1) 
\s will be noted from Fig. 2, now 
rB = csc 6, (2.19) 


where 26 is the angle subtended by the Mach circle r8 = 1 at the point P,in question, 
and Eq. (2.18) yields 
2m 


w ow 


ese 6¢ See (2.20) 


ws — ¢( 


1 +2 cot 6 et# 











Fia. 2. 
; : . 7 
Now w — 6, w + 6 are the angles \, u made with the real axis by the tangents from P 
. x ag ae ‘ 3 a 
to the Mach circle. Replacing the argument ¢ in (2.5) by its function log e/t =°A, un, 


we now obtain in place of (2.6) 


SO) + glu), A= wo — 4, p=wt+é (2.21) 
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for the velocity components of a conical flow, where f, g are arbitrary real functions; 
their arguments \, uw are constant along the tangent lines to the Mach circle. 

According to Busemann’ the introduction of ¢« is due to Chaplygin, whom he also 
credits with Fig. 1, and Eq. (2.10). 

Recalling the alternative procedure of I, Sec. 3, one obtains in a similar fashion 
for homogeneous zero-degree solutions of Eq. (1.1), or through the substitution of (2.1) 


in I. (3.5), the following differential equations: 


Br 1 26 r' p 
(1 — 3 Joo Bs ( — 2 Je. oa Fue = O (2.22) 
| \r Z r 
or in (2, y)-coordinates 
Bx 2p xy By” 28° (xy, + yy,) eer 
(1 —"$)o, — HF ¢., + (1 — =? }e — 2. =0. (2.23) 
It is sufficient to solve these for z = 1, and extend the solution to other values of 
z by replacing r by r/z. 
For z = 1, Eq. (2.23) is elliptic for r < 1/8, hyperbolic for r > 1/8. Its characteristics 
are given by 
(1 — 2°B’)dy’ + 26"dx dy + (1 — y’B’)dx” = 0, (2.24) 
whose integration leads to the tangent lines to the circle r = 1/8 of Fig. 2. 


For conical obstacles lying inside the Mach cone, the disturbance field vanishes 
outside the Mach circle. For conical obstacles outside the Mach cone the disturbance 
field in the plane z = 1 extends to the region enclosed by the characteristics from every 
point of the obstacle section and the Mach circle itself. Thus for a yawed conical fin 
shown in Fig. 3, the disturbance in z = const. covers the horizontally shaded region 
consisting of the Mach circle and the region between the circle and two of its tangents 


through the fin end. 





Fic. 3. 


$y introducing s = log p and w as variables one converts Eq. (2.22) for z = | into 


the Laplace equation as in I. (3.9) for r 8 < 1, while for r8 > 1, by X, uw as in (2.22) 
one is led to 
(2.25) 
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3. Flow around cone. The simplest example of a conical flow is obtained by choosing 
W = A loge = A(log p + ww) (3.1) 

where A is a real constant. The real part of this, 
W = A log p, (3.2) 


represents an axially symmetric flow: it will be shown that this flow corresponds to the 
disturbance produced by a circular cone of small angle placed with its axis along the 


z-axis. Indeed, Eq. (2.10) yields 


= — 8A ( — 1) -84(1_ Je 
utw= 9 \¢ 7” SG pk’; (3.3) 


the factor e’* shows that u, v lies in an axial plane and the factor p — 1/p shows that 
its magnitude depends on p only. 
If y is the (small) cone semi-angle, then the condition 


u+w 
epee tany = 7 (3.4) 
holds along the cone boundary r = tan y (for z = 1) along which Eq. (2.11) yields for 
small 8) 
Br 6B tan + By — 
y= — eo iar (3.5) 
1 + (1 — Br) 2 2 
Substituting this value in (3.4), (3.5) and neglecting p in comparison with 1/p, one 
obtains 
BA2 ; se — 
oa, = Woy, A= Woy (3.6) 
“VY 
and 3] \ ields 
W = Woy’ log «, w = Woy log p. (3.7) 
Substitution from (3.7), (3.8) in (1.2) yields for the pressure rise over the conical surface 
»— . B 2 | — 
P— Bo, = 24° log = 2° log (=), (3.8) 
1/2 poW 2 BY 
which essentially agrees with the familiar solutions of this (linearized) problem. 
[t will be noted from (3.1), (3.3) that along the Mach cone through the vertex the 
conditions 


3.9) 


u=t w= p— p = O0forrB =1,p=1 


hold. These conditions are satisfied by the conical flow around any conical obstacle 
lving entirely inside the Mach cone. 

The above example is often used to introduce the variables log p, @ in a purely geo- 
metric and fluid-dynamic manner, namely by means of the pressure field developed 


around a circular cone. 
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4. Lift decrease near edge of wing. As a further example, consider the lift decrease 
near the edges of a wing in a supersonic stream placed (with its span 
ay from its edges, the well-known two-dimensional flow is set up with a 


normally to the 


main flow 
Prandtl-Meye1 
The approximate solution for small angles of attack a is obtained from 


form of (1.1) (see for instance, Durand*) and yields equal and opposite pressure 


expansion on the suction side and a pressure shock on the pressure side. 


the two-dimen- 


sional 


A . . . 1 e 
changes Ap on both'sides, confined to dihedral Mach angles cot ‘8 to each side of 
the wing, as shown in Fig. 4, and given by 

$1) 


+ Ap = +p.Woa/B. 


Z-AXIS 








‘WING 


+ Ap 


Fic. 4. 
Near each edge the flow ceases te 
the wing edge at the origin it is found that the disturbance due to the edge OE is conical. 
If the wing edge OF lies inside the Mach cone, the disturbance is confined to a Mach 


» be two-dimensional. Placing the forward corner of 


cone which is tangent to the dihedral angles. 

We shall consider the case of a rectangular wing first (see Busemann'). For this case 
Fig. 5 shows the Mach cones of disturbance at each edge, and Fig. 6 the section by the 
plane z = 1. 

Proceeding far enough down stream one comes to the end of the wing and the flow 
is further disturbed, but this disturbance is only felt down stream and does not affect 
the conical flow in the Mach cone originating from each forward vertex. 

The boundary values for the pressure p — pp are given in Fig. 6 and are as follows: 


(—]1 over AB, 


mer = ++1 over AD, (4.2) 


0 over BCD. 


Similar conditions apply to w which is related to p — po through (1.2). It will be noted 





4‘W. F. Durand, Aerodynamic theory, Vol. 3, Springer, Berlin, 1935, p. 235. 
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that (3.9) is violated, which is not surprising since the wing does not lie inside the 


Mach cone. 

In addition to (4.2) there is a condition over AO. To determine the latter note the 
boundary condition 
(4.3) 


v= —W, tana = —wW,a, 


which holds to each side of the wing surface, and which we apply over its projection 


on y = 0. Differentiation of (4.3) yields 

















Ov 
— =), 4.4 
Oz (4.4) 
y 
x 
y 
| 
Ip 
O E z t 
| 
Ap + , 
p> 
B | 
! a a oe 
0 
Ap 
D 
Fig. 5. Fia. 6. 
From the conditions of vanishing velocity curl, one now obtains from (4.4) 
ow - 
—— 2 (0, (4.5) 
oy 
This boundary condition enables one to continue w (and hence p — po which is pro- 


portional to w) from one side of the wing section AO to the other side by means of 
flection across AB. It must be kept in mind that the reflection refers to the 
analytic continuation of w across the boundary, and not to the flow actually existing 
on the other side. The condition (4.4) and the reflection is applied separately to the flow 


posit re 


above and below AO. 

Carrying out a positive reflection for w and p from above AO, and similar reflection 
from below, it is found that the values obtained for both images are the negatives of 
the actual values obtaining in the physical flow at the same points. Thus w, p are single- 
valued over a “two-sheeted Riemann surface’ which has a second-order branch point 
at O. The relation 
va (4.6) 


€, = € 
maps the two-sheeted interior of the unit circle of the «plane into the single-sheeted 
interior of the unit circle in the ¢;-plane with OABCDA of Fig. 6 going into the semi- 
circle OABCDA of Fig. 7. In the ¢,-plane the boundary values for (p — po)/Ap 
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are as follows 


| Oover DCB, 


a | —1 over BAB’, 
en = 3 (4.7) 


Ap ) 
I | Oover B’D’, 
| lover D’/AD. 
Discontinuities for w, p occur at B, e¢, = 7‘? and at D, ¢, = —i(7)'”’ as well as at their 
images in AOA, B’, e€, = i(i)'”*; D’, «, = —7'”*. In terms of ¢, we may set up (p — po) 


as follows: 


n= 2 
Ap 








where 6, , 62, -:: as shown in Fig. 9, are the arguments of the factors e, — 7' “+++ wis 


obtained by dividing p — po by —Wopo . In particular along AA, 6; simplify; thus 


p — 1/2 
6, = tan 9172 = tan” [(2p)'”* — 1], (4.9) 
and one obtains from (4.8) for the local lift coefficient along OA of Fig. 5, as a ratio to 
the lift along the wing proper, 


7 1 (2p) =< , 2p — : 
tan = —sin — : (4.10) 
T : =% T l + p 
Utilizing the relation 
9 
2p 
sr = 5 ( 1. . 1) 
l+p 


which follows from (2.11), one may replace (4.10) by 


_ \1/2 shes , | 
sin~’ (6r)'’? = = cos’ (1 — 6r)'”” = — cos‘ (1 — 26r). (4.12) 
: us T 


us 
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eplacing r by —z and averaging from z = 0 tox = —1/8, one obtains 


ef os”? =2 ff cos aa 
a Oe cos (1 + az) dz = 2 J, du = 5° (4.13) 
The net decrease in lift over the area within the Mach cone due to the edge effect is 


thus 50%. 
We turn now to a trapezoidal wing (see Fig. 8) whose edges make an angle r/2 + y 














x 
a Ww 
’ a oe 
” ol ‘i 
a “ 
~ 
a5 
y 
-| 
cot |g 
—_—— 
z 
Fia. 8.* 
with the forward edge, where 


0<y < cot 's. (4.14) 


is now replaced by Fig. 9 and in the e-plane the Mach circle is transformed 























Figure ( 
into the unit cirele of Fig. 10. «plane 
_ B B 
A 
B 
8 
cna 4 fC) 
0 E 
D 
Fic. 10. 


Fia. 9.* 


with the point Z going into « = a where 





*Nole added in proof: 6 in Figs. 8 and 9 should be replaced by v. 
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6 tan y . 
68 4.15 
1 — (1 — 6’ tan’ y)'” (4.15) 
The boundary conditions along p = | are still given by (4.7), while (4.5) now obtains 
our AOE. By means of 
= tS (4.16 
o._ : 
de — | ») 


one transforms the unit circle in the «plane into a unit circle in the ¢’-plane (see Fig. 11) 








with the slit AOE going into the radius along the negative real axis. The point B, « = 7 
of Fig. 9 goes into 


F a-?1 — 4 - 
a= —_—— = —74 -=e@e 
aa — |] a+ 
(4.17) 
yea —25 — «/2, 56 = cot’ a. 


Analytic continuation by means of positive reflection across AE of Figs. 10, 11 is possible 
and similarly leads to a single valued function in the 


€ = ” te (4.18) 





oO, 


€,-PLANE 





D’ D 


Fic. 12. 


points B, B’, D, D’ replaced by points making angles +v/2 with the directions of positive 
and negative e, . A solution analogous to (4.8) is now given by 


Fe. piv/2ye, pa tv/2 
os = ot Moar] 4.19) 





ae 0 l } 
P— Po _ - (6, — 6, — 6, + 6) = Re| - ‘ins rw =, 
Ta tee, +e”) 


Ap T 
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The solution just obtained (for p — p> and hence also for w), while it satisfies all the 
boundary conditions, will not do, however, for the following reason: since the wing 
possesses lift and builds up a pressure rise and fall Ap on the lower and upper sides 
respectively, there is a tendency to set up circulation around each edge, similar to the 
circulation of an incompressible fluid around the edges of an airfoil of finite span, and 
responsible for induced drag. Since the component of the main flow normal to the wing 
edge, W, sin y, is subsonic when (4.14) holds, this flow around the edge (indicated by 
arrows on Fig. 9) possesses a singularity at the edge similar to that of an incompressible 
fluid, leading to infinities of the form 


A(e — a)” (4.20) 
in the functions U, V, W. This type of singularity arises also in the case y = 0, but 
there it affects only u and v and not w*. For y > 0, the flow lines due to this singular 


flow lie essentially in planes normal to the edge, and W is also affected by it. 
To introduce the required singularity one adds to (4.18) 


rel | (2 _ «) | = c(t _ p,) sin w, : (4.21) 
1 €) Pi / 


This vanishes for p,; = 1 and hence does not interfere with (4.7); it also satisfies (4.5) 
over AOE of Fig. 9. Thus C in (4.20) cannot-be determined from considerations of p — po , 


T 
I+ 


w only. It may be found, however, by means of the relation 


aw 


dé |.=0 


= 
—_) 


(4.22) 


without which (4.19), (2.8), and (2.9) would yield singularities for U and W at « = 0. 
One is lead to 
4a a — cos (v/2 : 
C= - >> = 28 6/3) (4.23) 
1—aa + 2a cosrv+ |! 
This completes the determination of p — pp and w. 

For negative y with the wing edge still in the Mach cone, a Kutta-Joukowsky condi- 
tion holds at the edge and the solution (4.19) has to be properly modified. No singularity 
is introduced at the edge. 

—_ ~ 3 

The case y > cot 'B has been treated by Snow’. 


5. Non-conical flows. Harmonic functions of various degrees of homogeneity were 
considered in I, Sec. 4. By means of the substitutions (1.6), (2.1)—(2.3) they can lead to 
solutions of (1.1) of similar degree of homogeneity. 

As an example, familiar harmonic functions of degree n, — (n + 1) of the “product”’ 
type are given by 


[R", R~“"*” |P™%(cos 6)[cos m w, sin mw], (5.1) 


where P™ are the “associated”? Legendre functions, and cos 6 = 2/R. Carrying out 


*It is introduced by the integrations (2.8)—(2.10). 
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~ 


the indicated substitutions in (5.1), we obtain the following solutions of the wave 


equation (1.1) for the supersonic case (see Hayes’): 
2 2 2\ [n/2,—(n+1)/2] m 2 2.2 2 ° / 
(2? — pr’) ’ P'™(z/{z" — B’r’}'”*)(cos mw, sin mw); (5.2) 


additional solutions of the same product form may be obtained by utilizing the asso- 
ciated Legendre functions of the ‘second kind,”’ Q’", in place of P®. In particular, m = 0 
yields the axially symmetric solutions: 

2 2. 2\ ([n/2,—(n+1 2) 2.2 2\- 2 - 
(2” — p’r’) P[(l — pr’ /2)-*”"). (5.3) 
The velocity potential for the flow around an axially symmetric body with axis along 
the z-axis possesses this symmetry. 

As a special case for » = 0 and the second exponent in the first factor in (5.3), there 


results 


(2? =: 87?) 2 (5.4) 


This is the potential of a “point source’ at the origin; it is of degree —1. The velocity 
components for the same symmetric flow, as well as part of the velocity potential due 
to transverse flow around an axially symmetric body, possess the type of symmetry 
obtained from (5.3) for m = 1. 

General solutions of (1.1) of degree n are obtained from I, (4.21), (4.22) by intro- 
ducing e and the substitutions (2.1) to (2.3). 

In particular I.(4.21) vields 
sai) , ., — 

- de + zf(e) (5.5) 

for a general solution of (1.1) of degree 1. The real part of this could be used as a starting 
point for the velocity potential of conical flows. The requisite boundary condition is 
given by 


A, T A, + y.A T W A; = 0, 0.60) 


where (A, , A, , A.) are the direction cosines of the normal n to the conical boundary. 


Neglecting ¢, = w in comparison with W, , one may transform the above into 


eM tem = (Ext ¥u)., 


where (A! , Aj , 0) are the direction cosines of the normal n’ to the boundary section by 
9 


a plane z = const., this normal lying in that plane (see Fig. 13). Equation (5.7) can 


be put in the form’ 


®‘W. D. Hayes, Linearized supersonic flows with axial symmetry, Q. Appl. Math. 4, 255-261 (1946). 
*In principle, the singularity of the flow field of Section 8 at A, A’ could be obtained by replacing 
the wing with its sharp edges by a boundary of finite curvature, solving flow subject to (5.8), then passing 


to the limit as the original wing with its sharp edge is restored. 
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Fia. 13. 


a W, a(r°) Vo, ; 
SP = 2S = —* [x cos (n’x) + ysin (n’2)], (5.8) 
on 2z On zZ 
where (n, x) is the angle between the normal n’ and the x-axis. 
By superposition of conical flows and of solutions of other degrees of homogeneity, 
flows around a great variety of profiles may be obtained. 
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ON THE CHARACTERIZATION OF FIELDS OF DIABATIC FLOW* 


BY 
B. L. HICKS 


Ballistic Research Laboratories, Aberdeen Proving Ground 


INTRODUCTION 


In an earlier paper (Reference 1) the steady flow of a compressible fluid containing 
distributed heat sources was described in terms of the velocity vector V and alternately 
two other vectors W and M. For rotational flow the W language was found to lead to 
more simple differential equations than the other two languages. The physical content 
of the equations expressed in each of the three languages was, of course, the same. When 
each of the vector fields V, W, M, was in turn assumed to be irrotational, however, the 
corresponding fields of diabatic flow were radically different in their physical charac- 
teristics. 

These results suggest several questions concerning the characterization of fields of 
steady diabatic flow: 

(7) is W the most:'convenient vector to use in formulating the equations for rotational 
flow? 

(77) how do the restrictions upon the heat source function Q depend upon the char- 
acter of the vector function N chosen to represent an irrotational field of flow? 

(777) can a relation be established between the character of a vector representation 
N of an irrotational flow field and the form of the partial differential equation for the 
potential function gy ? 

We answer these questions by introducing into the flow equations a vector N pro- 
portional to V, the scalar proportionality factor (V/N) not at first being explicitly 
restricted. From the form alone of the derived equations in N, it is possible to answer 
the first question. The assumption of irrotationality of the N field, V xX N = 0, then 
leads to the required information about irrotational fields of flow. 

Through proper choice of the V, N relation and subsequent inspection of the form 
of the N equations, a number of new types of irrotational diabatic flow can be dis- 
cussed by giving their character, specified by the function g( VN) = V/( N’?RT)'” and its 
form, specified by an arbitrary function /(gy) that enters the partial differential equation 
for the potential function gy . Four general types of flow are chosen for discussion to 
illustrate the wide variety of diabatic irrotational flows, as compared to adiabatic 
irrotational flows, that are possible even without explicit formulation of boundary 
conditions. We thereby emphasize the need for preliminary and simultaneous investi- 
gation of all steady diabatic flows at the formulational level before the attempt is made 
to construct special solutions. 


*Received April 5, 1948. Summaries of this paper and of related studies of diabatic flow have been 
presented at meetings of the American Physical Society. [Cf. Phys. Rev. 71, 476(A), 72, 179(A) (1947), 
73, 636 (L), 74, 1230 (A) (1948).] 
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BASIC EQUATIONS 


1. Derivation of N equations. We use essentially the same fundamental equations as 
in Reference 1, namely those governing the steady motion of an inviscid compressible 
fluid that contains heat sources. 


; Vp + V-VV=0, (1.1) 

V-pV =0, (1.2) 

p — RoeT = 0, (1.3) 

coV-VT, = TV-VS = Te,V-V log (Tp"‘7~"’") = Q. (1.4) 


The quantity Q is the heat added locally per unit mass of fluid and unit time. In Refer- 
ence 1 these equations were studied in the V, W, and M languages, W and M being 
related to V by the equations 


V = V.W = aM (1.5) 
in which V, = (2c,7,)'” is the limiting velocity corresponding to local stagnation 
temperature 7, and a = (yRT)'” will be called the local velocity of sound.* In the 


present paper we make the more general substitution 

V = (gRT)'N (1.6) 
which includes the W and M transformations when g is appropriately specialized. For 
the present g can be any function of the coordinates and other variables and will not 


be restricted until later. Both g and N = | N | are dimensionless. 
Substitution of (1.6) in (1.1, 2, 4) and elimination of V and p gives 


V logp + gN-VN + gNN-VY log gT = 0, (1.7) 
V-N+N-¥V logp +5 N-V log (g, T) = 0, (1.8) 
—% —t N-V logp + N-V log T = Q/[c,T(gRT)'”’] = 2qy . (1.9) 


y 


In (1.9) we have defined a quantity gy which represents in dimensionless form the local 
heating of the fluid. We adjoin the equation obtained from (1.7) by scalar multiplication 


of N 


N-V logp + * N°gN-V log T = — * N°gN-V log g — gN-V7 E N’ (1.10) 


— — 


*In a gas undergoing dissociation or chemical reactions, the velocity of sound is a function of the 
local thermochemical properties of the gas. 
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and solve (1.9, 10) for N-V log p, N-V log T in terms of N-V (N’/2), N-VY log g 


and q, 


N-V logp = -(N-v 3 N* +3 N’N-Y logg +N ‘aw (1 +7 N*9), (1.11) 


N-V log T 





- | 20. - ar (N-v 4 vin? +5 l AN. 7 log a) |r +2 aad ‘a) (1.12) 


Substitution in (1.7,8) gives the equations of motion and continuity in N language. 


7 ly — 1 -_ 1 
V log p = AC += N ‘a) (St on-v da — SN. log 9 — ay) 
(1.13) 


— 9V N? + gN X(V XN), 


rab 
7-N = (1+ Dy N“g 5 \l - “ N-VY log g 


(1.14) 
+ 1 2 
“> oN-Y 5 74 (1 + N’g)qy |. 
For completeness we list the expressions for N- VY log T, , N- V log p, , quantities pro- 
portional to the rate of variation of stagnation temperature and of stagnation pres- 
sure in the direction of flow. 


«gp SD ws o 3 . 
1+ a V%g)N-V log 7, al aie thee = 2qy , (1.15) 
“yy Pp 
— a —_-_ 2 
N-VY log p, = —N*gqx/\1 + 2 N°g|}. (1.16) 
As we might expect, these variations are proportional to gy and vanish for gy = 0. 


Furthermore 7’, increases and p, decreases in the direction of flow for positive gy . The 
decrease of p, owing to heat addition is known as “momentum pressure drop” in one- 
dimensional theory (compare Eqs. (15) and (60) in Reference 2 and references cited 


there 


2. Significance of choice of the function g. We can now observe how the character 
of the function g affects the es of the N equation. The factor [1 + (y — 1)N’g/2y] 
appears in both Eqs. (1.18, 14) together with g, Vg and qy . If either g or qy is an ex- 
plicit function of 7 (ef. the special case of V language described in Sec. 3) then it may 
not be especially profitable to reduce the equations to the N form because 7 cannot 
be determined and eliminated from (1.13, 14) prior to complete solution of the set 
(1.1, 2, 3, 4) or the set (1.12, 13, 14, 15). If, on the other hand, g and qy are functions 
of N, p and the coordinates x; or even of the individual components N; of N together 
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with p and the coordinates, one can regard (1.13, 14) as a set of four simultaneous 
non-linear partial differential equations in the four dependent variables N, p. 
We choose to study a less. general function g, namely 


g = g(N) (2.1) 

corresponding to the form for qy 
Qn = Qn(N; , 2;). (2.2) 
When g and qy are thus restricted, the N equation of motion (1.13) contains as dependent 
variables only N and p and the N equation of continuity (1.14) only N. Equation (1.14) 


can then be rewritten in a more compact form which exhibits the “sources” proportional 
to gN in the gN/H, field, letting g’ = dg/dN, 


— -1 

V-H,N = (1 + v*o(1 7 N*g) Ries (2.3) 

where 

ae a - Na) ¢ | Sal | *9) ly 9g 
log H,(N) = + | ~ g (1 me a 1+ Dy N°g daa : (2.4) 
Equation (2.1) includes both W and M languages for which g has the forms 
2 — : . 2 
gw =—C(1- Ww’), (N=W), (2.5) 
eae | 

ju = 7; (N = M). 2.6) 
Since N*g = V’/RT has the same value for all g(N), the relations among V, V, W and 

M are 
+2 _ r2 2y 72 r2 1 2 my 
(V*/RT) = N’g = ——— w*(1 — W’)"' = y¥M (2.7) 

y=. 7 


These equations are often used in the form 
a 72\—1 : eed 2 9 2’ 
1 + + —N’g) = (1 — W?)"' = (1 + —’). (2.8) 

\ 2y 2 
It is only W language that a term in N-V N(or V-N) does not appear in the equa- 
tion of motion. The W language therefore yields more simple equations than any other 
N (or V) language, and W is thus the most convenient vector to use in formulating the 


equations for rotational flow. 


PROPERTIES OF IRROTATIONAL N FIELDS 


3. Integrability of the equation of motion. We now suppose that the N field is con- 
tinuous and irrotational and therefore admits a potential function gy 


The consequences of the assumption of irrotationality for an N field were first examined 





0 ES 


a 


ee 
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in Referer 3a, b for the adiabatic M and W cases. In Reference 1 we saw that irrota- 
tional « M and W flows were quite different from one anothe:. We may therefore 
now ex} . still wider variety of flows because of flexibility in choice of the function g. 
Witl < N 0, the equation of motion becomes 
cy — 1 a2 ; 
V logp + nVev + 9V5N =9 (3.2) 
in i 
] 
= gN-T7 log (qT) 
a (3.3) 
( y¥=~—1.4\" — 1 er 
gh 1 +> —N’g) | wv — G—ogN-V ZN? +5N-V log g |. 
2 - 2 2 2 
\Ve distinguish two cases, 7 = 0 and ny # 0. For the degenerate case 7 = 0 we suppose 
iat g isn estricted to be a function of N alone. Then from (3.2, 3) 
-_ _ 
V logp + 9gV5N =0, (3.4) 
ek ee Pee ee — 
>, 9N-V SN —5N-V logg = w- (3.5) 
2) 2 2 
Equation (3.4) implies that p and g are functions of N alone. Also Eq. (3.3) shows that 
gT) and therefore, by (1.6), that (V/N) are constant on each streamline. This type 
of flow therefore corresponds to irrotational V flow (see Reference 1) that has been 
generalized to allow variation of (g7) between streamlines. 
When 7» + 0, it is sufficient to consider only functions of N alone as was shown to 
be appropriate in See. 2. If the integrability condition with g’ = dg/dN 
l - y y q’ y ae l ( » Bivens ] 72 ee 7 ‘ . 
n = 35-9| 2@v + NN-VAN here i+“ I's = F(yy) (8.6) 
2°L Ng Y \ ZY / 
for the equation of motion (3.2) is satisfied then the latter integrates to 
log Py + | F(¢gx) dgy = constant, (3.7) 
wher 
Py = pH.(N) = p exp [ gd N’. (3.8) 


Integration of (3.3) also yields a relation among N, p and T along each streamline.) 
Equation (3.7) is the analogue of the barotropic condition p = f(p) for irrotational V 
fields of flow. It states that Py is constant on potential surfaces and expresses the varia- 
tion of Py along streamlines, i.e. 8 log Py/ds = NF(gy). The function Py becomes 
equal to the stagnation pressure Pw = p, in W language. In general the behavior of 
Py is similar to p, in that as N (or W) increases, p decreases for constant Py or p; . 
It is noted that in terms of H, , the continuity equation (1.14) becomes (compare (2.3)) 


V-(gN/H.) = Hz‘(1 + N’g)F(¢) (3.9) 


showing that “sources” in the (gN/H,) field are proportional to F(gy). 
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The integrability condition can be rewritten in two other useful forms by com- 
bining it with the continuity equation (1.14). Then 


afy+1 wf _ N°) ¢ | , 


(3.10) 
= (£ -1=1 pvp, (,—£) 
= (¥ y g V YN + 2 MY] Na qn 
and 
N’ 1 Ng’\, 
qJ (1 =~ (1 + 5 oy (yy) 
7 s g 
(3.11) 


tele ¢=—} ) ! a A( = ) | 
“= a 7 ert Se a 


in which s = N/N is a unit vector in the direction of flow. According to Reference 3a, 


4 1 OAA 
a = omnia poh n ates 
vs lim AA 0s 


4A-—0 


can be interpreted as the fractional rate of change of stream tube area AA with respect to 
are length ds along the stream tube. Equations (3.6, 10, 11) are alternative statements of 
the relation between the character of the function g(N) and the mode of variation of 
Y -s and qy throughout the field of flow. Since the coefficients in these equations depend 
only upon N and g(N), through special choice of the function g(N) we will later be 
enabled to study flows characterized by special relations among the functions qy , 
V-s and F(gy). 

4. The partial differential equation satisfied by the potential function ¢, . li N = 
Vew is substituted directly into (1.14), the resulting equation contains gy as well as 
derivatives of gy . Because in general the function gy is not open to arbitrary specification 
in an iirotational N field, it is desirable to consider (gy) that appears in the integrability 
condition as the basic function describing the form of the flow. For a given g(N) and 
F(gy) the behavior of gy follows from the integrability condition. It is therefore best 
3.9) which does not contain gy and which can be expanded to read 


to use Eq. (3 
i 


+2 : ra | 
7-N=g (1+ Ng)F(en) + (g me JN-V, N°. 1.1) 
\ Nq 2 


Substitute 


’ — Oy oe —~ (Ipy\” Ayr dgy Igy Igy 


\ 02, 


The desired differential equation for the potential function is then 


—~ O On [ Oex ‘ Og, Igy Og = = 
hes. rt 8 oP ae a od Ge od 
2B OX: | ‘ Ox, i 2 Ox; OX; OX; OX; q is 


OX 
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in which 


, 


y g 
G=g- =. 4.3 
ont (4.3) 
The type of this partial differential equation is determined by the sign of the de- 
terminant D = det | a;; | where a;; is the coefficient of d°¢y/dx,0x; in the partial differ- 
ential equation. Calculation of D gives 
q’ 
D=1-NG=1+4N 7 — Ng. (4.4) 
tf 


The differential equation (4.2) is therefore of 


elliptic (G < N7?, 
parabolic typefor 4G =N™, (4.5) 
hyperbolic) iG > N77. 


Since G is a function of N alone, the equation can change type within the field of flow 
where NV attains the values specified in (4.5). This behavior is similar to that of irro- 
tational adiabatic flow in the V field which is of elliptic, parabolic or hyperbolic type 
for M <, =, >1. There is however, a major difference. The character of the function 
g(N) alone determines the values of N, if there be any, for which each of the expressions 
in (4.5) is satisfied. Appropriate choice of g(N) can then lead to flows whose type is 
fixed (e.g. elliptic for all values of V)) as well as to flows which change type for specified 


values of N. Because the character of g likewise determines the mutual behavior of 


V-s, gy and F(gy) (ef. Eq. (3.11)), the required connection between character of g 
and behavior of gy has thereby been demonstrated. Whenever, for example, D = 0, 
the relation between V-s and (gy) does not depend upon qy (ef. Eq. 3.11). In a later 
section, we shall study some flows illustrating these remarks. 

There may be a question at first sight as to whether F(gy) and gy can be found such 
that ¢y will satisfy both (3.10) and (4.2). We can analyze the situation as follows. Let 
N V,], 7 = 1, 2, 3. For rotational flow there are five dependent variables N; , p 
and gy that satisfy the four equations (1.13, 14) giving an underdetermined system. 
The variation of one dependent variable, e.g. p, gv (or one N,) can then be specified 
= dgy/ dx; we have seven equations in the six variables N; , 


arbitrarily. Now when JN, 
P, Gy and gy and the system is overdetermined. The equation of motion (3.2) however 
reduces to the two equations (3.6, 7) which contain the new (arbitrary) function F(¢gy). 
As this function is not determinable from the basic equation we regard its specification 
as part of the formulation of our problem and can consider then (3.7) to be an integral 
of the motion which eliminates p from further consideration. We now take (3.10) and 
(1.14 or 4.2) to form the basic system that must be solved. This system contains qy , 
Fy , gy and its derivatives and is accordingly underdetermined. If we choose gy to be 
the known function (1.14) with (V = Vey) can presumably be solved for gy . There 
would then be, however, no reason to suppose that these functions gy and gy would 
satisfy (3.10). We shall consequently regard (4.2) as the basic equation. Choice of the 
functions g and F (gy) will specify the character and form of the flow respectively. Pre- 
sumably specification of both form and character will work a corresponding restriction 
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“ 


upon the nature of the boundary conditions that can be imposed on the function gy 
For example we have already seen that the type and change of type (with corresponding 
adjustment of the nature of the appropriate boundary conditions) of the partial differ- 
ential equation (4.2) is entirely governed by the character of the flow 

The solution of a problem in diabatic irrotational flow can then be visualized as 
consisting of the following steps. With character and form of the flow specified together 
with matching boundary conditions, Eq. (4.2) can be solved for gy . Substitution of gy 
thus determined in the integrability condition (3.6 or 10) and the “Bernoulli” type 
equation (3.7) yields the functions gy and p. The system consisting of Eqs. (1.13, 14) 
‘an then be regarded as solved. Determination of the temperature follows by integration 
of (1.12) if T is known at one point on each streamline. The velocity V and heat function 

( 


Q can also be calculated subsequently from (1.6, 


The indirect nature of this solution process must be noted. Although one might 


normally expect Q to be known rather than F(¢y), it is the latter that, at present at least 
must be chosen first and Q thereafter determined. It is this peculiarity of irrotational 
diabatic flow which leads us to investigate in general fashion irrotational flows whose 
character only is specified, in order that we may derive more explicit equations con- 


taining q) 


THE RELATIONS AMONG CHARACTER, FORM, TYPE CHANGE AND HEAT 


SOURCES 
5. Methods of investigation. We now propose to illustrate the remarks in the last 
section by examining a number of flows of given character. There are several paths 


we may follow in proceeding to this specialization. We could for example try to find 
all functions for which the flow remains of a given type. This procedure is simple only 
for cases for which D is constant throughout the flow. We might also return to Eq 
(3.11), which connects the physically visualizable quantities V-s and qy , and impose 
conditions on the coefficients appearing there in order to insure, for example, that V-s 
be proportional to F(gy) throughout the field of flow. As a third possibility we could 
assume a priori convenient special forms for g(N) (such as rational functions) and 
examine the corresponding characteristics of the flow. Any of these procedures will lead 
to conditions upon g(N) and will supply partial answers to the questions implied by the 
other procedures. 

For most of the special flows to be discussed now the first procedure is convenient. 


Let us suppose throughout the flow 


D = D, constant. (5.1) 


Accordingly as D, is >, = or < 0 we shall term the flow field elliptic, parabolic or hyper- 
bolic. By (4.4) we are then led to the differential equation for g(N) 


or 
bo 
— 


N ; —N’%og+1-D,=0, (8. 


whose solution for D, # —1 is 


g = k-?(1 + D,)[(N/k)*~”? — (N/k)*7J"' (5.3) 
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and for D, = —1 
g = N-*(log k — log N)™’. (5.4) 


In both cases, (and also for D not a constant) 


1 =. Dy 
G= — yr (5.5) 


the partial differential equation for gy becomes 


PaaS oie doy \’ doy Opy doy 
> 2 O¢ N E we ( —_ iB (22x) | aan 1 — R On VEN a 
R U Dy) x Dy) p> Ox; OX; OX; OX; 


— 60%; Ox; 
(5.6) 
= ve + L\Ptex) 
g 
and the integrability condition in the form (3.11) simplifies to 
{ N’q 72 - , N’g ™ 
g \1 — —K1l + Do + N’g)F(ew) = 2Dogn + N\1 — Do — “aA V-s. (5.7) 
7 


The flows characterized by (5.6) do not resemble “incompressible” adiabatic irrota- 
tional flows for small N unless Dp, = 1 (ef. Sec. 6) because all terms in the coefficient of 
0° yy /Ax,0x; are of order N*. The Glauert-Prandtl type of approximation can be made, 
however, for suppose that throughout the flow 

N, N2 N; 


~ ] <3 mw 2 mw, 


N, , No No 


where JN is taken to differ but little from N, anywhere. Then (5.6) becomes 





p, 224 4 Hey 4 Fey. [nt + 1, [ron (5.8) 


dx; dx: x3 g(No) 


That the value of D, determines the type of the differential equation is now quite clear. 
Equation (5.7) may be developed in similar approximate fashion once the value of Do 
is chosen, and the function H, (Eq. (3.8)) can be regarded as a function of No . 

For some flows of interest (such as W and M flows) the type of the partial differential 
equation (4.2) may change in the flow field. A simple application of the second pro- 


cedure is then appropriate. Let 


g=g9+goN+g.N'+-:::. (5.9) 
Then for N < 1 
' , gi + 29.N 
G = g N- = ree, 5.10 
fe + ON — NG + gN) * iis 
Only for go # 0 will Eq. (4.2) admit an “incompressible” approximation for N < 1, 


The Glauert-Prandtl approximation can be written down for almost any choice of g(N). 
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6. Elliptic flow with D, = 1. The general formula (5.3) for g corresponding to flow 

of constant type reduces for Dp = 1 to G = 0 and 
g = H, = 2/(k’ — N”). (6.1) 

The equation for gy (compare (2.3) and (3.10)) is 


l 


V'¢ev = 9 (k* + N*)F(¢w) (6.2) 

which in the Glauert-Prandtl approximation (cf. (5.8)) is 
2 - r2\ 77 5 
Vier =35 (k° + No)F (ey). (6.3) 


- 


The heating parameter gy is given by 


qv = H(t = << nv") — N*)"F@y) + : N*(ke — N*)'V-s. (6.4) 


- 


The “‘Bernoulli” equation becomes 


log (2) + | F(¢y) dgy = constant. (6.5) 


Since the relation between W and N is 
W? = (y — 1)N’*/(yk* — N’); (6.6) 
values of W* = 0, (y — 1)/(y + 1), 1 correspond to 0, y/(y + 2), 1 for (N/k)’. 
It is noted that flow discontinuities such as stationary shocks and combustion fronts 
cannot occur within this or any other field wholly of elliptic type. 
7. Parabolic flow: D, = 0. The only function g that gives wholly parabolic flow is 


g= 1/N(k — N) (7.1) 


obtained from (5.3) where D, = 0. In this case G = 1/N’, H, (k—wN 
partial differential equation for gy becomes 


72 Oon F don Oey Opn dyn r3n = —) 
ies (Sex —3;-2 eo a ee oe EN Ploy (7.2) 
dX | Ox; Ox; > Ox; OX; OX; OX; ex) 
of which the Glauert-Prandtl approximation is 
ay I py el 4 
fy 4 <2 = kNiF (en). (7.3) 
OX 0x3 


When D = D, = 0 the coefficient of gy in Eq. (3.11) vanishes and the equation expresses 
the proportionality of V-s and F (gy) (ef. also 2.3) 

V-s = kF (oy). 7.4) 
The condition for vanishing of the coefficient of gy is thus synonymous with the condi- 
tion for wholly parabolic flow. Equation (7.4) expresses a simple geometrical property 
of parabolic flow, namely that the fractional rate of change of stream tube area in the 
direction of flow is a function of gy only. Accordingly, the fractional rate of change of 
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area of every stream tube intersecting a given potential surface is the same at that 
surface. 
Because (3.11) and (7.4) do not contain gy we use (3.10) to express its variation 





or = 1N& - NY (v EL y? — ony + i) Flew) 
2 Y 
(7.5) 
] ses i 
+5 (k — N) (k ~ itt N) Ven 
ef 
The “Bernoulli” equation becomes 
log j a + | Fes dgy = constant, (7 6) 
and the relation between W and N 
— y¥+1 ‘) me 
W* = . L w(k a N (7.7) 
corresponding to W* = 0, (y — 1)/(y + 1), 1 for (N/k) = 0, y/(y + 1), 1 
8. Hyperbolic flow with D, = —(y + 1). The character of g which is sufficient if 
the coefficient of F(gy) in (3.11) is to vanish is expressed as 
g = y/N’. (8.1) 
This corresponds to a type of wholly hyperbolic flow specified by D) = —(y + 1), 
k—o and toG = +(y + 2)/N’, H, = N’. The equation for ¢y is 
m ie dex | dyn at don Ion on 
> | 1+ 2)( Ses 7 Mt > ax; ax; at, 2, 
(8.2) 
1 
= 1£* N*F(oy) 
¥ 
or in the condensed form of (3.9) 
ad 5 aoe 
V-(N-?N) = 1s "Flgy). (8.3) 
The Glauert-Prandtl approximation is 
d°e 3° 3° y¥+1,, 
-@ + 0+ += NoF (en). (8.4) 
v3 ¥ 


Since the coefficient of F(¢gy) in (3.11) is zero, gy is given by the simple expression 
has . 
Iw=5 NV-s. (8.5) 


In this case the area change of stream tubes is determined solely by gy and N. The 
Bernoulli equation is 


log (pN’) + | Fes) dygy = constant, (8.6) 
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and the relation between W and N 


, = 3 
W* = ——., (8.7) 
YT I 
This particular wholly hyperbolic flow is thus characterized by a constant value of W 
(or M = 1) throughout the field of flow. 


9. Mixed flow. It has been seen how proper choice of the character of the flow, as 
defined by the function g(N), can establish flow whose type is constant throughout the 
field of flow. That the type may change within the flow field is illustrated by reconsidera- 
tion of the W and M fields of irrotational diabatic flow (Reference 1). (The functions gw 
and gy are expressible in the form of (5.9) and accordingly admit the incompressible as 
well as the Glauert-Prandtl approximation; see Reference 1.) The partial differential 
equation for yy is of elliptic, parabolic or hyperbolic type accordingly as M(=N) <, =, 
or > (1/y)'”’. The partial differential equation for gy changes type similarly at W = 
N = [(y — 1)/w+ 1)”, (M = 1). It is the latter field which completes the set of three 
flows determined by the vanishing in turn of the coefficients of gy , F (gy) and V-s in 
equation (3.11), for an irrotational W field (Reference 1) possesses the property that 
— | 


Pm 
| Maia 


— F(py) (9.1) 
2y 


10. Concluding remarks. If, in the definition N = V[g(N)RT]'’’, g is taken to bea 
function of N alone, any N language yields a simpler set of fundamental differential 
equations for rotational diabatic flow than the V language. In particular the W language 
(derived from g = 2y/(1 — N*)(y — 1)) corresponds to the most simple expression of the 
equation of motion and hence is best suited to describe a general rotational diabatic flow. 

With proper choice of the function g(N), which is taken to determine the character of 
the flow, the partial differential equation for the potential fuaction gy of irrotational 
N flow can be made wholly of one type (e.g. elliptic) or of mixed type throughout the 
field of flow. The differential equation for gy always contains an arbitrary function F (gy) 
which fixes the form of the flow and also occurs in the analogue of the Bernoulli equation. 
Choice of the flow character alone leads to the form of the relation among F(gy), V-s 
(the fractional rate of change of stream tube area in the direction of flow) and qy (pro- 
portional to the local rate of heat addition per unit mass of fluid). Because of the universal 
occurrence of F (yy) in the differential equations for gy there is need for further study of 
such quasi-linear partial differential equations containing arbitrary functions of the 
dependent variable if irrotational diabatic flows are to be characterized more fully with 


any degree of generality. 
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ON THE PLASTIC BENDING OF CIRCULAR PLATES UNDER UNIFORM 
TRANSVERSE LOADS’ 


BY 
D. TRIFAN 


Brown University 


1. Introduction. The primary purpose of this paper is to investigate the mechanical 
thin plates subjected to distributed transverse loads which produce plastic 


behaviot r ¢ 
(hroughout this discussion the plates will be considered as being homo- 


deformatio1 
composed of isotropic materials exhibiting a gradual transition from the 


plastic state. Two alternative theories of plasticity, namely a theory of 
plastic flow and a theory of plastic deformation,’ will be used, and their predictions 
will be compared with each other. In particular, the case of a circular plate subjected 


geneous and 


elastic to the 


to a uniformly distributed transverse load will be considered in detail. 

\s applied to strain-hardening materials, the theory of plastic flow is based on the 
assumption that the rate of change of strain is uniquely determined by the existing 
he rate of change of stress. Stress-strain relations of this type have been 
considered by Prager [2, 3] and, more recently, by Handelman, Lin, and Prager [4]; 
they may be considered as a natural generalization of the stress-strain relations proposed 
for perfectly plastic materials by Prandtl [5] and Reuss [6]. On the other hand the theory 
of plastic deformation developed by Hencky [7] for perfectly plastic materials and 
applied by Nadai [8, p. 75] to materials with continuous transition from the elastic to 
the plastic state is based on the assumption that in the plastic state as in the elastic 
tate of strain is uniquely determined by the state of stress. These stress- 


strain relations [9] are simpler than those of the flow theory, and thus have received 


consideration in mathematical studies concerning permanent deformations in 


] 
Stress and 


state, the 


more 
structural elements. 

Until recently, comparatively little attention had been given to the plastic bending 
of plates; however, some notable contributions in this field have recently been made by 
Russian investigators, in particular, Ilyushin [10] and Sokolovsky [11]. Most of their 
results, however, do not apply to materials which exhibit gradual transition from the 
elastic to the plastic state. 

2. Stress-strain relations. Investigation of the bending of plates made of a strain- 
hardening material by means of the theory of plastic flow is extremely difficult if the 
customary form of the stress-strain relations for such materials is adopted [4, Eq. 27]. 
Indeed, from this form of the stress-strain relations, the rate of change of the stress 
tensor can not be obtained explicitly in terms of the rate of change of the strain tensor, 
the strain tensor itself, and the strain invariants, as is necessary in the study of the 
bending of plates. This difficulty is avoided if slightly modified stress-strain relations 
for loading are adopted. Let the tensors of stress and strain be denoted by p,; and «,; , 
their deviations by 


*Received Jan. 30, 1948. The conclusions presented in this paper were obtained in the course of 
research conducted under Contract N7onr-358 sponsored jointly by the Office of Naval Research and the 


Bureau of Ships. 
1This terminology is due to A. A. Ilyushin [1] (numbers in square brackets refer to the bibliography 


at the end of the paper). 
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1 
rT; = a 405; e; == 2S 6;; 
i Pij 3 PROG i i 3 Reo Gj 
(6,; = Kronecker delta), and the elastic shear modulus and the tangent shear modulus 
by G, and G, respectively, the latter being assumed as a function of the strain intensity 
2 
E = 3 C530 i;5 


The stress-strain relation of the flow theory used in this paper then has the form 


E* 
rt, = 2G e%; om (Go —- G@) E Ciz (2.1) 

Throughout the following discussion, the Latin subscripts run from 1 through 3, and 
rates of change of variables with respect to time will be indicated by asterisks. Stress- 
strain relations of the type (2.1) may be termed stress theories of plastic flow, while 
relations such as are treated in [4] may be called strain theories of plastic flow. 

The expression (2.1) constitutes five linearly independent stress-strain relations 
which take into account only the change of shape of a deformed body. To these, a sixth 
expression dealing with the change in volume must be added. For the purpose of simpli- 
fication, we shall assume the material to be incompressible: 

€é,, = 0. (2.2) 
(Here, repeated indices indicate summation in accordance with the customary summa- 
tion convention of tensor calculus.) Accordingly, Poisson’s ratio vy = 1/2, €,, = @,,; , 
and Young’s modulus equals three times the shear modulus G, . 

In order to determine the quantity G which is a function of the strain intensity F 
only and thus independent of the specific type of loading, it is sufficient to consider a 
simple state of stress, e.g., simple tension. Indicating the tensile stress by p and the 
corresponding unit extension by e, Eq. (2.1) reduces to the single equation 


d ‘ ae 
—P — 3G, (2.3) 
de 

and the strain intensity becomes E = e¢’. Hence the tangent shear modulus G for any 


material obeying (2.1) can be determined if the stress-strain diagram for simple tension 
is known. For the material being considered here, the function p = p(e) can be repre- 
sented with sufficient accuracy by the power series, 


p=aqetaet+ase+-:--, (2.4) 


where a, , a3 , 4; , etc. are constants. Thus the tangent shear modulus is given by 


3G = a, + 3a,F + 5a;k’ + ---, (2.5) 


where a, equals Young’s modulus, i.e., a, = 3G). 

It is more convenient from this point on to work with non-dimensional quantities. 
Accordingly, the following non-dimensional components of stress and stress deviation 
will be introduced: o;; = pi;/Go , 8i; = 1i;/Go . By substituting these quantities and 
the series expansion (2.5) into Eq. (2.1), we obtain a stress-strain relation for the stress 
theory of plastic flow for a material with a stress-strain curve in simple tension which 
is represented by the series (2.4), or by the equivalent general stress-strain relation 
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st, = 2e*, —_ [b, — b.E a bE” ot atthe -|E*e,; ’ (2.6) 


where 


b, = —(2n + 1)aa,.,/3G , n= 1,2,3,--- 
The stress-strain relation of the theory of plastic deformation is of the form 
r;; = 2Ge,; , (2.7) 


where G’ = G‘(E) is the secant shear modulus. For simple tension Eq. (2.7) reduces to 
the single expression 


p/e = 3G. 
Using the same series expansion (2.4) as in the plastic flow theory, we obtain 
3G. = a, + a,E + a,E? + --- (2.8) 
Substituting this into Eq. (2.7) and using the non-dimensional form, we find 
8;; = 2e,, — [c, + oF + ¢,k* + ---]Ee;; , (2.9) 
WW here 
Co = —Qlens:/3Go - 


3. Basic differential equations. We consider a thin plate subjected to distributed 
transverse loads, and assume that the resulting deflections are small in comparison to 
the thickness of the plate. We introduce the rectangular coordinates X; , the X, , X, 
plane coinciding with the plane horizontal middle surface of the undeformed plate, and 
the axis of Y, being directed vertically downward. As is usual in the theory of thin 
plates, the stress components o;3 , ¢23 , ¢33 Will be treated as small in comparison with 
the stresses parallel to the middle surface of the plate (¢,, , o22 and o,2). Restricting 
Greek subscripts to values 1 and 2, we may, therefore, write Eq. (2.6) in the form 


210%. + Saspl aot (an T Celt ys }[d, a b.K — -JE*, (3.1) 
where / 2/3(« alas + ¢ va€BB)e 
We now adopt the following notations: 
etl oh a a, = 2X./h, 
w = Wh/2R?, Mas = 4Mas/h’G, , (3.2) 


Ll = 4gRi/h'G, , 


where FR, is a fixed suitably chosen length, h the constant thickness of the plate, W = 
i’ (X.) the deflection, g the intensity of the load, and 2.» are the flexural and twisting 
moments (3? J i ey: 

It is a fundamental hypothesis in the theory of elastic plates that points on a normal 
to the middle surface of the plate before deformation remain on a normal to the deformed 
middle surface. It is assumed that this hypothesis remains valid for the plastic materials 
considered here. Accordingly, the relation between the strains ¢,s and the deflection 
W of the middle surface is [12, p. 323] 

€ap = —X3W as; 


where W ., = 0°W/dX,0X;5. 
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In non-dimensional form this expression becomes 
€ap = —X3W, ap - (3.3) 


In terms of the deflection w, Eq. (3.1) then becomes 
2 3 + 5 ’ 
aX, = —27,0% 4 op 3 b, 23, L4Q.8 oo 9 bot3 Bo EF Qa8 fe he dibel” (3.4) 


. , 7 , 2 ° ° 
where the quantity EZ, , related to E by means of FE = 2/3 23, , is written as 


E, = QapW, ap 9 
(3.5) 
Qap = Wap + SapW, 77 - 
The time rate of change of the non-dimensional bending moments 
1 
Mt, = | ots dr, 
-1 
can now be easily calculated. We obtain 
M*, = 4 O* 1 ( - bh. F*ER 26 
a 3 ap ~~ 5 bE ap — 21 Dole 0 Ho Qas ss ae (3.6) 


For transversely loaded plates with small deflections, the equation of equilibrium? 


[13, p. 87] becomes 


which, together with Eq. (3.6), furnishes 
] = 
Oty. ap — = b(EE Gap} 0p — Ze ba{E$B.Mep} ep — °° = SU. (3.8) 


Finally, referring to definitions (3.5) and carrying out the necessary steps, we arrive 
at the basic non-linear partial differential equation in w for the flow theory, namely 


b, 
) ) ) ) ) ) y* 
w*.0 — 5 {Wap W,ys W*apys + QW, aa Ws Wepys + Waa Wp Wy 758 


, . * On ‘ a 
+ QW as Ways Ways + 2W aa Wpys Ways + 2w ap W ayy Wigs 
+ 6,00 W,pyy W555 + 40,08 Wy 78 wW* 65 + Wap W,apys Ws 
(3.9) 


, , ’ * 
+ Waa W,ppy5 Wes + Wap W.apyy Ws + 3Waa W.ppyy Ws: 


+ 4W, app Ways ws + 4Wass, Wary w*s; + 2w, as W +788 w*.s} 


b. en as 3 
— 91! YE Qap}.ap — °° = 8 I*. 


The corresponding basic partial differential equation for the plastic deformation 
theory can be obtained in a similar manner. We note merely the final equation, 


*See Eq. (g): M: = M, ’ Man = M, ’ Me: = —Miy ’ Ma = Myz ° 
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, Says ; ) 
W aapsp ~ iW. as W +s u ,aBys + 2W aa W by W py ss + 20 aa W 6g Wy 78s 


+ Wap Wap W788 HF Was W,ays Wipys + Wap Wavy W988 


(3.10) 
+ 4W. ap W. aby W788 + Waa W678 W873 + SWaa Wipyy W,38} 
ere eee 3 
a 21 {FFE Qos} as Os g 


4. Example. We are now in a position to consider a specific problem in detail, namely 
the case where 

a) the thin plate is made of 24 S-T aluminum alloy, 

b) the contour of the plate is circular, 

c) the plate is built-in along its entire boundary, and 

d) the lateral load is uniformly distributed. 

In Fig. 1 the graph (full curve) of the stress p and the corresponding unit extension 
« is given for a tension test with a specimen of 24 S-T aluminum. The analytic approxi- 
mation for this curve most suitable for our purposes is 


wv.) 
= 36,(¢ — 1% #), , 
ath (. 128 £ 42) 
where 
3G, = 10.667 X 10° Ibs. /sq. in. 


This relation is represented by the dashed curve in Fig. 1. Referring to the expansion 
2.4), we find that the values of the coefficients for this material become 


\é. 


3-10° 
a, = 3G, . Ge ea Gy , a, = Q@,=>¢**' = 0. 


128 


From this it follows (see Eqs. (2.6) and (2.9)) that 


3-10° 
b=, b= b= = 0 (4.2) 
and 
2-10° 
= oer te 


With these values, the basic equations (3.9) and (3.10) are greatly simplified. 

Since the load intensity is assumed to be constant, i.e., 1 = 1, , the deflection w is 
independent of the polar angle 6. The transformation of Eq. (3.9) with the values (4.2) 
into the more suitable polar coordinate system gives a differential equation in w, now 
a function of the single non-dimensional radial coordinate r = R/R, , where R is the 
dimensional radial coordinate and R, the outer radius. The actual transformation is quite 
laborious but straightforward, and for that reason it will be omitted here. Expressing 
the differentiation with respect to r by primes, we find 


wt Pane — Legere g eae — Sap, a9 
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where F = F(r, t, w’, w”’, w’”, w’, w*’, u 
L(t). 

The complete representation of Eq. (4.3) shows that this expression is homogeneous 
in the time rates of change, and therefore the time variable ¢ can be changed to any 


ys 


, wr’, w**) represents 24 terms and |, = 




















Fig. 1. Stress-strain diagrams: (1) 24S-T Aluminum Alloy (full curve); (2) Approximation as given by 


Eq. (4.1) (dashed curve). 
other monotonically increasing variable. Since we assume the load to increase with time, 
t can be replaced by I, , thus, w* = I*(dw/dl,). 
It is found that if Eq. (4.3) is multiplied by the factor r then the differential equation 


can be integrated to a lower order expression, namely 


1 b l 
rwr!! + we! ” 4 w*’ me 5 Arw’’w'’ + 4w'w"’ +- 4 w'w' wt!’ 


oh (sru’’w’” o Rw’ w’’ 4w’w’"’ ot 2 wo’ ‘ai 4 w'w! Jw” (4.4) 


Matt? 4 tents 5. we 9 ee | w!)ur) = 39 
+ (40 + ww + oe ww 7p ww jw = 167 
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which serves as the basic differential equation of a circular plate under uniform load 
for the plastic flow theory. 

The corresponding basic differential equation for the plastic deformation theory of a 
circular plate under uniform load can be derived in a similar manner; one obtains 


ies l c 3 
wT, 4 w’ —_ w’ = <I 6rw’ ww!’ + 6w’ w/w!’ + ~ w'w'w'"’ 
r oO | r 


(4.5) 
3 
16 


‘7 , , 3 , , 3 4 
+ 4w’’w'’w"’ + -— www" — sw’w'w"” — = ww’? = lor 
r r 
The problem of solving the non-linear differential equations (4.4) and (4.5) for 
w(r, lb) presents great mathematical difficulties, and there is little hope for obtaining 
solutions in closed form. We therefore will attempt to approximate the solutions by 


means of series of the following type 
w(r, lo) = w,(rylo + we(r)li + ws(r)lo + +: , (4.6) 


where the coefficients w, , w; , w; , ete. of the odd powers of the parameter |, depend 
on r alone. Substituting this expansion into Eq. (4.4), and collecting terms involving 


like powers in /, , we obtain a sequence of differential equations in the unknown functions 
W., - . The first of these is 
3 
yyy rv , 7 2 
rw} t-wol—--w=—r 
“a 16 
(4.7) 


| d E d (rd) . By 
™ "adrtrar\’ adr ~ 16"? 


the next contains w, and w, ; only w; is unknown, however, since w, can be determined 
from Eq. (4.7). This procedure can be continued until the desired number of coefficients 
in (4.6) are found. 

In view of the boundary conditions for a built-in circular plate, 


w(i, ¢,) = 0, w’(1, 1) = 0, (4.8) 
and the requirement that w(0, /,) be finite, the solution of Eq. (4.7) is found to be 


W, a 1)’; (4.9) 


= 512° 


this is nothing but the elastic solution. The differential equation for w3(r) then becomes 


_adj1d dus) _ 2b, , S _ 1904 2 
"ar E dr (; dr 5(128)" {1437 129r° + 30r°}, (4.10) 
and that for w; 


a |? d ( dws) = 83-10 1350,493r'° — 622,668r° 
drrdr\ dr 25(128)’ 





(4.11) 
+ 419,112r° — 128,067r* + 15,270r*}. 


Integrating these equations and substituting the results into the series (4.6), we obtain 
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45 


{143r° — 344r° + 360r* 
4(12.8)* 


3 2 ’ 
wr, lo) = 512 (7 1)"ly + 


) 3 324 o~ 2 10 
— 260r° + 101}l, + ; - {584,155r'* — 1,868,004r . 
(12.8) i (4.12) 


+ 2,619,450r° — 2,134,450r" + 1,145,250r* 


— 529,860r° + 183,459}/3 +--+. 


This equation is the series development up to terms of the 5th order in /, of the solution 
w(r, ly) to the differential equation (4.4) of the plastic flow theory. 

Beginning with the differential Eq. (4.5), a similar expansion can be obtained for 
the deformation theory. The functions w,(r) and w,(r) in the series solution as given 
by the deformation theory are found to be identical with the corresponding functions 
obtained for the flow theory. This property can be shown to exist even under more 
general conditions; however, this fact will not be proved here. The first difference be- 
tween the two theories occurs for the function w, . Up to terms of this order, the series 
representation of the deflection furnished by the deformation theory is found to be: 

ir ° 45 


Gg? = 17h, + - {143r° — 344r° + 360r° ‘ 


wr, &b) = > o 
512 4(12.8) 


») 5 ( 2 0 
— 260r° + 101}/) + a - {352,495r'* — 1,130,796r' . 
(12.8) (4.13) 


+ 1,591,650r° — 1,300,750r° + 697,500r* 


— 320,340r° + 110,241}1, + ---. 


The stress distribution throughout the circular plate can be calculated in the case 
of the plastic flow theory by transforming expression (3.4) into polar coordinates, 
substituting Eq. (4.12), and then integrating. The resulting non-dimensional expressions 
for the radial stress component P and the circumferential stress component Q together 
with those similarly obtained for the plastic deformation theory are valid for all values 
of 1, < lj , where lj represents the maximum load which our assumed stress-strain law 
entitles us to consider. In order to calculate this load we return to Eq. (4.1). For uniaxial 
stress /{ would correspond to the maximum height of the dashed curve in Fig. 1; for the 
combined state of stress we are considering, ihe stress-strain curve expresses the relation 
between the intensities of stress and strain, where the stress intensity 

T = 3/2 7; ;73; ; 
reduces for simple tension to 7’ = p*. We have already seen that E = 2/3 e;,e;; reduces 
to ¢ for simple tension. Thus, for combined states of stress, the coordinate axes of the 
dashed curve of Fig. 1 are labeled T’”’ and E’”’. 
The non-dimensional stress intensity S = 7/G> = 3/2 s,;8,;; for the example being 


considered here is 


S=P?+ Q@- PQ. 
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Since the maximum value of S'”’ occurs at the boundary, i.e., zx; = +1, r = 1, we find 
that for l, = lf = .09 a value S'’” = .01312 is reached which takes us to the top of the 


stress-strain intensity curve in Fig. 1. 
A review of the stress distribution on the upper surface of the plate shows that 
for 1, = .02, see Fig. 2, there exists no appreciable difference between the data as given 


016 


Ol2 


-,008 


-.0l2 





Fic. 2. Curves for the non-dimensional radial and circumferential stresses in the upper surface of the 


circular plate for lo = .02 and ly = .09 as given by the plastic theories (full curves) and the elastic theory 


(dashed curves). 


by the elastic theory and the two plastic theories; to within five significant decimal 
places, the third order terms in |, have a negligible influence, and those of the fifth 
order have no influence at all. As the load increases to 1, = .04, there is a slight deviation 
which becomes more and more pronounced as I, increases; however, to within five decimal 
places there is at no time a difference between the values given by the theories of plastic 
flow and of plastic deformation. For 1, = .06 the third order terms have an influence 
ranging to 3.46% of the first order terms, while the fifth order terms reach only .53% 
of the first order terms. Third order terms range to an influence of 7.8%, and fifth 
order to 2.5% for the ultimate load |, = lf = .09. From this it is estimated that the 
solutions given by (4.12) and (4.13) approximate the real solutions with sufficient 


accuracy. 
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In order to obtain the work D done on the plate we calculate the following integral: 


Ro t a 
D = 2 [ [ Jo a dtRdR 


i _ was ae 990 
2r ' qgW — uv Ot dt |R dR. 


0 


In non-dimensional form—see (3.2)—the above reduces to 


; . 1 ot dl, 
D = rG,hRo J, E = / w 7 at dr 


or by changing the time parameter again to /, we obtain 


3,00 x10° 


200xI0° 


02 04 06 





08 


Fia. 3. Relation between non-dimensional loads I, versus non-dimensional work d as given by Eq. (4.14) 
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1 


D=35 


GohR d, 


..here the non-dimensional work d is given by 


1 le 
d = 2n / | tw _ [ w at} dr. (4.14) 
. ; 


0 


A graph of d/2z7 versus I, is given in Fig. 3. We note that in the purely elastic case, 
where w reduces to simply the first term of (4.12) or (4.13), the above formula (4.14) 
reduces to the familiar form 


d= r | l, wr dr. 


In conclusion we can state that though the theories of plastic flow and plastic 
deformation are founded on different hypotheses their predictions for circular plates 
under uniform loads and for the specific material considered here are identical for all 
practical purposes. Both plastic theories however, give a deviation from the elastic 
solution ranging to better than 10% depending on the load and the material point 
under consideration. 
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A MATHEMATICAL FORMULATION OF THE GENERAL CONTINUOUS 
DEFORMATION PROBLEM* 


BY 
A. GLEYZAL 
David Taylor Model Basin, Navy Department, Washington, D. C. 


1. Introduction. The purpose of this paper is to formulate precisely and simply the 
general stress-strain problem for continuous flow or continuous deformation of con- 
tinuous media such as viscous liquids or elastic or plastic solids. We shall use an embedded 
coordinate system which is considered to be embedded in the material so that each 
particle of the body carries with itself coordinates x = (x’, x’, 2*) which’ vary contin- 
uously from particle to particle but which do not vary with time for a given particle 
as the body deforms. 

The square of the are length for an embedded element dx’ = (dx’, dx’, dx*) is given 
by a quadratic form ds* = g;; dx dx’, summed’ on i = 1, 2, 3 and j = 1, 2, 3, where 
the fundamental tensor g;; is a continuous function not only of zx‘ but also of the time ¢. 
This changing fundamental tensor is of basic importance in the theory and will be called 
the defon mation tensor. 

In terms of the deformation tensor the following two principles may be stated: 

|. Every continuous flow or continuous deformation may be described by the deformation 
tensor g,;(x', t) where g,, ts the fundamental tensor of space referred to embedded coordinates 
in the body. 

2. All equations relating stress and finite strain must be expressible as tensor equations 
relating a stress tensor o;; and the deformation tensor g;; where o;; and g;; are referred to 
the same embedded coordinates in the body. 

Deformations or stress-strain laws which cannot be so expressed have an indefinite 
meaning since the deformation or stress-strain properties of the material would then 
depend upon the choice of coordinate system. The above principles immediately impart 
conceptual ease and clearness to the general precise stress-strain problem for continuous 
deformations. In many cases the writing of the complete equations for such a problem 
becomes mechanical. 

Elementary examples to illustrate the procedures and equations are given in this 
paper. A tensor generalization of ‘‘natural strain” is found and Young’s modulus and 
Poisson’s ratio for “‘finite strain’? are derived. Infinitesimal strain and strain velocity 
are defined for continuous deformations and related to the deformation tensor. 

The general problem of stress and finite strain for elastic bodies has been considered 
by Murnaghan (3), Biot (4), and Seth (5). The author believes it would be of interest 
to translate the stress-strain relations assumed in these papers into tensor equations in 


*Received March 6, 1948. 

‘Unless otherwise specified all superscripts or subscripts such as 7, 7, k, a, 8 take on the values 1, 2, 
and 3 

2An expression containing a letter which is repeated as a superscript and subscript is summed on 
that letter unless otherwise stated. 

’The reader is referred to Eisenhart [1] and McConnell [2], for accounts on tensors and their applica- 
tion in geometrical and physical problems. 

‘Numbers in brackets refer to the bibliography at the end of this report. 
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embedded coordinate systems. Definitions of finite strain based on geometrical quantities 
may also be so translated. 

This use of a moving and deforming coordinate system may be said to be an in- 
tensification of the Lagrangian viewpoint, where the position and velocity of a particle 
referred to space coordinates are expressed as functions of embedded coordinates and 
the time ¢. It is closely related to the intrinsic viewpoint in geometry, where quantities 
such as the curvature of a surface are expressed in terms of the fundamental tensor of 
the surface itself without reference to the space in which the surface is immersed. 

2. The general problem. Consider a body or continuous medium immersed in space. 
Let x", a = 1, 2, 3, be the coordinates of points of space referred to a fixed coordinate 
system. Let x’ be a set of coordinates embedded in the body so that x’ are the coordinates 
of a particle of the body. Any deformation may be uniquely described by three dis- 
placement functions: 

f = 2c, }) (1) 


which preseribe the position %* of the particle x’ at time ft. 
The square of the element of arc length in the space is given by a quadratic expression: 
ds’ = Gap ax dx’, 
where the coefficients g,3 are known functions of x*. Often Cartesian coordinates may 
be chosen so that 9,;; = goo = Js3 =1 and gi. = Qos Js, = 0. In the embedded co- 


ordinates 


lo " 3 
' O2 ,\f 0x 
ds = g;; dz’ dz’ = Gos dz ) — as ) 
Ox !\@2’ 
Therefore g;; , a function of x‘ and ¢, is given by 
Ox 02 (a 
Jii = Jas ; , “) 
Ox Ox 


so that g;; is expressible in terms of g.3; and Z“(z’, ¢). 
We shall make use of transformation equations such as 


where H, , Aas are the components of a tensor in Z* coordinates and yu, , A;; are com 
ponents of the same tensor in x' coordinates. The formulas for raising or lowering indices 
of tensors yu; or \;; are formulas such as 


w= 9 Ur, AH =GO Mi, Nez = Gir VN, 
where g'” is the “reciprocal” of the deformation tensor g,; , that is, 
9°'Gr5 = 5; , 
where 6; is the Kroneker delta: 
6; = | if i =4J, 5; = 0 if 1 j. 


It may be noted that the six components g*’ may each be written in terms of the six 





tet 
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components g;; by the rules for finding the reciprocal of a matrix. It follows that g"’ 
may also be expressed in terms of 9,3, and x“(z", t). Thus we may write 
g - 9 (Gas ’ z*). (3) 
In accordance with the principle stated in the introduction there exists for any 
material a stress-strain law of form 
03; = gilgis(2', D, x, d, (4) 


a functional’ of the six functions g,;(z', t) and may vary also with the 
numbers x’ and ¢. In many cases it suffices to choose as g,; a tensor formed from g;; 
and g"’ by the processes of addition, multiplication, differentiation and integration. For 


where ¢;; is 


example ¢;; may be formed as the sum of such tensors as” 


at at 


at t 
Gii | Jii dt, Jii | Y ij dt, | [ Jii dt, Jii » 9 'GriQei 


“oO “0 “0 “0 


ete., with coefficients which are functions or functionals of invariants such as 


t 


9° Gis 9" [ gi; At, 9g" Grud: 


etc. In short ¢g;,; is any covariant tensor of order two which may be formed from g;; 
and g’’ by tensor addition or tensor multiplication and by differentiation or integration. 
For a given material ¢;; is determined by experiment. 


The equilibrium conditions may be written: 


VW here 


{ i 6) 


and 


J t a , 1 (99:, O9:: 201) = 
° = ¢ [ k, Tis k, ri=- (<2 a a. «<a (7) 
J , OM Nae Pas a 


Equations (1) to (7) combine in a straightforward manner to form three equations in 
the three unknown functions * = %*(z', t) which may readily be written out. Solutions 
of these which satisfy the boundary conditions of the problems are sought. These 
boundary conditions may be, for example, specified stresses ¢*° referred to space co- 
ordinates. They may become boundary conditions on the three functions by use of the 
equations: o;; = i; = Fas OX" /dx' dx" /dx’. 

3. Relationship of G;; and infinitesimal strain. To translate Eqs. (2) into more 
familiar terms let us transform both the embedded and the space coordinates, in a given 


5A functional is an association of a number with one or more functions. 

6A dot placed over a quantity will indicate its time derivative when the x‘ are kept constant. 

7It may be readily verified that the derivative or integral of a tensor with respect to time is also a 
tensor with regard to any transformation of coordinates from one embedded system to another embedded 


system. 
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deformation, in such a way that both coordinate systems become Cartesian and have 
identical axes at the instant ¢ = ¢, . Then we may write 





=2* +1", 
where u' = u'(2’, t) are “displacements” and u'(x’, ¢,) = 0. The one-one component 
of the deformation tensor is, by Eq. (2), 
du' \° du’ \” du’ \* 
gi = \1l +5 + ' T — (4 
\ Ox 02 Ox 
and the one-two component is . 
* 1 n 1 ” 2 “ 2 * 3 4 3 
du\ou , Ow Ou ou ou 
ca = ( + = ) ae. eae (1 + :) + om fem : 
Ox / Ox Ox Ox Ox Ox | 
with similar expressions for the remaining components. It will be noted that these 


expressions are, except for a factor of two and the appearance of unity in some of the 
terms, identical with well known ‘‘exact”’ expressions for the “finite strain tensor’. 
If ¢ is infinitesimally near to ¢, then w’ is infinitesimal and 
Qi =1+ so. fiz = ~ ~ 
Ox Ox Ox 
with similar expressions for the remaining components. Now du'/dx', 1/2(du'/dx* 4 
du’ /dx') ete., are components of infinitesimal strain as usually defined* for the situation 


which exists at time ¢ = ¢, . The covariant components of strain velocity v;; may be 
defined to be the time derivatives of these infinitesimal strains at time ¢f = ¢, , so that 
Is Ou’ | - l ( Ou ~) 
qi =o? oF ae 44) Me * Uy = a 3 as + <- 7 ’ ete. 
Q 1 Ox" at 2% ‘2 2\ex’ at © dz’ at 


Alternatively, we may write 


1: Ov’ = | & Ov 
: — }, = . _S-YV.0=> —= — 5 »f - 
Q 9 ” dx’ Q 9 2 \da" si dx’ _ 
where v' = du'/dt (x; = const) is the velocity vector (in either coordinate system) at 
time ¢ = ¢, for the particle x’. In tensor form, therefore, the equations for strain velocity 
are 
l : l 
59 ayy. 5 T 0; i), 


; is the covariant derivative of v; defined by the equations 


av; _ fr 


Ui, aad — v;) . fi 
Ox’ vj 


{.";} being the Christoffel symbols as defined by Eqs. (7). 
Since g;; , v;; and de;; are each tensors with regard to any change of embedded co- 
above results do not depend upon the choice of embedded coordinates 


where v 


ordinates the 
and are valid at all times ¢. 


8Except, by some writers, for the factor 1/2. 
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The infinitesimal strain at time ¢ for the time increment di is simply 


de;; = v;,; dt = 5 9 dt = 5 ag, : 


Thus, assuming for convenience that g;; = 0 whent = —o, 


f de;; = 5 9:3 


20 


and we may say that one of the integrals of the infinitesimal strain tensor is one-half 
the deformation tensor in an embedded coordinate system. 

By the principles stated in the introduction it was to be expected that the infinites- 
imal strain tensor de;; should thus be expressible in terms of operations on g;; since 
surely the infinitesimal strain tensor is independent of any rigid motion of the body. 

It is instructive to consider “deformations” for which de;; = 0 at all times ¢. Such 
a condition may be considered to be the equivalent of a stress-strain law in a problem. 
If de,; = 0 then v;; = 1/2 g,;; = 0 and the deformation tensor g,; is a function of zx‘ 
alone and does not vary with time. Hence the arc length ds’ = g;; dx‘ dx’ of anembedded 
element dx’ is a function of x’ and dz‘ alone and does not vary with time since z’ and 
dx’ do not vary with time. Thus the “‘deformation” consists of any rigid motion. The 
solutions of Eqs. (2) for the displacement 7* = Z*(z', t) when x* denote Cartesian co- 
ordinates, are simply the general equations for transforming from one set of Cartesian 
coordinates to another such set. 

Example 1. We solve a problem of a body under uniform hydrostatic pressure for 
the case of finite strains using a set of generalized elasticity laws. We choose as stress- 
strain law referred to embedded coordinates, the law: 


Ci; = Gy; = Ae;; + Bog;; , (8) 
where A and B are “‘elastic constants’’, 
l ‘ re” ] ; rs 
5 = ir | TGs UFZ Gir | Gai Ul (9) 
and 
if . 
§= 5 [ g°'gi; at. (10) 


The tensor ¢;; as thus defined will be shown to be a generalization of ‘‘natural strain’’. 
It is easily verified that @ = V/V, where V is an element of volume of the body. Let 
us choose the space coordinate %* to be Cartesian, and the embedded coordinates x* to 
coincide with z* at t = 0. Therefore 


Gu = G22 = 933 = l, fie = 923 = J: = 0. 


As trial displacement functions we take 
t 


Z =arzr, 


where 


a = a(t), a(0O) = 1. 
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Using Eqs. (2) the deformation tensor is found to be 
911 = Jor = Jax = A, Giz = Jos = Jar = 0. 


The contravariant components of g;; are then found to be 


Equations (9) yield by direct computation 
Cin = Cog = Cos = a’ log a, €i2 = €23 = €3, = O. 


We note that raising one of the indices of ¢e;; by the equation e; = g*’e,; , we get 


1 2 3 1 2 3 
é: = & = ¢, = loga, &=G&= 6 = 0. 


It may be noted that e;; and €; are considered to be different components of the same 

strain tensor. It turns out for this example, that ¢} and ¢} have identical values as the 

components of the strain tensor or stress tensor, respectively, calculated in experimental 

work. The quantity @ is found to be, for this problem: @ = log a’. Since z' = az‘ and 

z are Cartesian coordinates it follows that a* = V/V» where V is the volume of the 

body at time ¢and V, is the volume at time t = 0. Therefore 6 = log (V/V) or 6 = V/V. 
Combining these equations we find 


o; = g; = Ae} + BOS; = (A + 3B)(log a)5; . 


All components g,;; are functions of time alone and do not depend on x’, x’, or x*. The 
Christoffel symbols [ij, k] and {;,} are linearly dependent on dg;;/dx* and consequently 
vanish identically. It follows that ¢}., = 0¢;/dx*. Since the components are functions 
of time alone the quantities dy;/dx° vanish identically and we get ¢;,, = 0. Thus the 
three equilibrium conditions ¢;,; = 0 are surely satisfied. 

The stress tensor at the surface of the body referred to Cartesian coordinates due to 
a hydrostatic pressure p is ¢3 = pés. Hence the mixed components of the stress tensor 


referred to embedded coordinates are 


ax‘ ax” hh i 
oi = pis —— = psi 5) as! = po}. 





Thus, the boundary conditions of the problem are satisfied if we relate the parameter a 


to the pressure p by the equations 


o, = pd; = (A + 3B)(log a) Si; 


Paar e 


or 


p = (A + 3B) log a. 


Therefore, 


l | 
p= 3 (A + 3B) log vy," 


This relation for pressure and volume of a body is algebraically simple because of 
the particular choice of definition of finite strain (9) and the stress-strain law (8). 
Equations of more complexity are obtained with other stress-strain laws. For example 
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we might have chosen e,;; = 1/2 [g;;(x‘, t) — gi;(x*, 0)] instead of Eqs. (9). In this case 
ls 

= i = 2 (a = 1), €i2 = &3 = és = 0 


and 
The invariant @ is 


The equilibrium conditions ¢;,; = 0 are again satisfied identically. At the surface of 
the material we have 


o} = Ae} + BOS; = pd; . 


Therefore, 


or, since a’ = V/V,, 


a 1 . “3 ‘Ad 
p=5(A + 3B)| 1 (¥ 


Example 2. We now apply the method to the tensile test where the stress is given 
by « = a(t), (0) = 0. This stress is applied, let us say, in a direction parallel to the 
z'-axis of a set of Cartesian coordinates Z', Z’, z*. Let us write for the displacement: 


1 1 2 2 =3 3 
<= az, x = br’, xz = be’, 


where xz’, x’, x* are embedded coordinates and a = a(t), b = b(t) with a(0) = 1, (0) = 1. 


Then, 
ds? = a*(dx')? + b?(dx’)? + b?(dx’*)’ 


and by inspection the deformation tensor is 
Ji. =a, Ges * Gas = b’, Gis @ Gn ™@ Fn = 0. 


The reciprocal matrix g"’ of g;; has components 


a”’ 
As stress-strain laws we choose again Eqs. (8), (9), and (10). The finite strains for this 
case are then: 


> 9 
é,, = a loga, €oo = €33 = 0° log D, 6:2 = 3 = 6, = 0 


or 


1 2 3 1 2 
e, = loga, e@ = «¢ = log b, é& = «6 = €¢, =.0. 
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Also, 
l =. 
7s [ g''gi; at = log ab’. 
“a Jo 
Therefore 


= A loga+ B log ab’, go = ¢, = A log b+ B log ab’, 


§ 


1 2 3 
=o = 4 = 0. 
All deformation components g;; are functions of time alone and do not depend on 
x’, x’, or x. Hence by using the same reasoning as in Example 1 it is found that the 
equilibrium conditions ¢;,; = 0 are satisfied. The stresses on the boundary of the speci- 
men referred to Cartesian space coordinates are 





—1 —2 —3 —1 -2 —3 
1 =a, do = 03 = 02 = 63 = go; = 0. 
U - I - sf . ati ls _ : 
sing the transtormation laws ; 
: i 
j cing MO OS. 
eS 
Ox Ox 
we find 
1 2 1 2 
1 =a, 02 = 03 = 62 = 673 =~ 0, = O 


as the mixed components of the stress tensor on the boundary of the specimen referred 
to the embedded coordinates. Combining these results we find 


o¢ = A loga+ B log ab’, 
0 A log b + B log ab’. 


II 


Solving for log b by means of the second of these equations there results 





B 
log b i+ 2B log a 
Therefore 
_ 4143/4, 
aiid. i Ga 


These equations are of the form of the elasticity laws for infinitesimal strains in a tensile 


test if we substitute for infinitesimal strains, the “finite strains” ¢, = log a, 2 = 6 = 
log b. Comparing coefficients we are led to write two possible definitions of a generalized 


Poisson’s ratio v. Either 
tte (14%) + (1+ 2) 


- B 
A + 2B 


or 


v 


These two definitions, it conveniently turns out, are equivalent. Hence, a generalized 
Young’s modulus is found to be 
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E = A(A + 3B) + (A + 2B). 


It is furthermore seen that the finite strains «| , 6 , €s , as defined by (9) are in the case 
of a tensile test the quantities called natural strains. 

The author has attempted to obtain the exact solution of the problem of the hollow 
cylinder under internal pressure or pure torsion and of a parallelepiped under pure shear 
by the general method given here. The precise conditions on the three displacement 
functions are readily obtained in each case for a number of stress-strain laws but their 
solutions require lengthy computations. The author hopes that they may be carried 
out in future work in this field. 

It is believed that the concepts given here may be developed further to give a general 
theory of continuous deformation in which the notion of energy is fundamental. The 
author hopes to pursue these and other questions in a later paper. 
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—NOTES— 
A PROBLEM OF THE TEMPERATURE DISTRIBUTION IN A MOVING MEDIUM* 


By A. E. BENFIELD (Cruft Laboratory, Harvard University) 


The thermal problem solved in this paper is primarily of geophysical interest. How- 
ever, it seems worth while to give an account of it here, as the boundary conditions are 
rather unusual and the solution, which makes use of the Laplace method, involves a 
transform not given in the usual references. It is thought that the problem may be of 
interest to those generally concerned with questions of thermal conductivity. The 
geophysical implications will be discussed in another paper, to be published elsewhere. f 

We wish to find the temperature 7'(x, ¢), where x is the space coordinate and ¢ is 
time, in a semi-infinite medium, having initially a constant temperature gradient, A, 
when the medium moves to the left with a positive constant speed, v, as shown in 


os 
ae 


TO 
TCX 


AAAS 


ANY 








AWA s 


—AT 4 


AANA 





cath ATS 


SSS SSS Sa. 





Fic. 1. Qualitative curves showing 7(z, t) at times? = 0, ti, fe = 2t; 
where 0 < 4; < tz. T(0,t) = —N. 





*Received June 15, 1948. 
tA. E. Benfield, J. Appl. Phys. (in press). 
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Fig. 1. No material passes the boundary plane x = 0, but instead, on reaching it, the 
material is annihilated without the evolution of heat. Furthermore, the plane x = 0 is 
maintained at a temperature which decreases linearly with time at the rate \. Therefore- 
the boundary conditions are as follows: when ¢t = 0, 7'(2,0) = Az; and atx = 0, T(0,t) = 
—t, where X is a constant. 
The application of the divergence theorem shows that the differential equation to 
be satisfied is 
K aT ae aT _ aT (1) 
Ox Ox ot 
where K is the thermal diffusivity. The medium moves in the negative x direction 


with this sign convention the speed v = —dx/dt > 0. 














—-CAv+a rf 


~CAV+AIT $ 





Fig. 2a. Curves for 71(z, t) = A(z + vt). T1(0, t) = Aot. 


Fic. 2b. Qualitative curves for T2(z, t). T2(0, t) = —(Av + A). 
T(z, t) = Ti(z, t) + Tr(z, t). 
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The problem may be split into two parts; thus 

T(x, t) = T,(2, }) + T,(z, 8 (2) 

where 
T,(a, t) = A(x + vt) (3) 
as shown in Fig. 2a, and T, , which is shown qualitatively in Fig. 2b, remains to be 
found. The boundary conditions for T.(x, t), however, are not hard to find; by considering 
Eq. (2), and the boundary conditions for T and T, , we see that they must be T,(0, t) = 
—(Av + d)t and T,(z, 0) = 0. As 7;(z, t) obviously satisfies Eq. (1), T.(z, t) must do 


the same. 
Calling 7’, the Laplace transform of T, , defined by 





T. = L{T,} = eT (x, t) dt, (4) 

Eq. (1) becomes 
wee dT, = 7 
K = ae ole pT, = 0 (5) 


where p, in Eqs. (4) and (5), is the usual constant, as used by Carslaw and Jaeger.’ 
The solution of Eq. (5), which is a second order differential equation with constant 
coefficients, is well-known and may be written as 


e~*°*1B, exp {a[v’/4K° + p/K]'”} + B, exp{—z[v’/4K* + p/K]'”}] 
where B, and B, are the arbitrary constants, whose values depend on boundary condi- 


tions. However, as we are interested in a solution of 7, which remains bounded, and in 
fact equals zero, as x ©, we shall set B, = 0 and write 





T, = B, e~**’** exp {—a[v?/4K? + p/K]"’}. (6) 
We must now find the value of B, , which is done in the conventional way by considering 
the situation when x = 0, for which value of z, T, = B, . Since T,(0, t) = —(Av + Abt, 
it follows that when x = 0, T, = B, = —(Av + A)L{t}; and so, since L{t} = 1/p’, 
we can now rewrite Eq. (6) as 
aa — a > /4K? 1/2 
fT, = —(Ao + Noo™ SRI-s ? + fy (7) 


In order to find 7, we now need to know L™ [exp {—2[v?/4K? + p/K]'”}/p’]. It is 
known’ that 
-1] €xp [ally + p)/KY"} = oe r\1/2) | = | 
L | ; = } exp {[+2(0/K)”’} erfe XK) + (vi) 
(8) 


1 coal r\1/2 = 1/2 
+ 34 exp {—2(»/K)”’’} erfe Frac (vt) 


1H. S. Carslaw and J. C. Jaeger, Conduction of heat in solids, Clarendon Press, Oxford, 1947, ch. XI. 
2Cf. Carslaw and Jaeger, loc. cit., p. 270, Eqs. (2) and (3). 
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‘ 1/2 ¢ w? . 2 r : ‘ 
where erfe y = 1 — 2 Joe “ du. If we set »y = v'/4K and make use of the theorem® 
that 


y Pa : ? 1 i 
1} | T(t’) ww} = 5 LIT} 


it simply remains for us to integrate the right-hand side of Eq. (8) with respect to time. 


Related problems have been considered by Jaeger,‘ Paterson® and Horenstein.° 
Proceeding now to the integration we find, integrating by parts, that 


at 9 f - 
| erfe f(t) dt = t erfe f + = [ te" df (9) 
“0 wT J © 
where f(t) = at” + bt’, a = 2/2K” and b = »” = 0/2K™*. Writing 
t = (1/b)(t'’*f — a), we find that 
nf 1/2 af 
-f? “a an wn f ] 1/2¢ (-f? e 
| te’ df = ob erfe f + fo t’“fe” df (10) 


The integral [2% ¢'”’ f e-” df in Eq. (10) may be evaluated by integrating by parts and 
further substitutions, but for the case in hand this is unnecessary. This is because Eq. 
(8) contains not only the expression erfe f, but also erfe g where g(t) = at™'”? — bt'”’; 
and as a result of the symmetry of f and g it turns out that the second term of the right- 
hand side of Eq. (10) cancels in the integration of Eq. (8). 

Therefore, combining Eqs. (9) and (10), we find that integrating Eq. (8) leads to 


xp {—al(b° 4 ibs aa | — <7 : 
p-| set mie va p)/K] 1] = ber 4 :) erfe f + 4e e( _ 3 erfeg (11) 
/ } 


from which it follows that 


a oe ) es Zt ts) en eee nats 9 
T(z, t) = {4 + Jt + vt) eric (Zt (a vie ert ¢ Lena |. (12) 


\ 


It can be shown that Eq. (12) satisfies Eq. (1) and the boundary conditions for 7',(z, 0). 
Furthermore, one can demonstrate that if we allow v ©, Eq. (12) leads to the ap- 
propriate expression for the temperature distribution in a stationary medium.’ 

The solution of our problem, therefore, is the sum of Eqs. (3) and (12). This is shown 
graphically in Fig. 1, which is the sum of Figs. 2a and 2b. 

Equation (12) may be derived by other means than the method used here. For 
instance, one may use the substitution® 

tT, #T."™ 

and make use of the fact that 7* satisfies the differential equation Ka°7* dx” = OT*/dt, 

3Cf. Carslaw and Jaeger, loc. cit., p. 242, th. IV. 

‘J. C. Jaeger, Quart. Appl. Math. 4, 100-103 (1946). 

5S. Paterson, Quart. Appl. Math. 4, 305-306 (1946). 

6*W. Horenstein, Quart. Appl. Math. 3, 183-184 (1945). 

7Cf. Carslaw and Jaeger, loc. cit., p. 45, eq. (4). 

8See,-for instance, P. Frank and R. von Mises, Differentialgleichungen der Physik, F. Vieweg & 
Sohn, Braunschweig, 1935, vol. 2, p. 605. 
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thus transforming the problem into one of a stationary, rather than a moving, medium. 
However, on applying the new boundary conditions, for 7*, there results another un- 
tabulated Laplace transform. In this case, the required transform may be found by a 
differentiation, rather than an integration, but the differentiation seems to involve a 
step of doubtful justifiability. 

In conclusion, I wish to thank Dr. P. E. LeCorbeiller for his interest in this problem 


and for several helpful discussions of it. 


ON A CLASS OF SINGULAR INTEGRAL EQUATIONS 
OCCURRING IN PHYSICS* 


By H. P. THIELMAN (Jowa State College) 


1. Introduction. Certain boundary value problems’ in electrodynamics can be formu- 


lated mathematically as Wiener-Hopf* integral equations. The equations are of the 
rorm 
f(z) = | K(| x — y |)g(y) dy, x>0o0 (1) 
“0 
where K( 2 — y |) is a real function of the absolute value of the difference between 


x and y. The Eq. (1) is an integral equation of the “first kind’, K(| x — y]|) and f(z) 
are given functions, and g(y) is unknown. One advantage of formulating a given prob- 
lem, if possible, as a Wiener-Hopf integral equation is that such equations are susceptible 
to the application of the theory of Fourier transforms. The application of the latter 
theory is especially advantageous when the Fourier transforms of the given functions 
are easily found. The object of the present paper is to show that for a certain class of 
Wiener-Hopf integral equations a more direct method yields the explicit solution in a 
more elementary and simpler way. The present method applies, however, only to those 
integral equations of type (1), and to integral equations of the second kind such as 


f(x) = g(x) — A [ K(| x — y |)g(y) dy, xz>0o0 (2) 
“0 
if K(x / |) satisfies the same linear homogeneous differential equations with constant 
coefficients in each of the regions 0 < y < 2, and 0 < zx < y. It should be mentioned 
here that the special case of Eq. (2) with f(z) = 0, and K(|z# — y|) = e~'*-"' was 
treated by Lalesco® by a method similar to the one given here. 
2. The general theory. We consider an integral equation of the form 


f(z) =x [ K(|x-—yl|)gy)dy 2«>0 (1’) 


*Received May 11, 1948. 

J. F. Carlson and A. E. Heins, The reflection of an electromagnetic plane wave by an infinite set of 
plates. I, Q. Appl. Math. 4, 313-329 (1947). 

2N. Wiener, E. Hopf, Uber eine Klasse singuldrer Integralgleichungen, Sitzungsber. Preuss. Ak. 
Wissensch. 696-706 (1931). 

*Trajan Lalesco, Théorie des équations intégrales, A. Hermann et Fils, Paris, 1912, p. 121. 
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where 
a) fo K(| x — y|) dy is assumed to exist, 
b) K(| x — y|}) is continuous for all positive x and y, 
ce) K(| zx — y|) asa function of x satisfies the same linear homogeneous differential 
equation with constant coefficients in each of the regions 0 < y < z, and 0 < 
zs y. 
Since K(| x — y |) is a function of the absolute value of the difference x — y only, 
we can write 
(F(a — y) £O<y<2 
K(|z—- y|) = 4 
[F(—a + y) if 0 


lA 


ry. 


By hypothesis F(x — y), and F(—x + y) as functions of x satisfy the same linear 
homogeneous differential equation with constant coefficients. Let this be the differential 
equation of lowest order which they satisfy. It must then be of the form 


( > a,D* )u = 0, (A) 


i=0 
where 2] = n, the order of the differential equation. This follows from the fact that for 
any m which is a root of the characteristic equation of the linear differential equation 
satisfied by both F(z — y) and F(—zx + y), the negative of m is also a root. Let us 
write (1’) in the form 


ja) =» | F@ — Woy) dy +r | F(—2 + yoy) dy. (1) 
“9 “a 
Taking derivatives of this function and indicating by F’(x — y) the derivative with 
respect to the argument | z — y |, 
-_ , = ce = 
— F(z ~ y) = F'@ — y), — F(-—zxz + y) = —F'(-x + y), 


dx dx 


we obtain 


f(z) = x | F’(x — y)g(y) dy — X | F’'(—x + y)gly) dy, 


/o z 


f'(z) = A F’'(a2 — y)g(y) dy + 2F’(0)g(x) + [ F’'(—x + y)g(y) dy, 


“0 


{'"(z) = dr | F’"'(a — y)g(y) dy + 2F’(0)g’(x) — A / F’'(—a + y)g(y) dy. 


By mathematical induction it is easily shown that 
f"'(z) = | F°? (2 — y)g(y) dy +r / F°?(—a + y)g(y) dy 


+ QF?'-?(O)g(x) + F°*~?(0)g’"(a) + +++ + F’O)g*~” (2)], 


where i = 1,2, --- , L. 





v 


Ae 
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Forming the sum 


i 


> af?" (2) 


1=0 


we obtain 


l l-1 I 
Dafne) = WDD ak OM Og”), (2’) 
i=0 t=O i=t+1 

and the integrals have all cancelled each other because F(x — y) and F(—z + y) satisfy 
Eq. (A). We have thus found that g(x) must satisfy a linear differential equation with 
constant coefficients of order n — 2. Hence if n = 2, g(x) is explicitly expressed in 
terms of f(x) and its second derivative. Equation (2’) gives a necessary condition to be 
a solution of (1’). In order to show that it is actually a soiution of (1’) we substitute 
the general solution of (2’) in (1’) and find the conditions on f(x) in order that g(x) 
may exist, and that the integral in (1’) may exist. 

It is obvious that the theory just discussed applies not only to equations of the type 


|’) but also to the homogeneous equation obtained from (1’) by letting f(x) = g(z), 
the unknown function. If instead of starting with (1’) we start with the integral equa- 
tion of the “second kind” (2), where K(|x — y') is of the type specified above, we 


obtain in an analogous manner the conditions on g(x) that it may be a solution of (2). 
It should be noticed that the kernels of the form 


K(|\z -—y|) = Doe*'* "P| 2 —y)), 


t=1 

where k; > 0 and P,(| x — y|) are polynomials in | x — y |, are of the required type. 

We shall illustrate this general theory by some examples which in themselves are 
of sufficient generality to warrant an independent investigation. For this reason the 
results obtained will be stated as theorems. The derivations of these results will illustrate 
the applicability of the method to problems in applied mathematics in which the reader 
may be interested. 

3. Application of the general theory. By an application of the above theory we 
prove the 

Theorem I. Necessary and sufficient conditions that g(2) be a solution of the integral 


equation 


f(z) = | ea) dy 2>0,k>0 (B) 


a) that kf(0) — f’(0) = 0, 

b) gc) = oe 

c) that the integral in (B) exist, namely that g(x) be of the order e“ where c < k. 

Proof. We note that the kernel satisfies the differential equation (k? — D*)u = 0. 
Writing (B) in the form 


f(z) = | e “= gly) dy + | e*@"* gly) dy (B’) 


“0 
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and differentiating with respect to x twice we obtain 


f'(z) = —\k [ e*=— o(y) dy + rg(x) + Ak | e*"- ay) dy — dg(z), 


f’"(2) = J e~*=- o(y) dy — 2dg(x) + rR’ | e**- a(y) dy. 


Forming k* f(x) — f’’(x) we get k’f(x) — f(x) = 2dkg(z), or 


k*f(x) — f’’(x) 
g(x) = ~ oe ° 


Substituting this value with x replaced by y into (B’) there results 


ja) =o f RW - SW dy 


ml 6 


ta | ct lRsf(y) — sl ay 


2k . 


k hed ee ee ‘ k ‘4 k(v—-z) ¢ 
=5 J e*=-” f(y) dy + a) ° f(y) dy 


“0 


a x [ eo I") ay — ” ] a Fy) dy. 


Integrating the integrals involving f(y) by parts twice by letting first dv = f’’(y) dy, 
and then dv = f’(y) dy we get 


kz 


f(z) = f(x) - = (kf(0) — f’(0)). 


If this last equation is to be an identity kf(0) — f’(0) = 0, as stated in the theorem. 
The function g(x) has to be of the order O(e*’) in order for the integral to exist. 

As a direct consequence of Theorem I we have 

Theorem ITI. All twice differentiable solutions of the homogeneous integral equation 


g(x) =X [ e*'*¥'a(y) dy (C) 
are 
a) g(x) = ¢,(k'” sin {[k(2\ — k)]'?x} + (2. — k)'” cos {(k(2 — k)]'’x}) 
trot 


b) g(x) = c.(1 + kz) if 1 = 


Nie 


c) g(x) = c,(k'” sinh [k(k — 2nr)’”*]a + (k — 2d)'” cosh [k(k — 2d)'”*]zx) 


if 0 <  < k/2, where the c;(i = 1, 2, 3) are arbitrary constants. 





PD TA Ms 


Vy 
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Proof. By setting f(z) in Theorem I equal to g(x) we obtain kg(0) — g’(0) = 0, 
g(x) = k’g(x) — g’'(x)/2dk. Rewriting the last equation for g(x) we obtain 


g(x) + k(2\ — k)g(x) = 0. 

We have now three possible cases: a) 2\ — k > 0, b) 2A — k = 0, andc) 2A-— k < 0. 

For case a) the general solution of the differential equation satisfied by g(x) is 

g(x) = Asin [k(2\ — k)]'?x + B cos [k(2 — k)]'z. 
Imposing the condition a) of Theorem I on g(z), i.e. kg(0) — g’(0) = 0, we find that 
B/A = [(2y — k)/k]', or 
g(x) = o,(k'” sin [k(2\ — k)]'?2 + (24 — k)'” cos [k(2X — k)]'”’2). 
For case b), when 2A = k, g’’(x) = 0, and g(x) = c,(1 + kx) where c, is an arbitrary 


constant. 
For case c) 2\ — k < 0, and 


g(x) = Asinh [k(k — 2d)]'?x + B cosh [k(k — 2a)]'’x. 
Since, however, kg(0) — g’(0) = 0, 
g(x) = c,(k'”? sinh [k(k — 2d)]'?2 + (&k — 2d)” cosh [k(k — 2d)]'”*2). 
g(x) has to be of order O(e**) with c < k, hence [k(k — 2d)]'” < k ord > 0. 

We see that the characteristic values, \, constitute a continuous spectrum, a phe- 
nomenon which is due entirely to the fact that our integral equation has an unbounded 
region of integration. For integral equations with bounded kernels, and bounded regions 
of integration it is well known that the number of characteristic values is at most de- 
numerable. If in Theorem II we set k = 1, we obtain the example treated by Lalesco’. 

By direct application of the general theory we obtain Theorem III. If g(x) is a 
solution of the integral equation 


ge) =r f eM" |x —y | oly) dy (D) 


then g(x) must satisfy the differential equation 


a) + Ko(z) = F@— AL @ + Bio 





Proof. The kernel satisfies the differential equation 
(D* — 2k?D + k‘)u(z) = 0, 


in each of the regions 0 < y < x, and0 < x < y. Writing (D) in the form 
fz) = | e*=Y(e — y)g(y) dy + | eX" (y — x)gly) dy, 
0 z 


differentiating four times and forming the expression on the left of the following equa- 
tion we obtain, 


f(a) — Qk? f(x) + k*f(x) = 2dk*g(x) + 2dg’"(z), 
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which can be written as 


(nr) — 2k f’’(x) k* f(x 
g’'(x) - k? g(x) — f{°@ —: i ) + f(x) 


If we set f(x) = g(x) in (D) we obtain the Theorem IV. If g(x) is a solution of the 
homogeneous equation 


g(x) = X [ e*'=-"! | x — y | gly) dy r>0,k>0 (FE) 
“0 


then g(x) must satisfy the differential equation 
g(x) — 2A + k’)g’"(x) — k(2N — k’)g(x) = 0. 


The real parts of the roots of the characteristic equation of the last differential 
equation have to be in absolute value less than k in order that the integral may exist 
with g(x) a solution of (E). The complete discussion of the spectrum of characteristic 
values of this equation, as well as the treatment of certain integral equations of the 
second kind is reserved for another paper. 


ON THE GENERAL EQUATIONS OF PROBLEMS OF AXIAL 
SYMMETRY IN THE THEORY OF PLASTICITY* 


By P. 8S. SYMONDS (Brown University) 


1. Introduction. The solutions of axially symmetric problems in plasticity which 
have been presented in the literature depend in general upon the assumption of the 
so-called ‘full plasticity” condition. According to this assumption, two principal stress 
differences are taken equal to the yield stress in pure tension. This implies that two of 
the principal stresses are equal, and thus reduces the problem to a statically determinate 
one. On this basis solutions have been given for a number of important problems, such 
as those of the indentation of a rigid surface of revolution into a plastic mass (see for 
example papers of Hencky’ and Ishlinsky’). 

This is a quite fictitious yield condition, with little physical or mathematical justi- 
fication. Hence it is of interest to examine the general equations for axial symmetry 
using the well confirmed Mises flow condition, and, in particular, to determine whether 
or not these equations form a system of hyperbolic type. If so, the techniques of nu- 
merical integration using the networks of characteristic curves, which have been de- 
veloped for supersonic fluid flow problems, could be applied to these problems of plas- 
ticity. Unfortunately the answer to this question, as shown in the present note, is that 
the equations are actually of elliptic type, the “characteristics” being real only under 
the special condition which reduces the problem to one of plane strain. 


*Received July 8, 1948. The results presented here were obtained in the course of research con- 
ducted under a contract sponsored jointly by the Office of Naval Research and the Bureau of Ships. 

1H. Hencky, Uber einige statisch bestimmte Fdlle des Gleichgewichts in plastischen Kérpern, Z. angew. 
Math. Mech. 3, 241-251 (1923). 

2A. Ishlinsky, The problem of plasticity with axial symmetry and Brinell’s test, Prikl. Mat. Mekh. 8, 
201-224 (1944). 
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2. The general equations. In problems of axial symmetry defined by an axial co- 
ordinate z and a radial coordinate r we are concerned with four stress components 
o,,0,,0¢,7-;, and with two velocity components u, w. The two equations of equilibrium 





are 
0c, OT;; ¢; — 
7? Ss"; | 
d OT rs 
eee I 
The Mises vield condition becomes 
(o, — o»)” + (a9 — o.)” + (6, — 0,) + 677, = 6K’, (2) 


where K is the yield stress in pure shear. The problem is evidently not statically de- 
terminate, so we must add the stress-strain relations to the foregoing. According to 
the Saint Venant-Levy-Mises flow theory* 


\ 
Ou 1 
é; = ar = 3 A(20, — Gg ~ C5), | 
u 1 | 
eo > * 3 A(2a% —o,—<4,), | 
(3) 
ow 1 
elk thes 3 N(2¢, — 6, — Gs), 





1 {ou Ow | 
~ 3 (2 + 2) Reis 


where @, , € , €s , Jrz are (tensor) strain-rate components (taken with respect to an 
arbitrary base, not necessarily time), and X(r, z) is a scalar point function depending 
on the material and on the state of stress at a generic point. We then have seven equa- 
tions for the seven unknowns (c, , 09 , 0: , Tr: , A, U, W). 

Further work is facilitated by changing notation as follows 





oo =wté, 
Oo, =~ wt 7, 
(4) 
o,=wot+f=w — €+ 7), 
Tre —_ Ty 
where w = (c, + a» + o,)/3, and &, 7, ¢ are the “stress deviation” components. In the 


following, all subscripts indicate partial differentiation. The equations become: 


%See for the instance William Prager, The stress-strain laws of the mathematical theory of plasticity— 
a survey of recent progress, J. Appl. Mech. 15, 226-234 (1948). 
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Be : 
a + é, + 7, = r (7 — §), (a) 
Equilibrium ) 
condition . 
1 ; P (5) 
;}@, — g, om a + _ iain r T; (b) 
Yield - r i 
S” 1 = “e y 
conditions ‘° Fa tints K; (¢) J 
u, = 2X, (a) 
. = rXn, ) 
Stress-strain | a (b f (6) 
relations ) ) 
relation: | «wet 6, (c) 
lu, + w, = 2dr. (d) | 


The four stress deviation components and X are determined by the five equations 
consisting of Eqs. (5) and two equations of compatibility derived from Eqs. (6). The 
compatibility equations may be written in the form 


(AE)ex = WAT) re HEL AD» | 


(7) 
r(An) zs _ 2(A7). + (AE + An); | 
as may be verified by substitution in (6). 
3. Investigation of characteristic curves. We now investigate the possibility of the 
existence of characteristic curves C, on which 


2s 2(r) (8) 


such that on C the second derivatives of the unknowns (a, &, 7, 7, \) are undefined 
although the functions themselves and all their first derivatives are defined. We first 
eliminate w = (¢, + o» + a,)/3 between Eqs. (5a) and (5b), obtaining 

] 


2-0 + Are — tre + tor : (>, ~i+=- : r). (9a) 


We next differentiate Eq. (5c) twice with respect to z, obtaining 

(2E + née + (E + n\n. + 277. = —2WEL + Em. + 1: + 72). (9b) 
The justification for this is that we are investigating possible discontinuities in second 
derivatives, hence we are interested only in how Eq. (5c) defines the relations between 
second derivatives. Finally, we write out Eqs. (7) and obtain equations on the left 
sides of which only second derivatives and the functions themselves appear: 
E,, ai NE.: + Aner + QT rs + (é + m)Arr + 2rh,: od X.. 

(9c) 
= — 2(A,€, = r.£, + Aer + ArTs + AsTr)) 


Mes + hee = —2Bdime + [2(Az). + OME + Ande]. (9d) 
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Noting that a function g(r, z) given on a curve C [defined as z = 2(r)] is a function 
of r only, on C we have 
dg dey 
= — + P. z ; ,s.* - r 
dr ¥ a ? g dr z 


where P = dz/dr. In particular if g = &, or &, we obtain, respectively: 


< — aals om Pe 

Szr dr S22) 

t a dé, _ dé, 225 
aia dr dr TP 


In like manner the second derivatives of all the functions appearing in Eqs. (9) may 
be written in terms of the double derivatives with respect to z, and the parameter 
P = dz/dr of the curve C. Equations (9) then take the form 


Aé,, + Bn. + Cte: + Dd,, = K, 


~ 


A Pes + B’.. + CO" te +. Dx. = age 


(10) 


At 


Szz 


+ B'' nz: + OC Fas + ag _ -'. 





AS + B'~,, 4. Cre. 4 Pt ids = =. 


where A, B, , D’” contain only the functions and the parameter P. The “char- 
acteristic condition” of the system is thus the equation: 








A B C - D | 
A’ B’ Cc’ D’ | 
| = 0. (11) 
| A” I tt Cc” I ww F 
AM Br qe pr 
Determining the coefficients A, B, --- , D’” for the system (9), we obtain the result: 
—2P —P —-P* +1 0 
| 2+ E+2n 2r 0 
| = 0. 
| \(P? — 1) \P? —2\P (¢ + »)P® — 27P -—¢ 
| 0 N 0 n 


This is a fourth degree algebraic equation in P, which reduces to 


[((P? — 1)(2 + ») — 47P]’ = —3[n(P* + I)’. 
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The four roots are 
a : 7) (12) 
dr §&€+ (1 + 13’) /2 ml 
where the + signs are independent; the yield equation in the original form of Eq. (5c) 
has been used in obtaining (12). 

4. Conclusion. It is clear that real characteristics do not exist in general, since P 
can be real only when 7 = o, — w = 0. But when 7 = 0 and u ¥ 0 the stress equations 
become identical with the three equations for plane strain, under the Mises yield condi- 
tion; the characteristic parameter P then reduces to 

i ve E 


fa 
S 


As may be easily verified, this is the standard result for the slope of the characteristic 
curves in problems of plane plastic flow*, which provides a useful check on the present 
results. 

An analogous but simpler investigation of Eqs. (6) leads to the conclusion that 
there are also no curves on which the velocity derivatives are undefined, even though 
the values of u, w and the stresses are defined. 


‘See, for example, William Prager, Plasticity for the aerodynamicist, J. Aero. Sci. 15, 253-262 (1948). 


EFFECT OF HYPERBOLIC NOTCHES ON THE STRESS DISTRIBUTION 
IN A WOOD PLATE* 


By C. BASSEL SMITH, University of Florida 


1. Introduction. If from an orthotropic material a flat plate is cut parallel to a 
plane of elastic symmetry, it will have two perpendicular axes of symmetry in the 
plane of the plate. Such a plate is said to be orthotropic." An example of this kind of 
plate is a plain-sawn board. In the discussion that follows, an orthotropic plate is 
assumed to be bounded on two sides by hyperbolic notches given by the equation 
y'/a — x°/b’ = 1, and indefinitely extended in the other direction. The z- and y-axes 
are taken parallel to the axes of symmetry. The plate is subjected only to forces directed 
parallel to the x-axis and applied in the plane of the plate. These forces are assumed to 
act at great distances from the y-axis, and in such a way that the traction over any 
cross-section perpendicular to the z-axis is statically equivalent to a single force of 
magnitude P directed along the x-axis. The problem will be treated as one of plane 
stress. 

2. The stress distribution in a wood plate with hyperbolic notches. The components 
of stress and strain in the orthotropic plate described are connected by the following 


relations’: 


*Received June 11, 1948. 

1A, T. Price, Phil. Trans. (A) 228, 1-62 (1928). H. W. March, Physics, 7, 32-41 (1936). U. S. Forest 
Products Laboratory Reports Nos. 1300, 1301, 1304, 1312, 1316. 

*U. S. Forest Products Laboratory Report No. 1503. Love’s notations for the stress and strain 
components are used, (A. E. H. Love, The mathematical theory of elasticity, 4th ed., Cambridge University 
Press, Cambridge, 1927). 
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= Pea _ Guz . 


re E, 
as Ss ply (1) 
a aa a 
ee 
Mery 


In these equations 2, and EZ, are Young’s moduli in the z- and y-directions, respectively. 
Poisson’s ratio o,, is the ratio of the contraction parallel to the y-axis to the extension 
parallel to the z-axis associated with a tension parallel to the x-axis and similarly for 
o,,. The quantity y,, is the modulus of rigidity associated with the directions of x and y. 
Since the equations of equilibrium are deduced with no reference to the law con- 
necting stresses and strains, they will be satisfied by a stress function F such that 


a°F : oF : a°F 
—; = —; y mt + : 2 
oy’ . Ox’ x, Ox OY (2) 





Substituting Eq. (2) into Eq. (1) and then making use of the compatibility equation 


Oe, , Vly Dery 
ay * ax" ~ ax ay’ 
we obtain the following differential equation, satisfied by the stress function F*: 
io atm in 
a" 4+ 2K a, (a = 0 (3) 
where® 
oF ws ee = 72:2) (4) 
and 
n= ey (5) 
with 
e = (E£,/E,)"™. (6) 
The solution of Eq. (3) is facilitated by taking F to be of the form 
F = R{F\(x + ian) + F(x + i8n)) (7) 
where 
= [K + (K’ — 1)”, 6 = [K — (K’ — 1)""]"", (8) 


i is the imaginary unit, and the letter R means that the real part of the expression in 
the bracket is to be taken. The functions F; and F, are any analytic functions of the 
complex variables x + ian and x + 78», respectively. 


‘ 
*For an orthotropic plate like the one discussed here the relation o:,Zy = ¢yzHz holds. See U. 8 


Forest Products Laboratory Report No. 1503, p. 7. 
‘For wood, K as defined here is probably always greater than 1. 
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The stress distribution in an orthotropic plate in a state of plane stress is then ob- 
tained by finding a solution to Eq. (3) in the form of Eq. (7), and in addition requiring 
that the function F satisfies the following relations on the boundary of the plate’: 


d | rc d (2) — 
ds (3 wee ds\ dx] — es 


/ 


where ds denotes an element of the boundary, and X, ds and Y, ds represent the z- 

and y-components of the force per unit thickness of the plate acting on the element of 

arc ds from the positive side as indicated by the direction of v, the normal to the boundary. 
A suitable choice of F for the problem considered is 





F = R\pBiz, log (2; + w,) — w, — 22 log (22 + w.) + wo} (9 
where 
4-24 1an, Z=xat 1Bn, 
lp 2\1/2 2 2\1/2 
WwW, = (4+ yi)’ 7» W. = (Zo + Yo)’, 
1=b+a'ea’, v2 = b° + Bea’, 
a a 
B= —— —————————— — =>, (10) 
€{@ arc tan (eX) — a are tan (aed)} b 
and p = P/2ad, the average stress on the cross-section x = 0, since the plate is taken 


to be of thickness d. To make the functions appearing in Eq. (9) single valued the fol- 
lowing relations must be assumed always to hold: 


R{w,] > 0, R{w.] > 0, 


T T Tv T 
— 5 < amp (@, + w) < 5, — 5 < amp (% + mw) <5 (11) 


a ~- — 


By means of Eqs. (2) and (9) the stress components are found to be 


a°F —— ; ia L) 2) 
Ox" i= r| pa( 2 We | ve 


oF | a r| pb(=2"¢ 4 ee) | (13) 





ay” - 1 We 
oF : ( —1a€ ie) | 
_ — . =— as oe —— ° l 4 
Ox Oy x, r| pe W We 4) 


The stress component X, as given by Eq. (13) was calculated at several points 
along the y-axis for a plain-sawn plate of Sitka spruce with hyperbolic notches, as in 





58. Timoshenko, Theory of elasticity, McGraw-Hill Co., New York, 1934, p. 179. U. S. Forest 
Products Laboratory Report No. 1510, p. 7. 





oe 
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the problem discussed here. The x-axis was taken parallel to the grain of the wood. 
The elastic constants were taken to be: 


| 


= 1.679 X 10° lb./in.’, uz, = 0.112 X 10° Ib./in.’, 


E, = 0.076 X 10° lb./in.’, o,, = 0.464 


The results of these calculations for several values of \ are given in the table below. 





VALUES OF X; FROM Ea. (13) 





= 10 A= 1.0 A = 0.1 
X-/p | y | X:/p y X./p 

— -—— — - —-— — - [/ — —E 
0.0a 0.6373 0.0a | 0.6989 0.0a 0.9784 
0.2a 0.6505 0.2a 0.7142 0.2a 0.9808 
0.4a 0.6955 0. 4a 0.7662 0.4a 0.9884 
0. 6a 0.7971 0. 6a 0.8832 0.6a 1.001 
0.8a 1.064 0. 8a 1.180 0.8a 1.020 
1.0a 30.84 1.0a 3.192 1.0a 1.046 


It is evident that the value of the stress component X, at x = 0, y = a is the maxi- 
mum stress occurring in the plate in each case. Setting x = 0, y = a in Eq. (13), it is 


found that 
= ——_—— 8, (15) 
b 


which gives the ratio of the maximum stress to p, the average stress on the cross-section 
a 0. 

3. Isotropic stress function as the limiting case of the orthotropic stress function. 
It is of interest to obtain from Eq. (9) the stress function for the corresponding problem 
for an isotropic material. For an isotropic material a = 8 = e = 1, as is shown by 
Eqs. (4), (6), and (8). For these values, the right side of Eq. (9) becomes indeterminate. 


Evaluating, it yields 


s (1 + \)ap atl ‘* — ili | 
Alo +™MyarctanAt+aAle+u log (@ + a | (16) 


where 


z=2+ wy, 
(17) 
w= (+7), y =a’ + bd’. 


To make F single-valued, relations similar to those of Eq. (11) must hold. By letting 


z= yz’ v(x’ + ty’) in Eq. (16), it is found that 
; ?F 1 OF ° (F 
Y, = : P => 357 = 5 (4) (18) 
Ox y Ox Ox” \y 
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with similar expressions for the other stress components. If, therefore, in Eq. (16) we 
replace z by yz and then divide the right side by y’, the resulting expression can be used 
as the stress function. Hence the change to elliptic coordinates, which is usually done 
by replacing z by y sinh (u + iv), can be accomplished in Eq. (16) by substituting sinh 
(u + iv) for z and then dividing the right side by y*. Thus we find the following expres- 
sion for the stress function in elliptic coordinates: 

SiN Vo 


F ————_"—— [wy + (sin’ v,)e™“ cosv (19) 
Vp + SiN Up COS Vo ly . ] 


where the hyperbola y*/a’ — x°/b*” = 1 is given by the particular value, vp , of the co- 
ordinate v. Equation (19) is the stress function for the corresponding problem involving 
the isotropic plate.° 


6See paper by H. Neuber, Zeitsch. f. Angew. Math. u. Mech. 13, 439 (1933). 


BOOK REVIEWS 


Eigenfunction expansions associated with second-order differential equations. By E. C. 
Titchmarsh. Oxford, at the Clarendon Press, 1946. 175 pp. $7.00. 


The subject matter of this book is composed of the analysis required to develop and justify the expan- 
sion of an arbitrary function f(z) in the eigenfunctions of the differential equation y"’ + [A — q(x)]y = 0. 
The classical theory is disposed of in Chapter 1 and the remainder of the book is largely devoted to those 
situations where q(x) is singular at an end-point of the region and/or those where f(r) is defined over an 
infinite or semi-infinite domain. In Chapter 4 one finds as examples the expansion formulae for Hermite, 
Bessel, Legendre, Sonine, Laguerre, and hypergeometric functions. The latter chapters deal with the dis- 
tribution of the eigenvalues, the nature of the spectrum, and with questions of convergence and summa- 
bility. 

The book is intended primarily for the mathematician, but its interest and value to the physicist 


and engineer should be great. 
G. F. CARRIER 


Modern operational calculus with applications in technical mathematics. By N. W. 
McLachlan. MacMillan and Co., Ltd., London, 1948. xiv + 218 pp. $5.00. 


This introduction to the theory of the Laplace Transform is intended as a guide to those interested 
in the application of this transform to engineering problems. However, the questions of rigor are treated 
in more detail than is customary in texts with this purpose. In fact, nearly one third of the book consists 
of appendices dealing with such items as convergence questions. In Chapters 1 and 2 the transform is 
defined and the fundamental theorems and identities are deduced. In Chapters 3 and 4, the solution by 
transform methods of ordinary and partial linear differential equations is discussed. These and the subse- 
quent chapters which, for the most part, are concerned with integral evaluation contain many examples 
drawn largely from electric and acoustic problems. It is regrettable that no reference has been made to the 
role of this transform in solving stability problems. 

The book should be especially useful to engineers whose background is essentially non-mathematical 
since the points of rigor discussed appear to have been selected for just such an audience. 


G. F. CARRIER 









































